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An Experimental Multichannel Pulse Code 
Modulation System of Toll Quality 


By L. A. MEACHAM and E. PETERSON 


Pulse Code Modulation offers attractive possibilities for multiplex telephony 
via such media as the microwave radio relay. The various problems involved 
in its use have been explored in terms of a 96-channel system designed to meet 
the transmission requirements commonly imposed upon commercial toll circuits 
Twenty-four of the 96 channels have been fully equipped in an experimental 
model of the system. Coding and decoding devices are described, along with 
other circuit details. The coder is based upon a new electron beam tube, and 
is characterized by speed and simplicity as well as accuracy of coding. Thes¢ 
qualities are matched in the decoder, which employs pulse excitation of a simple 
reactive network. 

I. INTRODUCTION 

N THE development of systems for transmitting telephonic speech, much 
effort has been directed toward minimizing the effects of noise picked up in 
the transmission medium. ‘The system described in this paper represents one 
method which has been successful in eliminating completely such effects under 
appropriate and practical conditions. In this method, known as Pulse Code 
Modulation'** (PCM), telephone waves are represented by sequences of 

on-off constant-amplitude pulses. 

Perfect reception of such pulses demands simply recognition of whether any 
pulse exists or not. Recognition can be carried out effectively in the presence 
of noise and interference amounting to a substantial fraction of the pulse 
amplitude. In contrast, telephonic waves carry information by subtle ampli- 
tude variations in the course of time. High quality telephone reception 
accordingly demands a much lower ratio of noise and interference—lower by 
as much as 50 decibels. 

The magnitude of this figure exhibits one good reason for exploring the 
possibilities of PCM. Another potent reason, which is of particular impor- 
tance in multi-link transmission, is that, with pulses involving just two standard 
values, regeneration can be used at repeater points and at the receiver to wipe 
out transmission impairments. The regeneration process consists of the pro- 

1A. H. Reeves, U. S. Patent 2,272,070. 

“Telephony by Pulse Code Modulation”, W. M. Goodall, Bell System Technical 
Journal, July, 1947. 

3 “Pulse Code Modulation”, H. S. Black and J. O. Edson; presented June 11, 1947 at 


the Montreal Summer Meeting of the American Institute of Electrical Engineers, and to 
be published in the A. J. E. E. Transactions. 
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duction of a properly formed standard pulse, free of noise, to correspond with 
each received pulse, even though the latter may be considerably misformed. 
The sole proviso here is that before regeneration the level of noise and distortion 
in each link be kept below the comparatively large threshold value at which a 
mark cannot be distinguished from a space. If this holds good throughout 
the transmission path then literally the received pulses can be made as good as 
new. In contrast it is impossible fully to repair or to regenerate signals not 
involving standard values of amplitude and of time. With such signals dis- 
tortion and noise in each span contribute to the total which therefore increases 
with the system length. 

To sum up, conversion of speech to a code of pulses and spaces permits tele- 
phony to assume certain new and desirable properties; ability to work with 
small signal-to-noise ratios, and ability to regenerate any number of times with 
no degradation of quality. These advantages do not accrue without certain 
penalties. Conversion of speech waves to pulse form and back to speech 
involves a certain degree of apparatus complexity at the terminals. This 
complexity is not decreased by the need to handle pulses at high speeds, of the 
order of a million per second. Here radar and television circuit techniques 
are helpful. Another characteristic is that a greater band width is occupied 
in the transmission medium. ‘This arises through the operation of two factors, 
of which one is the use of double sideband in pulse transmission (as against 
single sideband in carrier telephony), and the second invelves the number of 
pulses used in the code. The relatively wide band required can best be accom- 
modated in the microwave region and it happens that the properties of on-off 
pulse transmission can be used there to particular advantage. 

The PCM system to be described was set up to evaluate experimentally 
the problems involved in providing multichannel facilities of toll system 
quality. It was designed to accommodate 96 channels. For experimental 
studies of such things as crosstalk and methods of multiplexing channels, a 
fraction of the total number of channels is sufficient and only 24 of the 96 were 
built. These are arranged as two groups of 12 channels each. The channels 
of a group are assembled on a time division basis. Assembly of the groups is 
carried out on a frequency division basis, each group amplitude-modulating 
its own carrier. In a planned alternative arrangement the group pulses may 
be narrowed and interlaced to put all 96 channels in time division on a single 
carrier, but this alternative will not be explored here. 

The assignment of 12 channels per group fits in well with the present arrange- 
ment of carrier telephone circuits used in the Bell System plant, such as Types 
J, K,and L4 Use of time division for a group of this size involves pulses with 

4A Twelve-Channel Carrier Telephone System for Open-Wire Lines,’”’ B. W. Kendal! 


and H, A. Affel, Bell System Technical Journal, January, 1939. ‘Coaxial Systems in the 
United States,” M. E. Strieby, Signals, January-February, 1947. 
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repetition rates up to 672 kilocycles, and pulse durations as short as a quarter 
microsecond. These pulses, obtainable from more or less standard types of 
vacuum tube circuits, can be distributed from point to point in the equipment 
without too much difficulty. Amplitude modulators and demodulators at two 
neighboring carriers—65 and 66.5 megacycles—then serve to bring the PCM 
signals into the intermediate frequency range for transmission to and from the 
microwave equipment. 

The speech quality of the overall system in respect to such factors as band 
width, volume range, noise, distortion, and crosstalk more than meets the 
requirements generally imposed upon such systems. 

Figure 1 is a front view photograph of the experimental apparatus setup with 
covers removed from one bay to show typical construction. The two end bays 
contain intermediate frequency modulators and demodulators required for the 
two groups. In addition voice frequency terminating sets and jacks are 
mounted here, together with testing equipment. The center bay and the one 
to the right of it are identical; each includes apparatus for handling a group of 
twelve message channels. Transmitting equipment is mounted in the upper 
half, and receiving equipment in the lower half of each bay. The remaining 
bay, second from the left, holds all the timing equipment needed to furnish 
control pulses for operating eight of the message bays, a total of 96 channels. 
Included are circuits for synchronizing the receiver. Individual regulated 
power supplies are mounted near their loads on the several bays. 

Figure 2 is a rear view of the same equipment. Cables in the four horizontal 
ducts shown carry control pulses from the timing bay to the 12-channel group 
bays. These ducts are large enough in cross-section to handle all the cables 


required for a complete 96-channel terminal. 


II. FUNCTIONAL PROBLEMS INVOLVED 


The broad problems brought together in building this system may be con- 
sidered under the four classes following: 
1. The pulse code modulation problem; to convert signal waves to pulse pat- 
terns. 
2. The multiplex problem; to aggregate channels into groups and groups into 
a supergroup. 
3. The transmission problem; to fit the system into the minimum required 
band width, and to remove the effects of transmission impairments. 
4. The pulse code demodulation problem; to convert pulse patterns back 
to the original signal waves. 
These are to be discussed from a functional standpoint to provide background 
for discussion of the specific equipment. 
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Fig. 1. Front view of experimental PCM terminal equipment, with covers removed 
from a 12-channel group bay. 


Pulse Code Modulation 


Sampling. A basic premise of pulse mcdulation systems is that the infor 
mation content of a wave can be conveyed by samples taken at sufficiently 
frequent, equally spaced time intervals. The interval should be no greater 
than half the period of the highest frequency speech component to be repro: 
duced or, otherwise expressed, the sampling rate should be not less than twice 
the frequency of the highest speech component present. This provision insures 
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Fig. 2. Rear view. Cables in the horizontal ducts carry timing pulses leftward from 
the synchronizing bay. 


that sidebands produced by sampling do not overlap to introduce distortion, 
as illustrated in Fig. 3. For a speech band extending up to 340) cycles, a 
reasonable sampling rate is eight kilocycles. 

Samples may be intermittently-transmitted portions of the signal wave, of 
appreciable duration, or they may be essentially instantaneous amplitudes. 
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To afford time for coding, the instantaneous samples may be maintained at 
constant value for an appropriate interval—a process here referred to as 
“holding.” 

Quantization. The fundamental operation of PCM is the conversion of a 
signal sample into a code combination of on-off pulses. In any practical 
system a continuous range of signal values cannot be reproduced since only a 
finite number of combinations can be made available. Each combination 
stands for a specific value, of course, so that we wind up by representing a con- 
tinuous range of amplitudes by a finite number of discrete steps. This process 
is spoken of as quantization, a quantum being the difference between two 
adjacent discrete values. Graphically this means that a straight line repre- 
senting the relation between input and output samples in a linear continuous 
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Fig. 3. Spectrum of a sampled audio band, illustrating separation of components when 
the sampling frequency is at least twice the top audio frequency. 


system is here replaced by a flight of steps as in Fig. 4a. The midpoints of the 
treads fall on the straight line, and the height of the step is the quantum. 

Manifestly the greatest error inherent in quantization amounts to half a step. 
Hence the quality of reproduction may be measured by the size of that inter- 
val, which depends upon the total number of steps in the amplitude range 
covered. With ” pulses assigned to represent an amplitude range, the maxi- 
mum number of discrete steps is 2", and the size of each step is proportional to 
2—*" times the amplitude range. 

This error shows up as a noiselike form of distortion, affecting background 
noise in the absence of speech, and accompanying speech as well. The dis- 
tortion actually consists of a multiplicity of harmonics and high order modu- 
lation products between signal components and the sampling frequency scat- 
tered fairly evenly over the audio spectrum. If the audio signal is a simple 
sine wave, these many products may be identified individually ; but for speech 
or other complex signals they merge into an essentially flat band of noise that 
sounds much like thermal noise. Since the level of this distortion is fixed by 
the quantum size, an adequate number of steps must be provided for the lowest 
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MULTICHANNEL MODULATION SYSTEM 7 


amplitude sounds it is necessary to transmit. Considering that consonant 
levels may be of the order of 30 decibels below vowels, and that weak talkers 
may be of the order of 30 decibels down from loud talkers, it is clear that ampli 
tudes as far below maximum as 60 decibels require at least a few steps 
Ordinarily this would involve a large number of pulses for transmission, 


with a consequent increased complexity of terminal apparatus, and ar 


Fig. 4. Relation between input and quantized output, with quantization uniform in 
a) and tapered in (b). 


creased frequency band per channel. A more practical solution is to employ 
tapered steps rather than uniform ones. In this way a given number of steps 
can be assigned in greater proportion to the low amplitudes than to the highs, 
as shown in Fig. 4b. There results a degree of step subdivision sufficient to 
care adequately for the low-amplitude sounds including background noise 


A small penalty is paid at the upper end of the amplitude scale because of the 


proportionately smaller number of steps available there, bringing in higher 
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quantizing noise in this range. How large the effect may be is determined b) 
the degree of taper. To apportion the noise at different levels, the distributio 
of steps over the amplitude ragge can be varied. 

Preliminary studies of quantization, involving listening tests and _nois: 
measurements for various numbers of digits and various kinds of taper, led to 
the choice of a seven-digit code (128 steps) for the present system. The taper 
employed reduces the smallest steps 26 decibels telow the average size and 
increases the largest ones akout 6 decibels. 

Coding. ‘The coder is required to set up a pulse code combination for eac! 
quantized signal value. A great many codes are conceivable, but in practice 
a simple one in which the pulses correspond to digits of the binary number 
system allows greatest simplicity at the receiver. 

While coders may take a wide variety of forms, they can ke arranged in 
three categories according to the way in which they evaluate speech amplitudes 
In the first category an amplitude is measured by counting out, with a binary 
counter for example, the number of units contained in it one by one until the 
residue amounts to less than a unit. In the second, the amplitude is measured 
by comparison with one digit value after another, proceeding from the most 
significant digit to the least, and subtracting the digit value in question each 
time that value is found to be smaller than the amplitude (or its residue from 
the previous subtraction). In the third, amplitude is measured in toto by 
comparison with a set of scaled values. Modulators disclosed by Reeves' 
and by Black and Edson? are of the first category. That described by Goodall’ 
belongs to the second, and the one described in the present paper is in the third 
category. Generally speaking the number of operations and the time required 
for coding decrease in going from the first to the third. Rapid coding is obvi- 
ously desirable since it allows more channels to be handled in time division by 


common equipment. 


Multiplex 

Channels may be multiplexed by arranging them in time sequence, or by 
arranging them along the frequency scale. ‘These methods are known as time 
division and as frequency division, respectively. The first is accomplished 
by gating, or switching, at precisely fixed times. One way of doing this im- 
presses a more or less rectangular pulse on one of the grids of a gate tube, so 
that a signal wave may be transmitted during the gating pulse. ‘The second 
method here refers to the use of amplitude modulators, each supplied with an 
appropriate carrier, which translate the signals to their assigned positions on 
the frequency scale. To avoid crosstalk in a time-division system, operations 
in group equipment common to a number of channels must proceed without 
memory of the amplitudes of preceding channels, requiring build up and decay 
times to be held within limits. This implies a sufficiently wide pass band 
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together with phase linearity. For the same purpose in a frequency-division 
system, filters are used to select and to combine channels. To avoid crosstalk 
the filters must be sufhciently selective, and amplitude non-linearity must be 
held within limits. 

In the present system, the pulse code is delivered by the coder as on-off 
pulses in time sequence. It is therefore natural to organize the pulses of the 
different channels so that they appear in sequence, thus forming a time-division 
multiplex. Most types of coder require an appreciable length of time, after 
delivering the pulses of one channel, to prepare for coding the next. This 
preparation time may be afforded conveniently, without introducing gaps in 
the pulse train between assignments of consecutive channels, by providirg two 
coders, which take turns at the channe'!s in each time-division group. 

As the number of channels increases, evidently the time interval which can 
be assigned to each channel must be reduced since all of them must be fitted 
into one period of an 8-kilocycle wave. Similarly the allowable duration of a 
code or digit pulse becomes shorter as the number of time-division channels 
in a group is increased. Then too, pulses tend to become more difficult to 
generate and transmit as their duration decreases. For these reasons it is 
desirable, and eventually it becomes necessary, to restrict the number of chan- 
nels included within a time-division group. Frequency division may then be 
used to aggregate several time-division groups. 

For our purposes groups of 12 channels are multiplexed by time division. 
With seven digits per channel, each group has a capacity of 672,000 pulses per 
second. To combine eight of the groups for a 96-channel system we again have 
a choice between frequency division and time division. The equipment of 
Fig. 1 is laid out to accommodate either procedure. In the first case each 
group is assigned a carrier for amplitude modulation, as used in actual tests 
to be described. For the second case the pulse durations would be cut down 
by a factor of eight, and the pulses from the different groups interlaced. Here 
the supergroup would have a capacity of 5,376,0C0 pulses per second. 

In carrying out coding operations, and in multiplexing on a time-division 
basis, various control pulses are required which differ in repetition frequency, 
in time of occurrence, and in duration. These may be generated from a stable 
base frequency oscillator through the use of harmonics or, alternatively, of 
sub-harmonics. In a time-division system which requires a variety of flat- 
topped waves for switching operations, the use of sub-harmonics fits naturally. 
Frequency step-down circuits of the multivibrator type produce waves either 
approximating the desired forms directly, or requiring only simple circuits for 
reshaping. In contrast harmonic generation requires filters for component 
wave selection, more and more elaborate in structure as the order of the wanted 
harmonic goes up. Then, after selection, the harmonic has to be amplified 
and limited or otherwise shaped for switching purposes. Generally speaking 
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we need less apparatus and less space if we use multivibrator step-downs 
Another advantage is that noise in the base frequency supply produces propor- 
tionately less phase jitter with sub-harmonics. While multivibrators do not 
have high inherent stability, they are capable of great precision when suitably 
controlled. For these reasons frequency step-downs are used to generate ail 
the timing waves of the system. 

Timing and gating operations require accurate time alignment of pulses 
which may be accomplished by suitably delaying one set with respect to the 
other. For this purpose use is made of delay networks or cables, or delay 
multivibrators. Pulse durations may be varied through the use of interference 
effects between given pulses and their delayed replicas. Additional timing 
and gating wave forms are required for regeneration and for assembly in time 
division multiplex. All such control pulses are economically generated at a 
single common point rather than by a number of local generators, and can then 
be supplied to the message equipment by common power amplifiers via shielded 
cable. 


Transmission 


In this section we are to consider the genera! factors entering into satisfactory 
reception of on-off pulses including such limitations as those on band width 
and noise. These are to be viewed while keeping in mind the procedures 
available for pulse regeneration. 

lf we start with a rectangular pulse like one of those generated at the trans 
mitter, the corresponding frequency spectrum exhibits lobes extending indefi 
nitely on the frequency scale, with progressively decreasing amplitudes as 
indicated in Fig. 5a. When such a pulse is passed through a linear phase 
filter® which discriminates against frequencies beyond the first lobe, the re- 
sulting pulse is practically of the sinusoidal form shown in Fig. 5b. Reducing 
the high-frequency response to the extent prescribed rounds the corners of the 
transmitted pulse, its duration at half value remaining equal to that of the 
original input pulse. In practice, transmission characteristics depart from 
phase linearity and low frequency cutoffs exist. Both effects introduce irreg- 
ularities into the pulse response. While actual pulses therefore differ slightly 
in detail from the idealized picture given above, that picture will be retained 
for simplicity in discussion. 

The band width needed for good pulse transmission can be minimized by 
making pulses as wide as possible. But since a specified number of pulses have 
to be put into a given time interval (84 pulses in go'yq sec.) we must limit our 
broadening at a value where the presence or absence of a pulse may be clearly 

§ A suitable filter is one with a Gaussian characteristic, the loss in decibels varying as 


the square of frequency and having a value of 1 neper (8.68 decibels) at a frequency equa! 
to 1/" 
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determined in the presence of noise and interference. This spaces the sinus 
oidal pulses so that they overlap at their half-value points as illustrated in 
Fig. 5c. There it will be observed that no matter how many pulses occur in 
succession, the maximum amplitude of the pulse train is no different from that 
of a single pulse. The spectra of all such pulse trains have a common envelope, 
shown dashed in Fig. 5c. 
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(Cc) 
Fig. 5. Pulse forms and their associated amplitude spectra for: a, rectangular pulse; 
b, single lobe of a sinusoidally varying pulse; c, succession of pulses, each like that in (b) 


Further, the band needed to transmit any sequence of such pulses is the 
same as that needed for a single pulse. This can readily be shown by using 
the familiar relation between transient build-up time and band width. More 
over, this conclusion is consistent with statistical analysis of all possible pulse 
combinations, which yields a spectrum of the form of Fig. 5b. 

The relation between pulse duration and band width here described gives 
close to the optimum ratio of signal to noise and interference for the system 
considered. Narrowing the band would increase build-up and decay times, 
leading to reduced pulse amplitude and to increased pulse overlap. Thus, 
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although the narrower band would admit less extraneous noise, margins over 
noise and interference would be reduced. Widening the band, on the other 
hand, would allow pulses to build up and decay faster, but would not increase 
the pulse height. Thus the same signal would result, but the widened band 
would pass increased noise, and again the margins would be reduced. The 
optimum band width represents the most useful compromise in efforis to re- 
duce noise and interference and to increase signal. 

The filter characteristic we have been discussing is that of the entire link 
taking in all selectivity inserted between the practically rectangular pulses 
originally generated at the transmitter and the pulses delivered to the PCM 
receiver. Filters at both transmitting and receiving terminals of the link are 
included, the greater part of the overall selectivity being located at the receiver. 
With about 1.5 microsecond available per pulse at half amplitude, filters spaced 
1.5 megacycle apart accommodate the double sideband and keep the inter- 
ference between groups within tolerable limits. 

To establish the presence of pulses we can set up an amplitude threshold equa! 
to half the normal pulse amplitude, and test to see if that threshold is exceeded 
at a time near the center of the assigned pulse position. Selecting the threshold 
at half amplitude provides equal margin against the possibility of noise and 
interference bringing the full pulse amplitude, or mark, below threshold and 
bringing the nominal space above threshold. ‘Testing at the pulse position 
midpoint maximizes this margin. 

The amplitude threshold is set up by slicing a thin section horizontally out 
of the pulse at its half-amplitude level by means of an amplitude discriminator 
Evidently, this procedure restores the flat top of the pulse. To complete re- 
generation by restoring the pulse epoch to a standard value the sliced pulses 
are gated at the midpoints of their proper intervals with narrow pulses supplied 
by the timing equipment. By these two pulse regenerating processes—slicing 
and gating—noise and interference are made impotent to produce errors until 
they attain a substantial fraction of the pulse amplitude. With the effects 
of noise and distortion thus eliminated, the only noise inherent in the system 
is that of quantization. In long systems having many repeater points, re- 
generation has to be practiced at spans short enough to permit cleaning out 
noise and distortion before it piles up above threshold. ‘Thus in this system 
transmission impairments have to be considered only for the span between 
regeneration points; with their effects limited below threshold they are not 
carried over from one span to the next. 

Where PCM groups are multiplexed by frequency division, amplitude non- 
linearity of the system must be kept within limits. Otherwise intermodulation 
products may fall within transmission bands, adding to interference. Over- 
lapping of neighboring filter bands also must be kept within bounds to reduce 
direct crosstalk between the pulse trains. With pulses arranged exclusively 
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in time-division multiplex, however, amplitude non-linearity of itself is not a 
factor. In this case the limitation comes on the departure from a suitable 
attenuation characteristic and from constant delay with respect to frequency. 
Distortion from these sources may increase the pulse duration so that excessive 
overlap of adjacent pulses will leave less margin available for interference and 


noise. 
Pulse Code Demodulation 


At the receiver, the regenerated code pulses are operated upon to recreate 
as closely as possible the original signal by operations complementary to those 
at the transmitter. Unfortunately, however, no process at the receiver can 
undo the effects of quantization which remain as noise, so that the quanta must 
obviously be made small enough from the beginning to permit satisfactory 
speech quality. 

In each time-division group alternately working decoders may be used, in 


| the same way as the two coders at the transmitter. Routing is effected by 


} suitably timed gates. Conversion of a pulse code to amplitude may be ac- 


aie 


complished by causing each pulse of a code combination to contribute an ampli- 


tude corresponding to the binary digit it represents, and then summing the 


} contributions. When tapered steps are employed, due consideration must be 


4 


Ree SE 





given to overall linearity, discussed subsequently. The resulting output is a 
pulse-amplitude-modulated signal, which is then distributed to the channels of 
the group by an electronic commutatar. Reconstruction of the signal from 
these distributed pulses is accomplished by simple filtering, which serves to 
remove components extraneous to speech introduced by the sampling procedure 
and tied up with the sampling rate. 

Production of the necessary timing pulses at the receiving end proceeds in 
much the same manner as at the transmitter except for one thing. That is, 
instead of being initiated by a local oscillator, the receiver timing must be linked 
to the input pulses so that they may be properly routed. This invo'ves the 
problem of synchronizing or, to borrow a term from television, framing. Use 
of this term is based upon the similarity of the sequence of PCM digits within 
a single sampling period to the complete ordered array of television picture 
elements. Preferably framing should be done with a minimum of time interval 
and of band width. 

One method of synchronizing pulse systems employs a marker pulse which 
serves to initiate a timing sequence for each frame at the receiver. Here the 
marker pulse must differ sufficiently from the other pulses to permit its rapid 
and certain identification. This is ordinarily done by making the marker 


; several times as long as any message pulse. In PCM, however, where digit 


pulses are run together in many code combinations, the marker would have 
to be very long to be clearly distinguishable, thereby cutting into channel ca- 
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pacity. This inefficiency can be avoided by deriving the pulsing rate from the 
pulse train through the use of a narrow band filter. Framing may then be 
effected by using in each frame just one digit pulse, which is given a distinctive 
repetition rate. With this method, a certain amount of time is required to 
establish synchronism when the system is started up. In the system to be 
described the framing time is less than one-tenth second—a tolerable value. 


II]. EXPERIMENTAL SYSTEM 


In the light of the foregoing discussion, the block diagram of the experimental 
system shown in Fig. 6 is believed to be largely self-explanatory. It will be 
noted that for microwave transmission the modulated intermediate-fre- 
quency signals are simply translated in frequency to the 4000-megacycle 
band. The shaping filters and group selection filters shown have approxi- 
mately Gaussian characteristics, in accord with transmission considerations 
noted earlier. The band widths shown for these filters apply between points 
one neper down from the midband loss. Band widths given for the ampli- 
fiers, on the other hand, refer to their regions of essentially flat response. 

Overall measurements are facilitated and the amount of experimental 
apparatus minimized by looping the radio path through a non-regenerative 
microwave repeater 21 miles away in New York. Both ends of the 24 
channels are thus made available at Murray Hill, New Jersey, the location 
of the apparatus pictured in Fig. 1. With this arrangement, and using 
conventional 4-wire voice frequency terminating sets, 24 people are able to 
engage in 12 simultaneous conversations through the system. 

PCM Transmitter. The transmitting equipment for an individual 12- 
channel group is shown schematically in Fig. 7. Each audio input* is passed 
through a 3400-cycle low-pass filter and through a limiter which chops off 
the positive and negative peaks of any signal exceeding a.prescribed maxi- 
mum amplitude. This limit is chosen to penalize the loudest talkers to the 
degree customary in toll system practice. The inputs then enter a ‘‘col- 
lector” circuit, which assembles samples of the channels in time division 
multiplex on a common lead. Although it functions electronically under 
control of pulses from the timing bay, the circuit so resembles a mechanical 
commutator that this analogy has been used in the schematic. The period 
of rotation of the “‘contact arm” is 125 microseconds (8 kilocycles), and a 
conducting path is formed to the common multiplex lead from each channel 
circuit in turn for y's; of this period, or 104°; microseconds. It should be 

6In telephone terminology, these 4.wire inputs are normally at the —13 decibel Jevel 
point; i.e., 13 decibels below the transmitting level at the toll test board. Strap connec- 
tions are provided, as in the Western Electric A-2 channel bank, for adapting the system 
to inputs 3 decibels smaller. Similarly the final 4-wire outputs are delivered at the +4 


(or +7) decibel point. The normal gain through a link is thus 17 decibels but can be 
set as much as 6 decibels greater. 
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18 MULTICHANNEL MODULATION SYSTEM 


noted that the channel samples thus collected are not “held” at this stage; 
i.e., each sample does not remain constant in potential during its assigned 
interval, but rather changes to follow the wave form of the audio signal in 
the corresponding channel. 

The multiplex signa] is supplied to an instantaneous compressor, which 
employs silicon rectifier elements to give an input-output characteristic of 
the general form indicated within its block (Fig. 7). To understand the 
purpose of this device, we must recall the discussion of quantizing noise 
given earlier. There it was found desirable to provide a tapered distribution 
of step heights in the staircase-like quantizing characteristic, thus devoting 
a considerable number of small steps to the treatment of background noise 
and low-level signals. Although coders have been devised which inherently 
deal with signal amplitudes in this graded manner, it has been found more 
practicable in the present system to apply amplitude compression to the 
samples before coding, and to divide this compressed amplitude range into 
uniform steps in the coder. The result is a tapered step distribution with 
respect to the original uncompressed scale of amplitudes, details of the dis- 
tribution being determined by the shape of the compression characteristic. 

It may be well to note here that this method can be used in reverse at the 
receiver, with the decoding performed on an equal-step basis, and the re- 
sultant samples passed through a complementary instantaneous expandor. 
If the compression and expansion are truly complementary, the overall 
characteristic relating amplitudes of input and output samples will be linear 
except for the tapered array of quantizing steps (Fig. 4b). 

Incidentally, no added band width in the transmission path of this system 
is required to accommodate the instantaneous compandor action 

After compression, the multiplex signal is delivered at low impedance to 
the inputs of two coders. In Fig. 7 the switch analogy is called upon again 
to illustrate the routing of alternate samples to the ‘“‘odd” and ‘“‘even”’ coders, 
and concurrently the storage of these samples on “‘holding capacitors” to 
keep them unchanged during the coding operation. Here the switches 
rotate at 48,000 revolutions per second, each one closing six times in a com- 
plete 8-kilocycle frame. The contact segment is drawn as a short arc to 
indicate a brief closure, actually lasting about 5 microseconds and occurring 
while the switch of the collector is in contact with a single segment. When 
the circuit is thus completed from a particular channel to the holding ca- 
pacitor, the voltage on the latter very rapidly assumes and then follows, 
for the remainder of the 5-microsecond interval, the potential of the com- 
pressed version of the channel signal. When the circuit opens, the latest 
state of charge, which is essentially an “instantaneous” sample, is left on the 
capacitor, and is thus held for about 16 microseconds—until the next closure. 
These sampling operations occur alternately in the two coders. 
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By a process to be considered later, each coder produces a 7-digit PCM 
code representation of its set of samples. The two coders deliver their code 
groups alternately on a common output lead during the final 10 7’, micro- 
seconds of each 16-microsecond holding interval mentioned above. In- 
dividual pulses last about 1.5 microsecond, although as delivered from the 
coders they are somewhat irregular in timing and waveform. It will be 
noted that the interval allotted to the code group from each channel is just 
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Fig. 8. Waveforms of the pulse regenerator. 

fs of the 125-microsecond frame period and that a continuous train of 
code groups is thereby produced. 

The common output circuit of the coders goes to a pulse regenerator which 
standardizes the pulses in height by slicing, and in time by gating, as illus- 
trated in Fig. 8. The peaks of the coder output pulses (line A) are lined up 
in time with the gate control pulses (line B) supplied from the timing gear. 
The latter have a constant repetition frequency of 672 kilocycles and a 
uniform pulse length of 0.4 microsecond. Accordingly, the sliced and gated 
PCM pulses (line C) are also 0.4 microsecond long, and require lengthening 
to fill their allotted 1.5-microsecond intervals. This is accomplished by a 
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circuit in the pulse regenerator which first doubles the length of each puls 
by adding thereto its own delayed reflection obtained from a short-circuite 
delay cable, and then slightly less than doubles it again by a similar proces; 
using a longer cable. A final slicing, to eliminate amplitude irregularities 
acquired in the lengthening process, yields square pulses as shown in line D 
with adjacent pulses merged into a single longer pulse. This is the fina! 
output signal delivered to the intermediate-frequency modulator. Passag: 
of the modulator output through a shaping filter results in rounded pulses 
(line £) suitable for transmission over the radio relay path. 

In this regenerating apparatus, provision is also made for introducing 
‘framing control pulse,” supplied from the timing bay and normally applied 
only to Group 1, although any other group may be used if desired. This 
pulse is about 1.5 microsecond long, and occurs once in each 8-kilocycle 
frame, but has opposite polarity in successive frames. It is timed to syn 
chronize with the first digit of the Channel 1 code and is large enough in 
amplitude to override the pulse or space put out by the coder in that position. 
Hence in the final PCM output from Group 1, pulse 1 of Channel 1 is al- 
ternately present and absent regardless of the audio signal. This arrange- 
ment, used in automatic “framing” of the receiving timing gear as described 
in the following section, thus borrows the least significant digit from one 
channel of the system, leaving that channel usable, but with 6-digit instead 
of 7-digit quality. A 4-kilocycle tone of very low amplitude which it intro 
duces in that single channel is made inaudible by the low-pass filter in the 
final audio output. 

Synchronization. ‘The connection of a transmitting channel to its proper 
receiving circuit in the time-division part of the system requires the two 
terminals to be synchronized: timing operations at the receiver must follow 
closely those at the transmitter. In a broad and general way this timing 
matter amounts to getting a local clock to keep the same time as a distant 
standard clock. Here the criterion of good timekeeping might be thought 
fussy by some standards; we cannot work with a discrepancy as long asa 
microsecond for the very good reason that incorrect routing of pulses would 
then result, associated with intolerably large decoding errors. Three pro- 
visions are made to take care of this situation. First, the framing is auto- 
matically monitored at all times. Second, if the system is out of frame—as 
the monitor 





it may be after transmission has been temporarily interrupted 
circuit hunts for and establishes synchronism. Third, whenever the system 
is not properly synchronized and framed, all message circuits are cut off to 
avoid resulting noise and crosstalk. 

For the purpose of this description we can use a mechanical analogy once 
more and picture all the transmitting channels of a time-division group 
arranged in order around a circle (Fig. 9). This time, however, we let each 
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small division of the circle constitute a commutator bar which is activated by 
the information of a particular code digit. Thus as the brush is stepped 
around at a uniform rate it puts digit information on the line in proper se- 
quence. In each complete revolution there appears a set of on-off pulses 
constituting one frame. With twelve channels in the group, and seven 
digits to each channel, one revolution of the brush arm covers the frame of 
84 digit positions. The revolution period is 125 microseconds. Roughly 
one and a half microsecond is allotted to each digit, and the brush steps 
672,(CO times a second. The driving force for stepping is supplied by a 
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Fig. 9. Representation of the method used for framing, in relation to a sirgle 12-channel 
time-division group. 


stable oscillator, its frequency determined by a quartz crystal. One of the 
digit positions is taken away from its message duties and assigned for framing 
purposes, as described in the last section. 

Our clock now is a one-handed affair with the brush arm for the hand. 
Its scale has the customary twelve main divisions, but each main division 
is divided into seven parts, one for each digit. Hour and minute markings 
refer to channel and digit, respectively. 

At the receiving end we imagine a similar clock structure to be provided. 
Input pulses go directly to the clock hand which is stepped around to dis- 
tribute pulses to the twelve message channels, plus a framing channel which 
takes up but a single digit space. Here the clock has no independent driving 
source, since it must follow the transmitter. Accordingly, the basic 672- 
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kilocycle frequency is recovered from the input pulse train by means which 
include filtering in a narrow-band crystal filter. The filter band is narroy 
enough to select the base frequency with negligible modulation and noise 
That is, sidebands produced by pulse keying are attenuated to negligible 
proportions, as is the background noise. It is this output which serves to 
step the brush arm around the frame at precisely the transmitter frequency. 

With the transmitter frequency recovered, both hands are stepped around 
their dials at the same rate. While this is a necessary condition for running 
the receiver, it is not a sufficient one since most likely the system will not 
be framed. In fact the odds are 83 to 1 against the two clocks indicating 
the same time if connection to the receiver is initiated at random times. Ii 
we had to deal with ordinary clocks, both in view, the resetting procedure 
could be accomplished by moving one clock-hand to agree with the other 
at one fell swoop. But in the PCM case resetting has to be done more dis 
cretely since only one dial position per frame is viewed in the framing re- 
ceiver. Accordingly, an orderly procedure is set up for locating the framing 
pulse which consists in examining digit positions one by one until the framing 
pulse is reached. After any one position is viewed long enough to establish 
the absence of the framing pulse, the receiving clock is set back one digit 
position and the next position viewed. 

This resetting or framing procedure is governed by the framing receiver 
through its control of a switch which connects the recovered base frequency 
to the driving mechanism of the clock. If the channels are correctly routed, 
so that it is the framing pulse which is being viewed by the framing receiver, 
the switch is left closed, and the 672-kilocycle wave steps the clockhand 
around the dial without interruption. But if the system is not correctly 
framed the framing receiver does not get its distinctive pulses. In this 
case the switch is opened every little while for the duration of a single pulse 
interval, stopping the local receiver clock during that interval while the 
transmitter clock advances one digit position. In effect the receiving clock- 
hand is set back precisely one digit interval with respect to the standard. 
Thus the next digit pulse is brought into the framing receiver. If again the 
monitored input turns out to be other than the framing pulse, the stopping 
process is brought into play once more; this process is repeated until the 
system is framed. 

As pointed out in the last section, the framing pulse is alternately absent 
and present in successive frames, corresponding to a 4-kilocycle rate. This 
is readily distinguishable from any of the message pulses, which in practice 
are found to have little energy content at this frequency. The framing 
receiver accordingly includes a resonant circuit tuned to four kilocycles 
In the hunting process, eight frame periods are allowed between successive 
interruptions of the clock drive, to give the resonant circuit sufficient time 
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to build up above a threshold when the system is first framed. With this 
time interval thus fixed, each clock position is maintained for about one 
millisecond. Hence the time required to frame the receiver varies between 
one and eighty-three milliseconds, depending upon the epoch at which the 
system connection is established. 

PCM Receiver. The received PCM signals of a 12-channel group are 
filtered from the frequency multiplex while in the intermediate-frequency 
state and then detected to the original signal band, as indicated in Fig. 6. 
Ir 


They consist of rounded pulses, nominally sinusoidal in shape, but more 
less distorted by transmission defects and accompanied by noise and inter- 
ference. 

These signals are supplied as input to the PCM receiver shown in Fig. 10. 
They are first sliced in amplitude, the slice being taken at approximately 
half the average or noise-free pulse height. Code groups of seven pulses 
are then routed alternately to two decoders, which handle even and odd 
channels respectively. ‘The routing function is represented in the drawing 
by a two-segment commutator (A) rotating at 48,000 revolutions per second. 
Before entering the actual decoding circuit, the pulses are again sliced to 
secure very great uniformity, and are gated with 0.4-microsecond, 672 
kilocycle pulses from the receiver timing equipment. Immediately after 
the arrival of the seventh digit of each code group of these standardized 
pulses, the decoding circuit produces a voltage on its low-impedance output 
lead proportional to the quantized amplitude represented by the code. As 
in the case of coding, details of the decoding process will be reserved for a 
later section of this paper. 

The decoded amplitudes are available only momentarily; therefore it is 
desirable to sample each one at the proper time and store the result as a 
charge on a holding capacitor. This sampling process is represented by 
switch B in one decoder and B’ in the other. Here the switch closures last 
only two microseconds, and values are held for about 19 microseconds. 

The next step is to assemble the six samples from odd channels held 
successively by one capacitor and the six from even channels held by the 
other into a single time-division multiplex. Switch C performs this opera- 
tion, rotating at 48,000 revolutions per second, and making contact alter- 
nately with the output circuits of the odd and even decoders. 

This 12-channel multiplex signal is passed through an instantaneous 
expandor, the purpose of which has been noted in an earlier section. To 
simplify the problem of making the input-output characteristic of this cir- 
cuit accurately complementary to that of the compressor, the two devices are 


designed to use identical silicon units. In the expandor, however, the non- 
linear device is employed in the feedback path of a broadband amplifier 
rather than in the direct transmission path, thus giving the inverse character- 
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MULTICHANNEL MODULATION SYSTEM 
istic. To allow any compressor to be used with any expandor, all the silicon 
elements are matched to a chosen standard unit, using selection and 
resistance “‘padding’’ methods. By such means, and by use of sufficient 
loop gain in the feedback amplifier, very satisfactory overall linearity of the 
system has been attained. 

At the output of the expandor the waveform of the multiplex signal is 
essentially the same as that at the input of the compressor in the trans- 
mitting terminal. The samples are distributed to their respective channel 
destinations by a distributor D, resembling the collector described earlier. 
The rotation rate is 8CCO revolutions per second, but the duration of contact 
on any one segment is only five microseconds instead of the possible full 
twelfth of the 125-microsecond frame period. This effective narrowing of 
the contact segments is done to allow the closure to occur well within the 
interval in which the circuit is completed by switch C from the output of a 
particular decoder. Each of the 12 segments of the distributor is provided 
with a holding capacitor, which stores its allotted samples for the full 125- 
microsecond frame period. The potential on any one of these capacitors 
thus changes at 125-microsecond intervals from one quantized sample ampli- 
tude to the next derived from the same original speech wave. 

This potential is of sufficient magnitude to require use of only a simple 
single-stage triode amplifier for the output of each channel. Lengthening 
the samples by holding, as described, helps to make this possible by causing 
the amplifier to deliver useful power continuously instead of on a fractional 
time basis. 

The only disadvantage of using lengthened pulses arises from an effect, 
very similar to the ‘aperture effect” encountered in sound movies, which 
introduces a curving slope across the audio gain-frequency characteristic of 
In the present case the gain drops about three decibels as the 
frequency goes from the lowest to the highest value of interest. This slope 
can be corrected by a simple equalizing network, as shown in Fig. 11. In 
the present system the equalization is incorporated in the low-pass filter at 
This is preferable to equalizing at the 


the system. 


the input of each audio channel. 
output, where power is at a premium. 

The outputs from the channel amplifiers are passed through 3400-cycle 
low-pass filters, identical with the input filters except for omission of the 
equalization, and are delivered to standard voice-frequency circuits at the 
same levels® as are provided by a type J, A, or L carrier system. 


[V. COMPONENT CIRCUITS 


Many of the circuit techniques used in the experimental system are con- 

: 7 - *y: . . 
ventional, others are more or less unfamiliar, and still others are believed to 
In the following some of the more important building blocks are 


be novel. 
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described. These include sampling circuits, the instantaneous compressor 
slicers, and the PCM coder and decoder. 

Sampling Circuits. The function of opening and closing a circuit at pr 
scribed instants, represented by rotating switches in the block schematics 
is actually performed’ by one or the other of the devices shown in Fig. 12 
employing diodes and triodes, respectively. 


The diode type (Circuit a) is normally biased ‘‘open”’ by rectified charges 


stored on the two large capacitors. While in this condition it presents « 
extremely high series resistance between the low-impedance input and th: 
load. But when a flat-topped pulse is impressed upon the pulse transforme 
the aforesaid high impedance changes to a low value (of the order of 1 
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Fig. 11. Equalization for the aperture effect. 


ohms), for the duration of the pulse. The upward-sloping tops of the pulses 
sketched in Fig. 12 represent predistortion to allow for transmission throug]! 
the pulse transformer, for the purpose of obtaining rectangular pulses at the 
tube itself. 

In the other variety (Circuit b) the plate-cathode paths of the two triodes 
are connected directly in parallel, conducting in opposite directions. The 
grids are both arranged to be biased below cut-off during the ‘“‘open” con- 
dition of the sampler by grid rectification of pulses, and to be driven strongly 
positive by a pulse when a low-resistance conducting path is required be 
tween source and load. 

These two types have their respective advantages. Tor example, the low 
capacitance to ground which the diode type affords across its load makes it 


7A single exception is the case of switch A of Fig. 10. In this case two gated slicer 
circuits are used, as described subsequently. 
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preferable for use in the collector, where twelve samplers are multipled to a 
common load. On the other hand the triode type affords a d-c. path (without 
the series blocking capacitors that are required by circuit a) between source 
and load. In cases where a holding capacitor is to be charged to a succession 
of sample amplitudes which must be kept mutually independent to avoid 
crosstalk, the d-c. path is very desirable for it avoids “memory” effects 
associated with passage of the charging currents through the blocking 
capacitors. A further useful property of the triode circuit is its ability to 
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Fig. 12. Diode and triode sampling circuits. 


sample signals of greater amplitude than that of the control pulses. With 
the diode circuit, the signal amplitude must be kept less than half the pulse 
height. 

Instantaneous Compressor. The simple configuration of the non-linear 
circuit used in the compressor and in the feedback path of the expandor is 
shown in Fig. 13. Two selected silicon rectifiers are connected in parallel, 
but poled oppositely, with a small padding resistor (R; or Re) in series with 
each. A parallel padding resistor (R3) of large value is also provided. The 
direct-current resistance of this combination varies from about 6000 ohms 
at zero signal to about 190 ohms at the peak of a full-load signal. Input is 
applied as a current through the relatively high resistor Ry, and the voltage 








28 BELL SYSTEM TECHNICAL JOURNAL 

drop across the varistor unit constitutes the compressed output. Althou:| in Fig. 1 
the development of this compressor has been a problem of many interesting amplituc 
aspects, it must suffice here to point out that copper-oxide elements ar Slicers 


unsuitable at the speeds involved because of excessive capacitance, that ingle-tr 
one : ° a Z StS 
silicon is superior to at least the presently available types of germanium ently be 


made la 
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silicon elements appears to be obtainable by the selection and resistance tive d 
padding methods employed. Temperature control is used for the sake of tude : 
constancy of the compression characteristic. Aging has been found resiet 
negligible over a period of many months. sie tubes 
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in Fig. 14 in terms of compressed signal amplitude vs. uncompressed signa! 
amplitude, both in decibels. 

Slicers. The slicer circuit shown in Fig. 15 resembles a conventional! 
single-trip multivibrator in configuration, but functions somewhat differ 
ently because of the choice of parameters. In particular, capacitor C is 
made large enough so that the potential drop across it does not vary appreci- 
ably during normal operation, and plate resistor R; is given a small value 
such that the gain around the feedback path of the circuit is approximately 
unity when both triodes are in their active regions. Germanium varistors 
VR; and VR» maintain desired bias conditions regardless of the number of 
pulses or spaces in the input signal. Unlike the single-trip multivibrator, 


which when tripped remains so until the charge on C has had time to relax, 


+150V 





SLICING = 
POTENTIAL F =F 








VR, 


- 150V 


Fic. 15. Slicer circuit. 


this circuit trips whenever the input signal falls through a narrow potential 
range near the value £, and trips back again when the input rises through 
the same potential range. A square wave of constant amplitude repre- 
senting an accurate thin slice of the input is thus made available across 
resistor Ro. The circuit is capable of high speed, slicing pulses as narrow as 
().1 microsecond. 

Time gating may readily be included among the functions of this circuit, 
through the addition of a triode having its plate and cathode connected 
directly to the plate and cathode, respectively, of tube V;. The grid of this 
added triode is held normally at ground potential, and is pulsed in the nega- 
tive direction by approximately rectangular gating pulses, having an ampli- 
tude of about 2E. The total current passing through the common cathode 
resistor Rx is essentially constant, and if either or both of the paralleled 
tubes are conducting at a given instant this current is carried by one or 


shared by both of them. Hence the tripping action, involving transfer of 
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the cathode current to V2, cannot occur as long as the grid of either V; or t 
added tube is positive with respect to the critical slicing potential. Trippi 
does occur whenever this limitation is removed, and thus the desired gati 
and slicing functions are performed concurrently. 

In the experimental system this circuit is used to regenerate the PCM 
pulses at the common output of the two coders, and again at the input 
each decoder to slice the received pulses and to sort out code groups of 0: 
and even channels. 

Coders. The method of binary coding used in this system was original); 
suggested by F. B. Llewellyn. It employs a novel electron beam tub 
This device, pictured in Fig. 16, carries out simultaneously the two functions 
of quantizing and of coding. The tube is about 103 inches long ai 
2} inches in diameter. It has the 128 combinations of the 7-digit cod 
laid out permanently as holes in an “‘aperture plate,” and translates sampk 
amplitudes from the form of beam deflections directly into PCM pulse sym 
bols. Figure 17 shows one of these tubes in its socket on the rear of a coder 
panel, and above it a rectangular permalloy shield which covers the coding 
tube of the other coder serving the same 12-channel group. 

As shown schematically in Fig. 18 the coder includes, in addition to th 
coding tube, a sampling and holding circuit which sorts out the odd (or even 
channels from the input multiplex signal, push-pull amplifiers for vertical 
and horizontal beam deflections, and simple arrangements for blanking, 
focusing and centering. Within the tube are shown, from left to right, a 
conventional electron gun, vertical and horizontal deflection plates, a rec 
tangular “‘collector” for secondary electrons, a “quantizing grid,’ the 
“aperture plate” and finally a “‘pulse plate.” Figure 19 shows the target 
end of the tube, as viewed from a point near the gun. ‘‘Digit holes” in the 
aperture plate, laid out in accordance with the desired binary code, may be 
seen behind stretched parallel wires of the quantizing grid. One may count 
64 narrow holes separated by equally narrow bars of metal in the left-hand 
vertical column, 32 holes in the next column, 16 in the next, and so on for 
seven columns. There are 129 grid wires uniformly spaced and accurately 
aligned so as to mask the upper and lower edges of every one of these holes 
when viewed from the geometric ‘point of origin” of the beam. 

Stored audio samples from the sampling and holding circuit provide poten- 
tial for the vertical deflection, with zero at the center, positive amplitudes in 
the upper half, and negative amplitudes in the lower half of the target area. 
A sawtooth sweep provides the horizontal deflection. The beam is blanked 
while deflection potentials are being changed to move it upward or downward 
from one sample amplitude to the next. When first restored, the beam 
strikes in the left-hand unperforated region of the aperture plate, and is then 
swept linearly across from left to right. Electrons which pass through the 
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into pulse codes 


This new electron device transforms speech samples 


tube 


16. Coding 


Fig. 
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digit holes during the sweep are caught by the pulse plate, forming pulses | 2" 
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beam is blanked, and is simultaneous with the application of the succeeding 


audio sample. 

The wires of the quantizing grid are used to guide the beam so that it can 
illuminate only the particular row of apertures which correspond to the 
initial vertical deflection. Without this feature, erroneous codes would be 
produced when the beam straddled the edges of apertures or crossed from 
one amplitude level into another, as a result of electrode misalignment or 
possible slight drift in potential of an applied sample. The guiding action 
(basically proposed by W. A. Marrison and applied to the PCM coder at the 
suggestion of G. Hecht) is obtained by means of feedback from the quan 
tizing grid to the vertical deflection amplifier. 
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Fig. 18. Functional schematic of the coder. 


The feedback signal is actually a current taken from the positively biased 
collector, which draws to it secondary electrons from the grid wires. The 
portion of the beam current striking the grid varies as a cyclic function of 
the vertical deflection. It follows that for some spot positions the value 
43 in the feedback loop is positive, for others negative; hence with proper 
amplifier design there is a stable and an unstable region associated with each 
wire. 

The spot can come to rest (vertically) only within one of the stable regions 
In order to locate it consistently near the center of such a region, and thus 
gain equal margins against “hopping” upward or downward aeross a wire, 
a “quantizing bias” is introduced into the vertical deflection amplifier, along 
with the feedback and the signal samples. This bias is a current of opposite 
polarity from the unidirectional feedback current, and of magnitude equal 
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to the average between the two values of feedback current which exist wh« V 








the beam falls (1) directly on a grid wire and (2) midway between two wire upV 
One may regard this bias as pressing the beam upward against a wire, resiste posi 
by downward pressure associated with the feedback current. The latte mic 
get 

wire 

blar 

Fig. 19. Interior of the coder tube, viewed from the gun end. —_ 

Int 

increases as the beam approaches the wire, and equilibrium is reached when will 

the beam is about half way between positions (1) and (2) mentioned above. gric 
The feedback current is actually intermittent, turned off and on by the 
blanking pulse, and careful analysis shows it necessary to make the bias inter poe 

mittent also, with its wave fronts synchronous with those of the feedback lor 

current. This is readily accomplished by deriving the bias from the blanking _ 

signal itself. Cire 
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With this arrangement the beam, suddenly turned on, moves either 
upward or downward from its initial “unquantized” position to the nearest 
position of stable equilibrium. Quantization is completed in less than a 
microsecond. Thereafter, as the beam is swept horizontally across the tar 
get area, it remains pressed upward against the lower surface of its guiding 
wire. Quantization is thus maintained until the end of the sweep, when 


blanking occurs. The margins against hopping across wires, while quanti 





Fig. 20. ‘Television picture” of the aperture plate of a coder tube 





zation is in effect, are related to the amount of loop gain in the feedback path. 
In the present system feedback measuring about 20 decibels is provided, and 
will counteract changes which, without feedback, would move the beam two 
grid wires in either direction from the initial position. 


, 
; This coding process, based on the electron beam coding tube and making 
" use of a two-dimensional permanent layout of the code, is more straight- 
7 forward than the various coding processes which depend upon counting or 


sequential comparison. Accordingly it leads to higher speeds and greater 
circuit simplicity. 





Figure 20 illustrates the coding accuracy obtained. This photograph of 
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the screen of a test oscilloscope may be thought of as a television picture ¢ 
the aperture plate of the coding-tube. To produce the pattern an audi 
frequency sawtooth wave of full-load amplitude was applied to the samplin; 
and holding circuit at the input of a coder. The resulting sample ampli 
tudes, quantized by the coding tube and falling into all the possible 128 step 
of the quantizing characteristic with uniform regularity, were used to ener 
gize the vertical deflection of the test oscilloscope. An ordinary synchro 
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Fig. 21. Shannon decoding circuit and waveforms. 


nized sweep provided the horizontal deflection. The square PCM pulses, 
delivered by the coder via the pulse regenerator, were applied to the intensity 
control. 

Thus the code pulses corresponding to each quantized amplitude were 
made to appear as a row of blanks in a horizontal trace at the proper relative 
height. This pattern is very useful in studying coder performance. 

Decoders. The decoding method is an impressively simple one originally 
proposed by C. E. Shannon. In its basic form it employs a pulsed resistance- 
capacitor circuit as illustrated in Fig. 21. Upon arrival of each pulse of the 
code, an identical increment of charge is placed upon the capacitor of the 











easininaeinbieias Es 





lic 
ins 
ph 


ep 


iro 


ses, 


‘ity 


ere 


ive 


lly 


ce- 
the 
the 











MULTICHANNEL MODULATION SYSTEM 37 
device. The time constant ¢ = RC is such that, during any single pulse 
interval, whatever charge is on this capacitor decays precisely 50°; in ampli- 
tude. It follows that the charge remaining at some chosen instant after 
the arrival of a complete code group consists of contributions of all its pulses, 
weighted in a binary manner. That is, if we detine the contribution of a 
pulse in the final digit position as 3, then contributions of }, §, ye, 3'z, 
gy and yes, respectively, are made by pulses in successively earlier posi- 
tions. Any value from 0 to }3§, in steps of +}s, may thus be produced 


Of course the digit holes in the aperture plate of the coder are laid out to 
make this straight-forward scheme workable. Since samples of low-level 
audio signals are coded near the center of the aperture plate, the corre- 
sponding decoded values lie in the neighborhood of }. 

The basic Shannon decoder, then, comprises the resistance-capacitor cir- 
cuit, means for supplying it with precisely controlled units of charge at 
precisely determined times, and a sampling and holding circuit (repre- 
sented by switch B or B’ in Fig. 10) to measure and store the decoded 
potential which is fleetingly present across the capacitor at a regularly 
recurring instant 7, following the final pulse position. The scheme 
employed to inject the identical charges involves a regulated source of 
current and a gate to admit this current to the resistance-capacitor circuit 
under control of the regenerated PCM pulses. Two successive slicing opera- 
tions and careful gating, as described earlier, make these pulses more than 
adequately uniform. 

The wave form sketches of Fig. 21 show three typical decoding cycles. In 
the first, a single pulse in digit position 7 produces a decoded amplitude of 
+; in the second, a pulse in position 1 gives +35; and in the third, pulses at 
2, 5 and 7 provide a decoded value of §}. It may readily be verified that 
the provision of an idle channel interval between operations (following from 
the alternate use of two decoders) allows the residue of one decoding opera- 
tion to decay to a negligible value (never larger than y$g of a single step 
height or “‘quantum’’) by the time of the next consecutive sampling. 
Experimentally, interchannel crosstalk from this source is virtually non- 
existent. 

In the foregoing it has been emphasized that precise timing is required for 
this basic Shannon decoder. Although the necessary accuracy was actually 
obtained without great difficulty in early tests, a modification has also been 
introduced which eases the requirements to a very marked extent. This 
scheme, devised by A. J. Rack, employs a damped resonant circuit in con- 
junction with the resistance-capacitance elements in a manner illustrated 
by Fig. 22. The natural frequency of the resonant circuit L, C2, R2 is made 
equal to the PCM pulse rate, and the time constant of the damped oscillation 
(t=2R.C2) is matched to that of the circuit R;, C;. The same charging 
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pulses pass through both sections of the circuit; hence by proper choice ¢ 
C, and C, the amplitudes of the damped sine wave and the exponential ma 
be proportioned so that the rate of change of their combined potenti 
becomes zero at successive points one pulse period apart. In fact it has bee 
found possible to make both the first and the second time derivatives « 
potential equal to zero simultaneously at such points. 
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Fig. 22. Shannon-Rack decoder. 





Fig. 23. Output of Shannon-Rack decoder for a signal giving 100° modulation. 


By this modification, not only the times of application of the charges but 
the time of sampling is made much less critical. Of course this presumes 
that the sampling circuit is designed to complete its operation near the center 
of a level region. In Fig. 22, the voltage transient due to a typical pair of 
pulses is sketched, and Fig. 23 shows an oscilloscope screen, on which the 
waveforms delivered by the cathode follower of one of the Shannon-Rack 
decoders of the system are superimposed for a full set of 128 possible code 
combinations in sequence. 
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V. PERFORMANCI 


In general the behavior of the system has shown promise for toll plant 
application. Among other things the stability realized in the adjustments 
of the coder and decoder, and the apparent absence of aging or other drift 
in the compressor and expandor have been gratifying. A daily check of the 
focus and centering in the coders and of the time constants in the decoders 


appears adequate to keep them in optimum adjustment. The synchronizing 





Fig. 24. Output of decoder (vertical) vs. input to coder (horizontal). 


gear has been found equally easy to maintain. On the other side of the 
ledger, occasional breakdowns have pointed up the need for alarm and auto- 
matic replacement facilities in any version of the system which might be 
developed for commercial service. 

A few measured characteristics are given in this section in addition to the 
compression and audio gain-frequency characteristics already shown (Tigs. 
14 and 11, respectively). 

Input-Output Characteristics. The diagonal trace in the oscilloscope 
pattern of Fig. 24 shows the relationship between the input of a coder (hori- 
zontal deflection) and the output of the corresponding decoder (vertical), 


For this pattern a full-load audio signal was applied to the input of the odd 
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coder only (without passing through the compressor), and the output was 
taken directly from the common output of the two decoders. Thus the sis 
odd channels took turns transmitting the signal while the six even channe!s 
produced the horizontal center-line. Uniform quantizing steps may be se¢ 
along most of the trace, but are obscured near the ends by defocusing of t! 
test oscilloscope. 

A similar pattern, obtained with the compressor and expandor included 
in the transmission path, is shown in Fig. 25. Here tapered steps may | 
discerned, as well as a small amount of non-linearity due to residual impe 





Fig. 25. Output of expandor (vertical) vs. input to compressor (horizontal) 


fections in the companding. The effect of the channel peak choppers is not 
included in this pattern. 
Two measured overall input-output characteristics appear in Fig. 26, fo: 


the case of a typical single channel and for five channels patched in tandem 


through 17-decibel pads on a 4-wire basis. The latter simulates a possibl 


extreme case of a long circuit in which for some reason it is desired to decod« 


to audio at four junction points as well as at the final terminal. It should 


not be confused with a series of spans between which the PCM pulses ar 
amplified or regenerated without decoding. In the latter case, of course, th 
overall audio characteristics are independent of the number of spans. 


Audio Noise. (Quantizing was found to be the only significant source o! 
noise in the received audio signals. Noise levels measured in the absence oi! 


speech are shown in Fig. 27. The measurements are given for various num 
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» bers of channels from one to ten, connected in tandem as described in the 
preceding paragraph. Two scales of ordinates are shown in this figure. On 
( ¢ — 
| a" T 
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“f* s 
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Fig. 26. Input-output characteristics of PCM channels 
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Fig. 27. Noise measurements on idle PCM channels 


the left is a reference scale of weighted noise employed in the Western 


Electric 2B Noise Set; on the right, a scale which relates the corresponding 
unweighted root-mean-square noise to the nominal maximum undistorted 
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single-frequency output of the system (taken to be 9 decibels above a milli- 
watt at zero level). This output corresponds to an input just reaching the 
peak limiters. Thus in a single link the idle circuit noise is down about 68 
decibels from the full-load sine wave. For five links in tandem, the readings 
are at least 8 decibels below the accepted limit (29 decibels on the left-hand 
scale) for noise at the end of a 4000-mile circuit. Quantizing noise is found 
to increase approximately three decibels for each doubling of the number of 
links, as is generally the case with other forms of random noise. 

For a single channel, or a small number of channels in tandem, the id 
circuit noise varies considerably with the vertical centering of the coding 
tube. This may be understood by noting that if the zero-signal operating 
point effectively rests near the center of a ‘‘tread’’ on the quantizing stair- 
case, (fig. 4) a small amount of power hum or other disturbance may simp), 
move it back and forth on the same tread, in which case quantizing noise is 
entirely absent. On the other hand, if the operating point is near a “riser” 
the small disturbances may cause it to joggle from one tread to another, 


producing noise. The measurements given in Fig. 27 were obtained with a 
very quiet input circuit and with the centering adjusted for maximum noise. 
The “‘joggling”’ was produced largely by residual power hum. 

A-B tests to compare PCM transmission over a single link with direct 
transmission over a noise-free circuit of the same audio band were carried out 
using a wide range of talker volumes. In such tests only a few experienced 
observers were able to pick the PCM path consistently. When the PCM 
circuit included five tandem links most observers could tell the difference, 
but all judged the quality to be more than satisfactory for toll service. 

Crosstalk. It has been pointed out earlier in this paper that interchannel 
crosstalk in a PCM system can occur only in the terminal equipment. 
Considerable care was exercised, particularly in the design of the time- 
division parts of the system, to hold individual sources of crosstalk to 70) 
decibels or better. As a result, measurements using a single-frequency test 
tone and a current analyzer have shown the overall crosstalk from any one 
channel to any other to be down 66 decibels in the worst cases. 

Very severe tests have also been made in which a loud talker was connected 
to ten channels of a group simultaneously and crosstalk into either of the 
two remaining channels (one odd, one even) was listened to, and measured 
with the 2-B Noise Set. In such tests unintelligible crosstalk could be de- 
tected, which seemed to consist of changes in the quality of the quantizing 
noise occurring at a syllabic rate. The 2-B readings averaged about a 
decibel above the quantizing noise of a single idle channel with occasional 
peaks reaching the 17-decibel point on the reference scale. 

In tests involving 24 talkers in 12 simultaneous conversations, crosstalk 
was practically undetectable. 


ie PLE at Eee 


See ea Se a ae 





% 
: 
4 
a 
<I 
a 
i 
( 











the 
oul 


su] 
to 

the 
pe 
det 
wh 
dec 
onl 
wa 
of 

tra 
pri 
rou 
ide 


Wi 


pu 
an 
fer 
Kl. 
of 

me 


sta 


C( 
he 
ec 
le- 


ng 


lk 


eS. 











UULTICHANNEL MODULATION SYSTEM $3 


Radio Interference and Noise. Yo obtain experimental contirmation of 
the expected tolerance to high interference levels in the radio path, the 
output of an oscillator, tunable through the band near 65 megacycles, was 
superimposed upon the received intermediate-frequency signal at the input 
to the group selection filters. With this controllable interference tuned near 
the center of the Group 1 filter, and its amplitude set 6 decibels below the 
peak amplitude of the (noise-free) pulses, errors were so plentiful that the 
demodulator did not remain synchronized. Proper framing was restored 
when the amplitude difference was increased to 7 decibels, but enough 
decoding errors remained to give intolerable audio noise. At & decibels 
only an occasional crackle of noise was observed, and at 9 decibels reception 
was perfectly normal. Similar tests of Group 2 gave the same results except, 
of course, that synchronization was not affected. The fact that noise-free 
transmission was not maintained quite up to the ideal 6-decibel point is due 
principally to the width (0.4 microsecond) of the gate which is applied to the 
rounded PCM pulses entering the receiving equipment. If necessary the 
ideal could undoubtedly be approached more closely by reducing this 
width, thus admitting only the extreme peaks and troughs of the signal. 

The effects of actual fluctuation noise were studied by sending the PCM 
pulses over the radio path at reduced level. The boundary between good 
and bad transmission was not so sharp as with the continuous-wave inter- 
ference, as should be expected because of the random nature of the noise. 
Flawless reception occurred when the root-mean-square signal at the peak 
of a pulse was greater by 18 decibels than the root-mean-square noise, both 
measurements being made at the output of a group selection filter. 
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Electron Beam Deflection Tube for Pulse Code Modulation 
By R. W. SEARS 
INTRODUCTION 
ULSE code transmission systems! in which successive signal amplitude 
samples are transmitted by pulse code groups require special modulators 
The essential operational requirements of a pulse code modulator are: (1 
to quantize or measure the signal amplitude sample to the nearest step in 
the discrete amplitude scale transmitted by the pulse code system, and (2 
to generate the group of on-off pulses identifying the step. 


>» 4in which quantization and pulse 


Several methods have been proposed’: ‘ 
formation were performed with circuits employing conventional electron 
tubes. The circuits involved sequential and comparison operations and 
were not easily adapted to a multi-channel time division system because of 
limitations in coding speeds and the complexity of the equipment. An 
electron beam deflection tube has been developed which, together with as- 
sociated beam positioning and sweep circuits, performs the modulation 
rapidly, making possible the sequential modulation of a number of channels 
in time division multiplex. 

The electronic principles, design and characteristics of the experimental 
tube are described in the present paper. 


CONVERSION FROM SIGNAL INPUT TO PULSE CODE OUTPUT 


An electrical input voltage may be converted into an output code pulse 
group with the electron beam deflection tube shown in Fig. la. An aperture 
or code masking plate is arranged perpendicular to the axis of the electron 
gun at the focal point. The coordinates of the aperture plate are aligned 
with the deflection axes of the XY and Y deflector plate pairs. The electron 
beam strikes the output plate when it passes through an opening in the 
aperture plate. 

An input voltage of appropriate value applied to the V deflector plates 
will deflect the beam to point ‘‘a” of the aperture plate as indicated in Fig. 
la. A linear sweep voltage applied to the X deflection plates, while the 

1 An Experimental Multi-Channel Pulse Code Modulation System of Toll Quality, 
L. A. Meacham and E. Peterson, this issue. 

2 A. H. Reeves, U.S. Patent * 2,272,070, Feb. 3, 1942. 

3H. S. Black and J. O. Edson, paper presented June 11, 1947 at A. I. E. E. meeting; 


Montreal, Canada 
4W. M. Goodall, Bell System Technical Journal, July 1947. 
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input voltage on the Y deflection plates is held constant, causes the electron 
beam to sweep across the aperture plate along the dashed line a-6. A time 
sequence of output pulses is produced at the output plate when the beam 
passes through the apertures of the code pate along the path a-. 

A series of output pulses or a ‘“‘pulse group”’ is characterized by the pres- 
ence of pulses at time positions corresponding to the several vertical columns 
of apertures. The code plate shown in Fig. 1 is laid out in accordance with 

















| = 





Fig. 1—Electron beam deflection tube for coding. 


the binary number system pattern’ which was chosen for the present pulse 
code work because of the simplicity of decoding. The four-digit code plate 
shown in the tube of Fig. 1 provides for coding only 16 amplitude values 
and was used to facilitate the illustration. The tube developed for the ex- 


5 The four vertical columns of the four-digit code plate in Fig. 1b provide four positions 
in time for output pulses. Vertical columns marked 1, 2, 3 and 4 correspond to the first 
second, third and fourth digits, respectively, of the pulse code transmitted. Pulses lo 
cated in time in accordance with this notation are given “weights” of 1, 2, 4 and 8, respec 
tively. Sixteen pulse group combinations are indicated by the 16 dashed horizontal 
lines in Fig. 1b. The pulse group defined by beam sweep a corresponds to a total “weight” 
of 15; beam positions 8 and y correspond to total “weights” of 14 and 13, respectivel 
and so on, with the bottom horizontal dashed line in the figure corresponding to zero. 
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perimental pulse code system uses a seven-digit code plate which provices 
for coding 128 amplitude values. 

The vertical distances between the centers of successive code sweep px 
tions (horizontal dashed lines in Fig. 1b) are made equal so that the codes 
are spaced by equal input voltage steps. A continuous range of input sig- 
nal amplitudes will result in a continuous range of horizontal sweep positions 
With an infinitely small electron beam, input signal amplitudes in the range 
from o to y; will produce a single output pulse group and input signal sam- 
ple amplitudes from y; to ye will produce another output pulse group. This 
process of dividing the total input amplitude range into finite steps and ar- 
ranging that input voltages falling within each step produce one and only 
one output pulse group is called quantization. 

The tube of Fig. 1 will only quantize effectively if the electron beam is 
infinitely small and the sweep and aperture plate axes are aligned exactly 
With a finite beam size, there will be sweep positions for which the beam 
straddles and sweeps out a combination of two adjacent codes. 

Precise quantization and a uniform pulse output are required. The 
problems of quantizing, alignment and uniform pulse output have been 
solved by the use of a wire grid, called the quantizing grid, located in front 
of the aperture plate. 


QUANTIZATION OF BEAM POSITION BY FEEDBACK 


The quantizing grid consists of a horizontal array of grid wires aligned 
parallel to the code sweep or X axis of the aperture plate. The grid spacings 
and alignment are such that a wire lies between each adjacent pair of code 
groups as viewed from the various incident angles of the deflected electron 
beam. The quantizing grid, by means of an electrical feedback path to the 
signal deflection plates, divides the input signal range into a number of 
equal steps and positions the electron beam tc the proper level for the code 
corresponding to the voltage step within which the signal amplitude sample 
falls. The quantizing grid wires also constrain the electron beam during 
the formation of the output code pulses, so that it must sweep out the code 
initially selected. In general, wires or shaped electrodes of any sort located 
where the electron beam can impinge thereon and connected in feedback 
relation to the deflection system constrains the electron beam to move in 
patterns prescribed by these electrodes and are thus called beam guides. 

The coding tube with quantizing grid and feedback circuit is shown 
schematically in Fig. 2. The electrode line-up, reading from left to right 
in the figure, consists of an electron gun, deflection system, secondary elec 
tron collector, quantizing grid, aperture plate and output plate. 

For the present purpose, a consideration of the collector and output plate 
electrodes is omitted and it is assumed that the grid does not emit secondary 
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electrons when bombarded. Electrons which strike the grid produce a cur- 
rent in the grid circuit while electrons that miss the grid have no effect in 
the grid circuit. The electron beam current intercepted by the grid will be 
dependent on the y deflection or position of the electron beam. The cur- 
rent is a maximum when the beam is centered on a grid wire and a mini- 
mum when it is centered between two grid wires and varies with beam posi- 
tion as indicated by the curve in the lower portion of Fig. 3. The curve is 
constructed for the case in which the beam diameter is slightly greater than 
the space between two grid wires. The current to the grid never becomes 


zero for any beam position. It may be thought of as having a “‘d-c. com- 


ponent” B. Amplitier 2 in Fig. 2 introduces a bias which cancels this ‘‘d-c. 


component” so that the feedback voltage is symmetrical about zero. 


! INPUT aoa a. 
' VOLTAGE APERTURE OUTP 











AMPLIFIER ’ 
1 
Mm v 


\ ee oh Sn ee 7 £5 ee eee 
b — > , ( \ aa 
~~ | oma } ea —>" | ” ” 




















BEAM FEED-BACK VOLTAGE eo 


<«— 











J AMPLIFIER 
2 
2 


Fig. 2—Coding tube with quantizing grid and circuit schematic. 
The grid current is amplified and the voltage developed is applied in feed- 


back relation to the Y deflection plates as shown in Fig. 2. 


The case in 
which a positive feedback potential deflects the beam in the same direction 
as that for a positive signal voltage will be considered. 

The beam deflecting voltage is equal to the sum of signal and feedback 
voltages and the beam position is a linear function of the deflection voltage 
so we may write 


—et+ Dy=e (1) 


where ¢ is the input signal voltage, D the deflection constant, y the beam 
deflection or position and e; the feedback voltage. The feedback voltage 
ey is a periodic. function of beam position y. Equation 1 therefore defines 
equilibrium beam positions for input signal voltage e. Equilibrium beam 


positions in accordance with Equation 1 are determined graphically in the 
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top portion of Fig. 3. The feedback voltage representing the right-h: posit 
side of Equation 1 is plotted as a function of beam position. The electro: 7 conti 
beam will have several possible positions of equilibrium at points (p1, )2 3 Th 
3, p4, p5, p6 and p7) where the deflection line D erected from —e, repre. samp 
senting the left-hand portion of Equation 1, crosses the feedback curve ; slope 
These are the only beam positions for which the deflection potential (signa S +ium 
plus feedback) attains correct values for corresponding beam positions = This 
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Fig. 3—Graphical representation of quantization. ; age é 
the \ 
However, only positions 1, p3, p5 and 97 will be in a true state of equilib orde 
rium. This can be seen as follows: Consider the beam at point p3; if the | grou 
beam is perturbed toward the right, the feedback voltage tending to deflect 
the beam to the left increases. The opposite action ensues if the beam is 
perturbed in the left-hand direction. Point 3 is a true equilibrium point. 
On the other hand, consider the beam at position p4. If the beam is per ry 
turbed to the right, the feedback voltage tending to deflect it to the right ; syste 
increases and the beam continues to move until it reaches the equilibrium } Pict 
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position p5. If the beam is perturbed to the left from position p4 it will 
continue to move to the left until it reaches position p3. 

The number of possible equilibrium beam positions for an input signal 
sample depends on the maximum values of the feedback voltage and the 
slope of the deflection characteristic. It is necessary that only one equilib- 
rium beam position be available for a particular small range of signal voltage 
This can be achieved as follows: if a signal voltage e is established with 
the feedback circuit inoperative, the beam will be at a position p, Fig. 3. 
When the feedback circuit is activated with the signal voltage held at e, the 
beam will move from p to ~5. With this procedure, signals in the range 
from e; to é2 will result in equilibrium beam positions between points / and 
ron the curve. ‘Thus, for signal voltages within the range from e; to és the 
beam will fall in the small spacial interval m, whereas beam positions for 
input signals in the same range without feedback would vary from grid 
wire 5 to grid wire 6. Likewise, signal voltages between e, and e; will cause 
the beam to assume positions between s and / in the spacial interval ». 
The electron beam may be thought of as “leaning” on one side of a grid 
wire for a finite signal voltage range and on the same side of the next grid 
wire for an adjacent signal voltage range. 

If the feedback voltage is of the opposite polarity to that assumed above, 
the quantizing action proceeds in the same manner except that the quan- 
tized beam positions lie at the left of the wires. ‘The beam may be thought 
of as “leaning” on the opposite side of the grid wire. 

The proper quantizing action is obtained by establishing and holding the 
signal voltage with the feedback circuit inoperative and then activating this 
circuit. The feedback circuit may be deactivated and activated by either 
(a) blanking and deblanking the electron beam, or (b) defocusing and 
focusing the electron beam by applying the proper voltage change to the 
beam control or focusing electrodes of the gun, respectively. 

Since the grid wires are parallel to the horizontal rows of aperture holes, 
the feedback action constrains the beam to sweep out the code group initially 
selected even though the sweep axis is tilted slightly with respect to the 
grid wires and aperture plate. The maximum swing of the feedback volt- 
age at the deflection plates should be about three or four times the value of 
the voltage required to deflect the beam from one code group to the next in 
order to provide ample protection against the beam jumping from one code 
group to the next code group during the sweep. 


THE EXPERIMENTAL CODING TUBI 


The seven-digit experimental tube developed for pulse code transmission 
system trials utilizing the electrode lineup shown schematically in Fig. 2 is 
pictured in Fig. 4. The electron gun assembly and the target plate assembly 
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Experimental seven-digit coding tube. 
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are sealed in at opposite ends of the tube envelope. The over-all length of 
the tube is 11}” with a maximum bulb diameter of 2}” 

The electron gun operates at a tinal anode potential of 1000 volts with a 
beam current of approximately 10 microamperes. A potential of about 100 
volts applied to the signal input deflection plates deflects the beam from the 
center to the top of the aperture plate. This corresponds to a maximum 
deflection angle of 103 

The four electrodes of the target assembly, secondary collector, quantizing 
grid, aperture plate and output plate are shown in the photograph of Fig. 5 
from bottom to top respectively. The secondary collector is a simple rec- 
tangular shaped electrode. The quantizing grid consists of a circular frame 
with a parallel array of grid wires stretched across a rectangular opening. 
The aperture plate is a thin disc with apertures arranged in a binary pattern 
which provides for a seven-digit code. The output plate is a thin circular 
disc. Both aperture and output plates are coated with a carbon layer to 
suppress secondary electron emission from their surfaces. 

The parts of the target assembly are held in accurate alignment in a jig 
and cemented and held in position on four ceramic rods. The entire as- 
sembly is held rigidly in the tube envelope by means of spacers attached to 
the quantizing grid and output plate. 

The target assembly is aligned with the electron gun and deflection plate 
axes by means of lineup tools in the glass lathe at the time the final seal is 
made at the center of the glass envelope. It has been possible to hold the 
alignment of the deflection axes with the aperture plate to within slightly 
less than 1° with this construction. 

The construction of the quantizing grid may be seen more clearly in the 
photograph of Fig. 6. The grid frame has raised portions on two sides of 
the rectangular opening. These are milled with a series of grooves for each 
grid wire. The grid laterals are affixed in the grooves by brazing and are 
thus accurately spaced with respect to each other and to assembly lineup 
holes which can be seen spaced around the edge of the grid frame. 
The wires are held taut by means of a flat spring which is welded to the 
grid frame and supplies tension to stretch the lateral wires. The grid wires 
are 4.0 mils in diameter, processed to have a secondary emission coefficient 
of about 3. The laterals are spaced 11.6 mils between centers. 

The openings in the aperture plate are made by a punching operation. 
The area of the aperture plate covered by the first and second digit columns 
is milled to a thickness of 5 mils in order to facilitate accurate punching of 
the smallest apertures. The apertures in the first digit column (bottom 
horizontal row in Fig. 5) are rectangular .012” x .062”. The seven-digit 
code pattern provides 128 different output pulse groups. The wide openings 
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Fig. 5—Target electrodes. 
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wn 
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at the left hand side of the aperture plate in Fig. 5 provide a peak amplitude 
range for which the output pulse group consists of all seven pulses. This in 
effect provides a peak limiting action. 

Leads from the four electrodes of the target assembly are brought out 
directly to terminal caps on the side and end of the tube envelope to decrease 
the interelectrode capacitances and to facilitate direct connection to external 
circuits. 








Fig. 6—Quantizing grid. 


OPERATIONAL AND DESIGN CONSIDERATIONS 


The quantizing action depends on the periodic variation of electron beam 
current intercepted by the wires of the quantizing grid with beam position. 
Secondary electron emission from the grid and other electrocles was assumed 
to be negligible with the feedback circuit connected directly to the quantiz- 
ing grid as shown in Fig. 2. The uniform suppression of secondary electrons 
is difficult to achieve even though the grid is coated and processed for a low 
secondary ratio. It is also difficult to prevent secondary electrons produced 
at the aperture plate from being collected by the grid. 

A preferred method of operation utilizes the secondary electrons produced 
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at the grid by the impinging primary beam for the quantizing action rai her 
than the directly intercepted electron current as heretofore assumed. r 
secondary electron collector located in front of the quantizing grid is n 
tained at a positive potential and collects most of the secondaries from ‘hx 
grid. There is, of course, a correspondence between the secondary electro: 
current and the fraction of the beam current intercepted by the grid wires 
The quantizing circuit is made by connecting the feedback path to the 
secondary collector and the quantizing action proceeds as described hereto 
fore. This method has the following advantages over the direct prim 
current method: (1) the collector current as a function of beam position js 
much more regular; (2) the swing between maximum and minimum current 
is considerably larger because of the secondary emission multiplication at 
the grid surface; and (3) the capacitance between collector and ground is 
lower than the capacitance between the closely spaced grid and aperture 
plate. 

With secondary electron current feedback, the aperture plate is operated 
at a positive potential relative to the grid to suppress secondary electrons 
from the aperture plate. The proportion of the secondary emission from 
the grid collected by the aperture plate is small compared with that collected 
by the secondary collector. High velocity secondaries originating at the 
aperture plate are, however, able to penetrate the retarding field and strike 
the grid. These energetic secondaries produce low-velocity secondaries at 
the grid which flow to the secondary collector electrode. This alters the 
character of the secondary collector current in accordance with the spacial! 
pattern of the apertures in the code plate. The surface of the aperture plate 
is carbonized to reduce the emission of high-velocity secondaries. The spa- 
cial variation of the quantizing current is reduced to less than 10° of 
the total quantizing current swing with this treatment. 

The periodic nature of the quantizing current with beam position must be 
as uniform as possible as regards both the ‘“‘a-c” and ‘d-c’? components. 
The maximum swing should also be as large as possible to permit the use 
of the lowest possible impedance in the feedback path to obtain a wideband 
characteristic for fast coding. 

The factors which determine the maximum current swing are electroa 
beam current, secondary emission coefficient of the grid wires, electron 
beam spot size, grid wire diameter and spacing. The quantized beam falls 
approximately halfway between the center of a grid wire and the midpoint 
between wires. The beam must be small enough that its edge does not 
extend into the region beyond the wire on which it “leans” by an appre 
ciable amount. This is obtained with the electron beam focused to a spot 
slightly smaller than that for maximum quantizing current amplitude. 
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Optimum performance has been obtained with the electron beam focused 
to a radius® of about 5 mils for the grid spacings prev iously specitied 

he principal irregularities in the periodic quantizing current are caused 
by variations in secondary emission coefficient of the grid and aperture 
pte and deflection defocusing of the electron beam. The secondary 
emission from the grid surface can be made sufficiently uniform by careful 
processing. This is illustrated in Fig. 7 which is a trace of an oscilloscope 
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Fig. 8—Deflection defocussing 


presentation of the secondary collector current as the beam is deflected by a 
linear sweep at right angles to the grid laterals. The trace was limited to 
cover a sufficiently small number of wires to show details of its shape. 
Figure 8a illustrates the effect of deflection defocusing. The sweep was 
expanded to cover all of the grid wires and the electron beam focused for 
maximum amplitude in the center region. The curve is so compressed that 


Distance from the center of the beam to a point at which the electron current density 


is 5°% of its value at the center. 
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the individual oscillations of the current between maximum and mini: 
values can hardly be resolved. The envelope of the curve indicates 
extent to which the swing is reduced at the two ends by the increa 
electron beam spot size which results from deflection defocusing. In 
8b, the focusing voltage has been changed by 12 volts and it can be 
that the beam is in focus at the maximum deflection angles and ou 
focus in the center region. With an intermediate or compromise focus \ 


age, the tube will quantize satisfactorily over the entire code range. bes 


results have been obtained with the tube in the experimental pulse « 
system by the use of a simple circuit which changes the focus voltage 

linear manner with the rectified or absolute value of the input deflect 
signal thus compensating for the deflection defocusing of the tube. 





Fig. 9—Typical pulse code outputs. 


The output plate is usually operated at a positive potential relative to t! 
aperture plate with attendant suppression of secondary emission from t! 
former. Output pulses of opposite polarity may be obtained by operati 
the output plate at a negative potential and processing this plate to have 
secondary emission ratio greater than unity. 

The several groups of output pulses shown in Fig. 9 are illustrative of t! 
tube output. Successive digit pulse positions occur from left to right 


the figures. The small kink in the base line at the first and second digi 


e 


+ 


pulse positions in codes a and b respectively are small error pulses caused 


by a small portion of the edge of the beam overlapping a grid wire when th: 


beam is quantized. Error pulses of this magnitude are readily eliminat: 
by slicing when the tube is used in the pulse code transmission system. 
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Coding tubes have operated satisfactorily for long periods of time in the 
experimental multi-channel pulse code system and have required minor 
adjustments of potentials no more than once a day. The tube operates in 


this equipment with a code sweep time of 10 microseconds. 
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Metallic Delay Lenses 


By WINSTON E. KOCK 


\ metallic lens antenna is described in which the focussing action is obtains 
by a reduction of the phase velocity of radio waves passing through the lens rath 
than by increasing it as in the original metal plate lens. The lens shape accor: 
ingly corresponds to that of a glass optical lens, being thick at the center an: 
thin at the edges. The reduced velocity or ‘‘delay” is caused by the presence o 
conducting elements whose length in the direction of the electric vector of th: 
impressed field is small compared to the wavelength; these act as small dipoles 
similar to the molecular dipoles set up in non-polar dielectrics by an impress« 
field. The lens possesses the relatively broad band characteristics of a solic 
dielectric lens, and since the conducting element can be made quite light, the 
weight advantage of the metal lens is retained. Various types of lenses ar 
described and a theoretical discussion of the expected dielectric constants is given 
An antenna design which is especially suitable for microwave repeater applicatior 
is described in some detail 
INTRODUCTION 
HE metal lens antennas described by the writer elsewhere! comprised 
rows of conducting plates which acted as wave guides; a focussing effect 
was achieved by virtue of the higher phase velocity of electromagneti 
waves passing between the plates. Higher phase velocity connotes an 
effective index of refraction less than unity, and a converging lens therefore 
assumes a concave shape. The relation between the index of refraction », 
the plate spacing, a, and the wavelength \ 


n= Y1— (X/2a)?, l 


indicates that the refractive index varies with wave length. As a cons 
quence, such lenses exhibit ‘chromatic abberation”’; i.e. the band of fre- 
quencies over which they will satisfactorily operate is limited. Although 
some of these lenses may have ample bandwidth (15° to 207) for most 
microwave applications, others, having large apertures in wavelengths, may 
have objectionable bandwidth limitations. For example, the lens of Fig. 1, 


having an aperture diameter of 96 wavelengths, has a useful bandwidth « 
only 5%. 


As a means for overcoming these band limitations, the metallic lenses « 


this paper were developed. They are light in weight and possess an index « 


refraction which can be made sensibly constant over any desired band of 


microwave frequencies. They avoid the weight disadvantages of glass or 


4 


plastic lenses, and retain the tolerance and shielding advantages of the le: 


IW. E. Kock, Bell Laboratories Record, May 1940, p. 193; Proc. I. R. E., Nov. 194 
p. 828. 
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ver the reflector antenna. Because electromagnetic waves passing through 
them are slowed down or “delayed” (as in the glass lenses of optics), they are 


alled delay lenses; and since the elements in the lenses which produce the 


delay are purely metallic, they are called metallic delay lenses. 
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Fig. 1—Waveguide metallic lens having an aperture of 96 wavelengths and a useful 
ol bandwidth of 5°%. 





PART I—EXPERIMENTAL 


FUNDAMENTAL PRINCIPLES 





The artificial dielectric material which constitutes the delay lens was 
arrived at by reproducing, on a much larger scale, those processes occurring 
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in the molecules of a true dielectric which produce the observed del: 
electromagnetic waves in such dielectrics. This involved arranging me 
elements in a three-dimensional array or lattice structure to simulat 
crystalline lattices of the dielectric material. Such an array respon: 
radio waves just as a molecular lattice responds to light waves; the 
electrons in the metal elements flow back and forth under the action of 1 
alternating electric field, causing the elements to become oscillating dij 
similar to the oscillating molecular dipoles of the dielectric. In both cases 
the relation between the effective dielectric constant e of the medium, th: 
density of the elements .V (number per unit volume) and the dipole stre: 
(polarizability a of each element) is approximately given by 


€= €¢ + Na ) 


where €o is the dielectric constant of free space. 

There are two requirements which are imposed on the lattice structure 
First, the spacing of the elements must be somewhat less than one wave 
length of the shortest radio wave length to be transmitted, otherwise dii- 
fraction effects will occur as in ordinary dielectrics when the wavelength is 
shorter than the lattice spacing (X-ray diffraction by crystalline substances 
Secondly, the size of the elements must be small relative to the minimum 
wavelength so that resonance effects are avoided. The first resonance 
occurs when the element size is approximately one half wavelength, and for 
frequencies in the vicinity of this resonance frequency the polarizability a of 
the element is not independent of frequency. If the element size is made 
equal to or less than one quarter wavelength at the smallest operating 
wavelength, it is found that a and hence ¢ in equation 2 is substantially 
constant for all longer wavelengths. 

Since lenses of this type will effect an equal amount of wave delay at all 
wavelengths which are long compared to the size and spacing of the objects, 
they can be designed to operate over any desired wavelength band. For 
large operating bandwidths, the stepping process? is to be avoided, since the 
step design is correct only at one particular wavelength. Such unstepped 
lenses are thicker, but the diffraction at the steps is eliminated and a som« 
what higher gain and superior pattern compared toa stepped lens is achieved 
By tilting the lens a small amount, energy reflected from it is prevented from 
entering the feed line and a good impedance match between the antenna and 
the line can be maintained over a large band of frequencies. 

Another way of looking at the wave delay produced by lattices of sma 
conductors is to consider them as capacitative elements which “load” free 


2 The lens of Fig. 1 has 12 steps. 
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space, just as parallel capacitors on a transmission line act as loading ele- 
nents to reduce the wave velocity. Consider a charged parallel plate air 
ondenser with its electric lines of force perpendicular te the plates. Its 
capacity can be increased either by the insertion of dielectric material or by 
the insertion of insulated conducting objects between the plates if the ob- 
‘ects have some length in the direction of the electric lines of force. This is 
because such objects will cause a rearrangement of the lines of force (with a 
consequent increase in their number) similar to that produced by the shift, 
due to an applied field, of the oppositely charged particles comprising the 
molecules of the dielectric material. The conducting elements in the lens 
may thus be considered either as portions of individual condensers, or as 
objects which, under the action of the applied field, act as dipoles and pro- 
duce a dielectric polarization, similar to that formed by the rearrangement of 
the charged particles comprising a non-polar dielectric. Either viewpoint 
leads to the delay mechanism observed in the focussing action of the arti- 


ficial dielectric lenses to be described. 
EXPERIMENTAL MODELS 


We turn now to experimental exemplifications of lenses built in accordance 
with the principles outlined above. 
(a) Sphere Array 

One of the simpler shapes of conducting elements to be tried was the 
sphere. Figure 2 is a sketch of an array of conducting spheres arranged 
approximately in the shape of a convex lens. The spheres are mounted on 
insulated supporting rods; the microwave feed horn and receiver are shown 
at the right. The focal length is /, the radius of the lens “‘aperture’’ is y, 
the maximum thickness is x and not only the spacings s; and s2 but also the 
size of the spheres are small compared to the wavelength. Rays A and B 
are of equal electrical length because ray A is slowed down or delayed in 
passing through the lens. Figure 3 is a photograph of the lens of Fig. 2; 
it also portrays a similar sphere array lens made of steel ball bearings sup- 
ported by sheets of polystyrene foam‘ which have holes drilled in them to 
accept the spheres. In both cases the balls are arranged in a symmetrical 
lattice. It will be shown below that the polarizability a of a conducting 

3 Polar dielectrics have arrangements of charged particles which are electric dipoles 
even before an external electric field is applied; the field tends to align these and the amount 
of polarization (and hence the dielectric constant) that they exhibit depends upon tempera 
ture, since collisions tend to destroy the alignment. Non-polar (or hetero-polar) molecules 
have no dipole moment until an electric field is applied; the polarization of such materials 
(and of the artificial dielectrics we,are discussing) is accordingly independent of tempera 
ture. See, for example, Debye, “Polar Molecules”, Chap. III. 

‘Styrofoam (Dow). Because of its low density (1 to 2 pounds per cu. ft.), its contribu 
tion to the wave delay is negligible (e, = 1.02 
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Fig. 3—Left: The sphere lens of Fig. 2. Right: Sphere lattice supported by foam 
sheets. t 
sphere of radius a is 47e,7*, so that, for static fields, the relative dielectri: Fig 
constant is, from (2), — 
€é- = €/e9 = 1+ 41Na', (3 S Toa 
: magi 
where .V = number of spheres per unit volume. dines 
For most dielectrics, the index of refraction is simply the square root of (b) | 
the relative dielectric constant. However, in the case of the sphere lens at Or 
microwaves eddy currents on the surface of the spheres prevent the magneti: sphe 
4 
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lines of forces from penetrating them and it will be seen later that this effect 
= causes the expected value of the index of refraction to be smaller than that 


® determined by the usual equation 








Fig. 4—Lattice of conducting disks arranged to form a plano-convex lens. Polystyrene 


foam sheets support the disks 


To avoid this effect the elements should be shaped so as not to alter the 
magnetic lines of force. This can be done by making their dimension in the 
direction of propagation of the impressed waves negligibly small. 
(b) Disk Array 

One accordingly arrives at the lens design shown in Fig. 4 in which the 


spheres are replaced by copper foil disks lying in planes parallel to the 
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impressed & and H vectors. The supporting sheets are again polyst 
foam. The foil disks have negligible thickness so that the magnetic | ines 
are unaffected as shown in Fig. 5. Equation (4) is therevore valid ev: 
radio frequencies and the index of refraction of this material is obtaii 
; pate oh 16 €9 a° ‘ ; 
from (4) and (2) with a equal to — ~ — , as shown in the last section. 
* ] 

(c) Strip Array 

Both the sphere and disk type lenses have the advantage that they 
perform equally well on horizontally or vertically polarized waves. Ii 
lens is required to focus only one type of wave polarization the disks can ly 
replaced by thin, flat, conducting strips extending in length in the direct 





Fig. 5—Disturbance of the magnetic field is avoided by the use of disks instead 
spheres. 
of the magnetic vector of the applied field. A simple method of constructing 
such a lens is shown in Fig. 6. Slabs of dielectric foam are slotted to a 
depth equal to the strip width and each slab is marked with the protile 
contour necessary to produce the final plano-convex lens. The metal 
strips are then cut to the length indicated on the profile and inserted in the 
slots. With the strips in place the unit slabs are stacked on top of one 
another and held in a mounting frame to form the complete lens. Figure 7 
shows one slab of a 10-foot strip type lens being assembled and Fig. 8 shows 
a six-foot square lens half assembled. Figure 9 shows directional patterns 
of a 3-foot diameter lens of this type made up of } inch x .005 inch strips 
spaced }” in the slabs and the slabs 13” thick. The three patterns were 
taken over a 12% band of frequencies with the lens purposely illuminated b 
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low directivity feed at the focal point. This produced an almost uniform 
lumination across the aperture at all three wavelengths so that the side 
obes were not too well suppressed. However, the deep minima in all three 
directional patterns indicate the relative absence of curvature of the emerg- 
ing phase front; this signifies that the strip dielectric has a negligible varia- 
tion of refractive index over the indicated wavelength band. 
d) Sprayed Sheet Lenses 
The disk lens or strip lens can be constructed in the manner indicated in 
a Fig. 10, which shows two lenses made by spraying conducting paint on thin 
dielectric sheets through masks. This results in round dot or square dot 


CONVEX f- SLOTS f- METAL STRIPS CUT 
CONTOUR /\ TO PROPER LENGTH 








CIRCULAR — 
OUTLINE 





Fig. 6—Construction of a delay lens employing metallic strips as the delay elements 
£ » “ 5 


ting 

to a patterns on the sheets and the size circle used on each sheet determines the 
ofile three-dimensional contour when the sheets are stacked up to form the lens 
etal Those in the photograph are spaced by wooden spacers as shown in the sketch 
the of Fig. 11; for large lenses it would probably be preferable to cement the 
one sheets to polystyrene foam spacers, thereby making a solid foam lens. 
re 7 Because of the small size of the elements, the lens on the left in Fig. 10, was 
— effective at wavelengths as short as 1.25 cm, in addition to longer wave- 
rns lengths (3, 7 and 10 cm). The lens of Fig. 12 was made by metal spraying 
rips metallic tin directly on circular foam slabs through masks, and the foam 
vere disks then cemented together. 





Strip lenses can also be constructed in this way. The lens in Fig. 13, 
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Fig. 8—A six-foot square strip lens half assembled. 
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Fig. 9—Directional patterns of a 3-foot diameter lens at several frequencies. 
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Fig. 10—Lenses constructed by spraying conducting paint on cellophane sheets 
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was made by affixing copper foil strips to cellophane sheets and stacking th 
sheets with no spacing other than the sheets themselves between them 
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is extremely close spacing and the staggered arrangement of strips cor- 


spond to a heavy capacitative loading and introduced so much delay per 


nit length that the measured effective dielectric constant of this lens proved 
to be 225 (i.e. the index of refraction was 15). Because of such a high 


dielectric constant the reflection losses at the surface of this lens are high. 





Fig. 12—Lens formed by metal-spraying tin directly onto polystyrene foam sheets 


so that, for efficient operation, tapered or quarter wave matching sections on 
each surface would be necessary. 
(e) Frequency Sensitive Delay Lensus 

When the conducting elements approach a half wavelength in their length 
parallel to the electric vector, resonance effects occur and the artificial 
dielectric behaves like an ordinary dielectric near its region of anomalous 
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dispersion.” The index of refraction of an artificial dielectric usin, 

















metallic elements would increase rapidly as the wavelength approa: 
13” until, at \ = 1)” the dielectric would be opaque. At still lower 
lengths, the material would appear to have an index of refraction less 


one. 





Fig. 13—Closely spaced strip construction comprising copper foil strips affixed t 
cellophane sheets. Juxtaposition of the sheets yielded an effective index of refraction 


15 


Elements, such as rods, which are \/ 2 long, have a very broad resonance 
band, and the region of anomalous dispersion in a dielectric utilizing suc! 


5 Anomalous dispersion occurs in optical substances when the frequency of the inciden' 
light wave approximates one of the vibrational resonance frequencies of the molecule. 01 
the long wavelength side of this resonance region the index of refraction is greater tha 
one and increases rapidly as the resonance wavelength is approached. Dispersion, whi: 
is the change of index of refraction with frequency, ‘s therefore very high, but it is the 
“normal” type of dispersion. At resonance, the abscrption of the wave is high and th: 
substance becomes opaque. At still shorter wavelengths, the resonance phenomer 
acts to make the index of refraction less than its long wavelength value, often less than on 
and the index again varies rapidly with wavelength, but because this is an “abnormal” 
type of dispersion, it is called the region of anomalous dispersion. 
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ing 3 rods is very large, i.e., the dispersion is not very great. If it is desired to 
Oar he ave a highly dispersive material, this band can be considerably reduced 
; Seed 
Wa ve p oby the process of tilting the rods so that they are more nearly perpendicular 
3S + to the electric vector. They t-ereby become “loosely coupled”? to the 
; incident wave and acquire a higher (0. Some unsymmetrical arrangement 
such as that shown in Fig. 14 is necessary to insure that the radiation 
damping of the elements is actually reduced, since a uniform tilt of all the 
elements would allow the array to radiate, unhindered, a wave polarized 
parallel to the elements. Measurements of the index of refraction of a 
dielectric made up of successive arrays of rods arranged as in Fig. 14 are 
a E 
/2 eon © 
a eel ee | 
tik te 
Fig. 14—Array of metal rods to produce a narrow dispersion band 
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Fig. 15—Measured index of refraction of a metallic dielectric comprising the tilted 
ced rods of Fig. 14 
‘tion 


given in Fig. 15, and one observes a marked similarity to the behavior of the 
index of refraction of dielectrics in the region of their anomalous dispersion. 


peated Because of the rapid change of with wavelength, such material may be 


y slic? - - . . - 
a useful as a means of separating narrow band radio channels by the use of 
neident prisms or by lenses having several feed horns located at the optimum focal 
le. On points for the various frequency bands involved. 
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met For radio relay applications, there are three electrical characteristics 
an one, which antennas should possess. The first is high gain (effective area), as 
orma 


this will reduce the path loss and accordingly the requirements on trans- 
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mitter power. The second is good directional qualities so as to min 
interference from outside sources and also interferences between adj 


antennas. The third is a good impedance match so that reflections bet \ e 
the antenna and the repeater equipment will not distort the transm)'t 
signal. These characteristics should preferably be attainable withou 
imposition of severe mechanical or constructional requirements. 

The shielded lens type of antenna (lens in the aperture of a horn) |i 
itself well to repeater work because of its moderate tolerance requirements 
its good directional properties associated with the shielding, and the desirab| 
impedance characteristic obtainable by tilting the lens. The delay lens 


offers the additional advantage of broad band performance with the 
sequent possibility of operating on several bands widely separated in wave 
length. In this section, construction details and performance of a 6-fo0t 
square aperture strip type delay lens will be discussed. 
(a) Destgn of the Artificial Dielectric 

The operating frequency band envisioned for this antenna was 3700 to 
4200 megacycles (A = 7.15 to 8.1 cm), and to keep the element lengt! 
sufficiently well removed from the half-wave resonant length a value of {” 
for the strip width was chosen. The index of refraction of solid polystyren« 
is approximately 1.61 and this introduces a reflection loss (mismatch loss 
at each surface of 0.225 decibels. To reduce this loss to 0.18 db. the arti 
ficial dielectric was designed to have an n of 1.5 as this still did not impart 
too great a thickness to the lens. A combination of strip spacings whic! 
vields an » of 1.5 was found to be 3” in the horizontal direction and 1,°,” 
center to center spacing in the vertical direction as shown in Fig. 16. The 
construction method of Fig. 6 was used which involved inserting .002” 
copper strip into slots cut in foam sheets. 
(b) Lens Design 

Several lens shapes were possible: (1) bi-convex, (2) plano-convex with 
the flat side toward the feed, and (3) plano-convex with the curved side 
toward the feed. For a given thickness and therefore weight of lens, the 
third possibility produces the shortest over-all structure of lens plus horn 
feed, and it was accordingly selected. The curved side is then a hyperboloid 
of revolution as shown in Fig. 17 and for the chosen focal length of 69”, the 
profile, as calculated from the equation shown for » = 1.5, reaches a maxi- 
mum thickness of 16”. ‘To eliminate the reflection from the lens into the 
feed, a lens tilt could have been employed. It was found that a quarter 
wave offset of one half of the lens relative to the other half could accomplish 
this same purpose, because reflected rays from one half of the lens then 
undergo a one half wavelength longer path in returning to the feed and the 
reflections from the two halves cancel. As this process, however, is com 
pletely effective only at one frequency, the final lens design employed both a 
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Construction of a six-foot square shielded delay lens for 


Ihe wooden horn has a metal interior surface 
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tilt in one plane and a quarter wave offset in the other plane. To \ 
most efficiently the space afforded the antenna on the top of the relay t 
the lens aperture was made square. Since an unstepped lens has, | 
nature, a circular aperture, the four corners of the horn aperture mu 
filled in with lens material as sketched in Fig. 18. The step heig 
designed for midband wavelength and in the present case follows the « 
tion of Fig. 17 with the focal length reduced from the value used fo: 
main lens section by AX, (nm — 1) where K is an integer. That integer A 
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Fig. 18—Filling in the four corners of a square horn aperture with lens material 


selected which brings the step profile nearest the horn corners as indicated 
in Fig. 18. 
(c) Horn Shield 

Lenses can be energized by means of a small feed horn placed at the focal 
point, but to obtain the best directional properties a full metallic shield 
extending from the wave guide feed up to the sides of the lens should be 
employed. The use of foam slabs in the construction of the delay lens 
prevents this shield from extending completely up to the lens as indicated 
by the extension of A-.1,; (shown dotted) in Fig. 16. Optical formulas 








; indici 
a spre: 
W ill b 
m cond 
= assoc 
§ d) / 
: Th 
The 
e area 
area. 
expel 
40 
| 
i 39 
=) 
Ww 
. © 
% U 
w38 
(a) 
z 
Zz 
a7 
rt) 
r Ww 
7 re 
; ye 
: +36 
- Oo 
: a 
, o 
, << 
H 35 
$ 
5 34 
: 
s The 
§ 
b indi 
H of re 
f freq 
) freq 
i (1 < 
m3 
f . 
* gain 
; 4 
> hor 
; Fig. 
; 
by 
| tole 
| 
3 
z 1 


i bona actor 








Bord 


a Base 


Ait SOAR eee ee Tk 


ee ee 


1 Ahn Nd a a 


SUF il asic nse 


WETALLIC DELAY LENSES 5 


indicate, however, that for such large apertures, the amount of diffraction 
spreading of the waves outside of the pyramid formed by the dotted lines) 
will be small, and it is only necessary that the sides 4.1,, 4,8 and BC all be 
conducting to insure good back lobe suppression and other desired properties 
associated with a horn shield. 
d) Performance 

The gain of this antenna over an isotropic radiator is plotted in Fig. 19. 
Che top curve is the theoretical gain of a uniform current sheet of the same 
area (6’ x 6’), the lower curve the gain of a 6’ x 6’ area having 6° effective 
area. The points, which fall approximately on the lower curve, are the 


experimentally measured gain values of this antenna at the frequencies 
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Fig. 19—-Measured gain characteristics of the six-foot square shielded lens of Fig. 16 
The lower line indicates 60°, effective area and the circles are experimental points 


indicated. The constant percentage effective area indicates that the index 
of refraction of the lens material remains quite constant over the indicated 
frequency band. In contrast to this, the 10’ x 10’ metal plate lenses! 
(n < 1) exhibit, at the band edges, a falling off of 1} decibels from midband 
gain for a 10°% wavelength band. 

The magnetic plane pattern of the lens when fed by a 6-inch square feed 
horn is shown in Fig. 20. The remarkable symmetry of the minor lobes in 
Fig. 20 and also in Fig. 9 shows that the phase fronts of the waves radiated 
by these antennas are very accurately flat. This result emphasizes the 


tolerance advantages of the lens over the reflector. It is believed that the 


1 Loc. cit. 
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further improvement in pattern of this lens over previous wave guide 
metal lenses is attributable to the absence of steps, as they tend to intr 
diffraction effects. The symmetry of the pattern also indicates a high di 
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rig. 20-— Directional pattern of the lens of Figs. 8 and 16 when fed with a six 
square electromagnetic horn 
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Fig. 21— Directional pattern of the same lens when enclosed with a full horn shield 


of homogeneity of the dielectric (the lens is 16 inches thick); this is a prop 
erty not always shared by ordinary dielectrics such as polystyrene. 


When a full horn shield is used, the illumination across the aperture has 


very strong taper (somewhat stronger than a cosine taper because of 1 
wide flare angle of the horn). This results in a very effective suppression 
the close-in side lobes as shown in Fig. 21. In the vertical plane, the sic 
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bes are not very well suppressed because the illumination is only slightly 
tapered. However, for repeater work, lobes in the vertical plane are not as 
objectionable as lobes in the horizontal plane because interfering signals 
generally originate from other stations lying in the horizon plane. 

The impedance match of the shielded lens antenna is affected by the 
discontinuity of the wave guide at the expanding horn throat and by what- 
ever energy is reflected back into the feed line from the lens. By tuning 


means, the throat mismatch can be held to less than 0.2 db SWR (1.02 


V.S.W.R.) over the 500 me band and by a combination of lens tilt and 
quarter-wave step in the lens the SWR due to the energy reflected from the 


lens was also held to less than 0.2 db over the band. 
PART II—THEORETICAL CONSIDERATIONS 

We turn now to a consideration of the electromagnetic theory underlying 
the operation of artificial dielectrics. If the polarizability of the individual 
conducting element employed is known, equation (2) will permit a calcula- 
tion of the effective dielectric constant. Before proceeding to this, however, 
a brief review of dipole moment, dipole potential, and dielectric polarization 
will be given (MKS units). 


DiroLE MOMENT, POTENTIAL AND POLARIZATION 
Two charges, +y and —g, displaced a small distance from one another, 
constitute an electric dipole. If the vector joining them is called ds, the 
dipole moment is defined as 
q ds. 
The potential 1’, at any point, due to a point-charge gq is defined as 
y= q 4rrer, 


where r is the distance from the point in question to the charge y. 
potential V due to a dipole of moment m is 


m 
COS A. 


trrer* 
where @ is the angle between r and m. 
A conducting object becomes polarized when placed in an electric tield. 
Its dipole moment m depends upon the field strength E and upon its own 
polarizability a: 


m = aE. &) 


If there are V elements per unit volume, the polarization P of the artificial 
dielectric is 


P = NaE., 
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But P is related to the displacement vector D and the dielectric consta 


follows: 
D cE = o£ + RP, 


so that € = e9 + Va which is equation (1). A knowledge of @ thus per 
a determination of e«. We obtain @ from (8) by first finding the di 
moment m of the particular shape of element when immersed in a uni 


held E. 


CALCULATION OF DIELECTRIC CONSTANTS OF ARTIFICIAL DIELECTRICS 


(1) Conducting Sphere 
Consider a perfectly conducting sphere immersed in an originally unil 
field of potent al 
| Ey = Er cos 6. 1] 
The free charges on the sphere are displaced by the applied tield and 
thereby becomes a dipole whose moment m we wish to determine. Thi 
external potential field is the sum of the applied potential and the dipol 
potential, and from (7) and (11) we have 
, m cos @ 
V ow Er cos 6 4 — l 
firey r” 
The internal field is zero because the sphere is conducting. At a boundary 
between two dielectrics, there is the requirement! 


V awtatas wos V inside: (1 3 
Equation (13) gives, at r = a (the radius of the sphere), 
; m cos @ 
Ea cos 6 + = =x (); (14 
418) a” 
or 
m = 4re, Ea, (15 


the dipole moment of the sphere. From (8) we see that the polarizability of 


the conducting sphere is accordingly 47re92°, from which equation (3) follows 
(2) Magnetic Effects of a conducting sphere array 

The above calculations on a conducting sphere assume an electrostati 
field. At microwaves, the rapidly varying fields induce eddy currents 01 
the surface of the sphere which prevent the magnetic lines of force fron 


penetrating the sphere. The magnetic lines are perturbed as shown in 


6 The author is indebted to Dr. S. A. Schelkunoff for the polarizability formulas give: 
in this memorandum. 
7 Smythe, “Static and Dynamic Electricity”, McGraw-Hill, 1939, p. 19. 
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Fig. 5. Now a conducting sphere in an alternating magnetic field is equi 
valent to an oscillating magnetic dipole,* and we should observe an effective 


permeability for a sphere array.* The magnetic dipole tield is, however, 


ypposed to the inducing field, in other words, the dipole moment Is negative 


Smythe’ shows that the magnetic polarizability of a conducting sphere of 
radius a in a high frequency field is 
a = 20a’. 16 


Phe effective relative permeability of an array of .V spheres per unit 


volume is therefore 


up = 1 — 2eoN. 


Phe index of refraction was given in equation (+) as the square root of the 
dielectric constant. This is strictly true only when the permeability of the 
die’ectric is unity. We have seen above, however, that the sphere array at 
microwaves possesses an effective permeability given by (17), and there 


fore (4) is not valid. The correct expression for 7 is 
n= WV Ure, 18 
and the effective refractive index of a sphere array is accordingly 
n= V (1 — 2nrNa*)(1 + 427Na’), 19 


which is smaller than that given by (3) and (4). The disk and strip arrays, 
besides being lighter, avoid the diminishing effect on the refractive index 
caused by the perturbing of the magnetic lines. 
(3) Conducting Circular disk 

The determination of the dipole moment of a disk involves the use of 
ellipsoidal coordinates and will not be carried through here. Smythe’ gives 
an expression for the torque on a flat disk of radius a in a uniform field in 
Gaussian units. In M.K.S. units, his formula becomes 


2a" E sin 20 2a" E 2 sin 6 cos 8 
T = 4re - = tire = 
<7 <T 
l6a'e5,,.. ' 
: (E sin 6) || E cos 9}. 20 
J 


The first bracket represents the dipole moment, the second the field, and 
the product gives the torque. When the plane of the disks is parallel to & 
sin @ is one, and from m = aE, we have 


- 16€) a 1) 


> 


Be 


eo 


8’ T.S. E. Thomas, Wireless World, Dec. 1946, }. 322 
*L. Lewin, Jour. J. E. E., Part I, Jan. 1947, p. 65 
7 Loc. Cit. Eq. 14, p. 397. 

7 Ibid, eq. 7, p. 163 
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so that, from (2), 
é-=1+48Na 
(4) Strips 
The calculation of the dipole moment of a thin conducting strip as sed 
in the strip lens of Fig. 15 involves two dimensional elliptic coordinates 


will also be omitted here. Again, the torque is given in Smythe? for 
elliptic dielectric cylinder, from which we obtain, 


= TE) S 


4 


’ 


where s is the strip width, so that 


where » is the number of strips per sq. unit area looking erd on at the strips 
(5) Validity of the polarizability equations 

Equation (2), which expresses the dielectric constant to be expected from 
an array of V elements each having a polarizability of a, was derived (by 
(8), (9) and (10)) by assuming that the fied acting on an element, and 
tending to polarize it, was the impressed field E alone. This isa satisfactory 
assumption when the separation between the objects is so large that the 
elements themselves do not distort the field acting on the neighboring 
elements. Such is not the case when the value of e, exceeds 1.5 or there 
abouts. For the usually desired values of €, of 2 or 3 it is thus seen that the 
above formulas such as equations (22) and (24) will yield only qualitative 
results and that the exact spacings of the elements to produce a desired 
refractive index will have to be determined by experimental methods. 

For lattices having 3-dimensional symmetry, an improvement over 
equation (2), which takes into account not only the impressed field E but 
also the field due to the surrounding elements, is the so-called Clausius 
Mosotti equation: 


é = 5 Na 
etle 3e, 
This, along with a similar expression for the permeability [replacing (17)], 


would permit a fairly accurate determination from (18) of 7 for the conduct 
ing sphere array. 


* Thid, eq. 6, p- 97, 
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RESONANCE EFFECTS 

It was stated earlier that the size of the elements should be small relative 
to a half wavelength in order for the refractive index to be independent of 
frequency. A qualitative idea of this criterion can be obtained by an 
elementary analysis of forced oscillations of dipoles. It is known that a 
dielectric medium which is composed of elements that resonate under the 
action of an alternating electric field, such as atoms having bound electrons, 
will exhibit a dielectric constant which varies with frequency: 


where fo is the frequency of resonance of the element, f the frequency of the 
incident radiation and & a proportionality constant. Thus when / is small 
relative to fo, €, is practically independent of /. 

As a means of estimating the change in refractive index with frequency of 


metal delay lenses, we consider a specific example: Let » = 1.50 when the 
elements are \ 4 in length, i.e. when /*? = 1/5, then the last term of (26 


equals 1.25. Decreasing f by 20°; reduces n from 1.50 to 1.46. Thus the 
change in n from midband to the edges of a + 10°; band is about .02. 
From this, at 7 cm, the phase front, even for a lens 30)” thick, should remain 
plane to within + ’5 wavelength over this 26°; band of wavelengths. If 
the members had been made } wavelength long, the variation in n from the 
design frequency all the way down to D.C. would have amounted to only 


ons 
1.2%. 
SUMMARY 


A metallic dielectric is constructed by arraying conducting elements in a 
three-dimensional lattice structure. For electromagnetic waves whose 
wavelength is long compared to the size and spacing of the elements, this 
structure displays an effective dielectric constant and index of refraction 
which is sensibly constant over wide frequency bands. Lenses can be 
designed according to these principles which will focus microwaves and 
longer radio waves as a glass lens focusses light waves. Such lenses have 
the advantage of broad-band performance over the earlier waveguide type 
metal lenses and they retain the advantages of light weight over dielectri: 
lenses. As microwave antennas, they are superior to parabolic dish reflec- 
tors from the standpoint of warping and twisting tolerance, profile tolerance, 
directional properties and impedance match. By eliminating the steps in 


the lens, the directive patterns are made cleaner and an increase in absolute 


10 See for example, Joos, Theoretical Physics, Blackie & Son, Book 4. Chapte> 4 
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gain results. Because of the broad-band properties of the artificial diele: 
this improved gain can be maintained over a very wide band of frequenc es 
The lenses can be built to focus waves of any polarization and, if desi 
the dielectric can be designed to exhibit strong dispersion. Theoret 
calculations of the expected dielectric constant are in fairly good agreen 
with experiment for values less than 1.5; for higher values an accu: 
determination of the true value must be obtained experimentally. 
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A Non-reflecting Branching Filter for Microwaves 
By W. D. LEWIS and L. C. TILLOTSON 


Microwave branching filters are required as integral parts of multi-channe 
microwave radio relay systems. These filters must have characteristics which 
are difficult to attain if one attempts to extend familiar lower frequency tech 
niques to the microwave region. A novel network configuration, through which 
currently anticipated requirements can be met without excessive difficulty, is 
described in this paper 

In this configuration individual constant resistance channel dropping units are 
formed of appropriate assemblies of two hybrid circuits, two band reflection 
filters and two quarter wavelengths of line. An assembly of N channel dropping 
units in cascade then forms an N channel constant resistance branching network 

The mechanical and electrical characteristics of a practical five channel brancl 
ing filter of this type are described. Asa result of experience with this prototype 
filter it can be stated with some safety that these requirements can be fulfilled with 
a network of this type. Experimentally observed impedance, insertion loss and 
phase characteristics were fully satisfactory. In addition the circuit appears to 
be flexible enough both electrically and mechanically to fulfill the various types 
of systems needs which may be encountered at branch points or when channels 
must be added or interchanged 


INTRODUCTION 


RESENT plans for point-to-point communication by means of micro 

wave radio relay systems call for the operation of several radio channels 
between each pair of repeaters. A proposed frequency plan for the 40.0 m« 
common carrier band (3700 to 4200 mc) specifies channels 20 me wide spaced 
40 me center to center. A possible arrangment for a five-channel radio 
repeater station is illustrated in Fig. 1. This arrangement is calculated to 
utilize the available frequency space in an efficient and technically sound 
manner. 

If this channel disposition is to be achieved without a costly increase in the 
number of antennas and the size of the supporting towers, radio frequency 
branching networks must be provided which connect the individual trans- 
mitting or receiving circuits at each repeater point to a common antenna 
(Fig. 1). If this connection is to be made without excessive loss of power 
these branching devices must have adequate adjacent channel rejection, low 
ohmic loss, and good impedance match in the channel bands. An excellent 
impedance match is especially desirable if circuit disturbances resulting from 
echoes in the long waveguide lines which lead from the filter assemblies to 
the antennas are to be minimized (Fig. 1). 

Since the type of microwave radio repeater now planned! obtains most of 
its gain at intermediate frequencies, the IF amplifiers will reject all spurious 

'See H. T. Friis, “Microwave Repeater Research’’, to appear in the April 1948 issue 
Oo.S.T. J. 
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signals entering the receiver except those in the vic inity of the receiver i: 
bands. Because of this, suppression requirements on the branching fi) er 
except possibly in the vicinity of receiver image bands, are not severe. 
Problems connected with the design of suitable microwave bran: 
filters naturally differ considerably from previous filter problems. The 
crimination and band utilization requirements are readily met, but 
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Fig. 1—Possible five-channel radio repeater station schematic diagram. 


impedance control required is not easy to obtain by familiar filter techniques 
Not more than about 5°%, variation in impedance or about 0.5 db standing 
wave ratio in the channels could be tolerated. At lower frequencies channe! 
passing networks which can be connected in series or parallel to form a 
channel branching filter can be designed on the basis of lumped circuit theory 
and built of coils, condensers and resistances, but in the microwave region 
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VON-REFLECTING BRANCHING FILTER 


simple parallel or series connections and simple lumped circuit elements do 
not exist. Although in the interests of flexibility it would be desirable to 
add or substitute individual channels without affecting other channels ina 
branching filter, the convenient possibility of doing so at a high impedance 
level on vacuum tube grids is not yet available in the microwave region. 

A satisfactory two-channel waveguide branching filter has been constructed 
following partially ‘classical’ methods. This filter, designed for the New 
York-Boston experimental radio relay system, is composed of two channel- 
passing ¢ avity filters each connected to a common input line through one arm 
of an E plane Y junction, the waveguide analogue of a series connection 
(Fig. 2). This solution, although relatively simple where only two channels 
are required, becomes quite complex when more than two are involved, since 
in every channel the sum of the interactions of all the inactive filters on 
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Fig. 2--Branching filter for New York—Boston experimental radio relay system 


transmission through the active filter must be zero. It is evidently not easy 
to satisfy this condition, particularly since in doing so one must take accurate 
account of the change with frequency of the effective phase shift of all wave- 
guide connecting lines. And even if such a solution were found it would be 
valid for only one set of channels, so that the problem must be solved all over 
again for every change in channel arrangement. 

As a result of these difficulties and after a few attempts to overcome them, 
it became apparent that a more flexible method of microwave filter con- 
struction should be found. Constant resistance filters, which provide dis- 
crimination by absorbing or diverting the unwanted incident waves rather 
than by reflecting them, can be useful in any frequency range. In the 
microwave region, where the shortest connecting pipe may be many wave- 
lengths long, constant resistance devices are particularly helpful. Accord- 
ingly a constant resistance channel-dropping network was devised which 
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could extract one channel from the line, while permitting others to 


through it without disturbance. Several of these networks were then pla 


in cascade to make up the required filter. 


THE Hyprip CHANNEL-DrRopPING UNII 


After several possibilities were considered the constant resistance chan: 
dropping circuit illustrated schematically in Fig. 3 was selected 1 
circuit is made up of two hybrid junctions, two identical channel reflect 
filters tuned to the dropped channel, and two quarter wavelength sections 
line. 

In order to understand the operation of this circuit, the properties of 
hybrid circuit, Fig. 4(a), must first be understood. This circuit, which | 
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CIRCUIT FLEXIBILITY 


Possible constant impedance channel dropping filter. 


been embodied at voice frequencies as the hybrid coil and at microwaves as 
the hybrid junction (E-H plane T junction)’ can be represented schematicall) 
as in Fig. 4(b). Let us connect four transmission lines, A, B, C, D, eact 
terminated in its characteristic impedance to the four pairs of terminals of 
the hybrid circuit. Then if each of these lines is matched to the pair of ter 
minals to which it is connected, the following characteristics will result 
Power in transmission line C flowing towards the junction will divide equall) 
into lines A and B, flow away from the junction and be absorbed in loads A 
and B. None of this power will appear in line D or be reflected back into 

* Lumped element networks with properties similar to those of this circuit have bee: 
devised by Vos and Laurent, U. S. Patent 1,920,041, and Bobis, U. S. Patent 2,044,04 
\. G. Fox of these jaboratories has irdependently devised similar microwave circuits. 


}W. A. Tyrrell, Hybrid Circuits for Microwaves, Proc. I. R. E., Vol. 35, pp. 1294 
1306, November 1947. 
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NON-REFLECTING BRANCHING FILTER 


eC. Similarly, power in line D flowing towards the junction, will appear 

jually in A and B but not in C or back in D. If these characteristics hold, 

en by the principle*of reciprocity similar characteristics must hold for lines 
\ and B. If proper planes of reference are chosen this behavior can he 
lescribed in a slightly different, but equivalent, manner. If waves in both 
\ and B flow towards the junction the vector sum of the voltages of these 
times a constant (0.707) appears in C and the vector difference times the 
same constant appears in D, but nothing is reflected back into A or B. - An 
equivalent statement can be made if the waves start in C and D. 

With the properties of the hybrid in mind, and if it is assumed that the 
hybrids, the identical reflection filters and the quarter wave lines are perfect 
and free of ohmic loss, the operation of the circuit of Fig. 3 is easy to under- 
stand, and is as follows: A wave entering from arm C of the input hybrid is 
divided equally into the two arms A and B. None of the power in this wave 


V 
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~w =—\ 
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Fig. 4—a. Classical hybrid coil 
b. Hybrid junction schematic diagram 


is reflected back into the arm C or appears initially inarm D. The two equal 
components of the wave now travel along the lines which are connected to the 
arms A and B of the input hybrid. If the frequency lies outside the band of 
the reflection filters the waves travel through these filters and appear in phase 
in the arms A and B of the output hybrid. The vector sum of these two 
waves appears in the arm C of the output hybrid and has an amplitude equal 
to that of the original input wave. Consequently all energy in the input line 
incident on this network, except that lying in the band of the reflection filters, 
will pass through it to the output line. 

If now the frequency of the input wave lies within the band of the reflection 
filters, the two equal components of the wave traveling away from the input 
hybrid wil! be reflected at the filters and will travel back towards the input 
hybrid. One of these components must, however, travel twice through an 


extra quarter wavelength of line, and will therefore be reversed in phase; with 
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respect to the other component when it reappears at the input hybrid. 
two comporents will consequently combine in the fourth or difference arn 
the input hybrid to form a wave equal to the original input wave. 

The circuit of Fig. 3 is only one of a general class of hybrid filter circu: 
From one viewpoint these circuits resemble spectroscopes, and, from anot} 
ordinary lumped circuits. We could, for example, replace the band reflect 
circuits of Fig. 3 with any two identical four-terminal networks. This circ 
would still retain its constant resistance character. One of its branches 
would contain all energy reflected by the four terminal networks; the ot! 
would contain all energy passed by them. 

In particular if the two identical filters are of the channel passing ty). 
input waves within this channel will be transmitted by both filters and w 
combine at the output hybrid and appear inarm C. All of the other channe 
will be reflected by the filters and thus will appear in arm D of the input 
hybrid, provided that the assumed 90° phase shift holds over a band whi 
includes all of the channels. 

The particular configuration of Fig. 3 was chosen to minimize the effect oi 
practical limitations. The dropped channel width (20 mc) is only a sma 
fraction of its midband frequency (about 4000 mc). Consequently when band 
reflection filters are used in Fig. 3, the change with frequency of the nominall) 
quarter wave sections of guide does not seriously affect the performance of the 
filter and the impedance match of arm D of the hybrids needs to be good only 
in the vicinity of the dropped channel. 

The circuit shown in Fig. 3 is a constant resistance network which drops 
the channel corresponding to the reflection filter. Several in sequence as 
indicated in Fig. 5 constitute a constant resistance channel branching filter 
In the sections to follow we will give an account of the physical configuration 
and electrical performance of a branching network designed according to this 
pattern to meet the radio frequency requirements of a typical practical radio 
relay system containing five 20 mc radio frequency channels, spaced 80 m« 
center to center. 


THE WAVEGUIDE HyBRID 


In choosing a hybrid configuration which could be successfully used in the 
network of Fig. 5, both electrical and mechanical requirements were considered. 
Since frequency f, passes through 2n —1 hybrids it is evident that if accept 
able overall performance js to be obtained, the balance and impedance char 
acteristics of each hybrid must be excellent. A broad-band balance can be 
obtained with relative ease by attaching the ‘driven’ arms (A and B on 
Fig. 5) symmetrically. Fortunately the strict impedance requirement applies 
to only one of the two ‘driving’ arms (C and D), the other being required to 
transmit only a single channel. 
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Because of the reentrant nature of the circuit of Fig. 5 it is evidently desir- 
ble, if not essential, to employ a hybrid with a convenient mechanical layout 
if, for example, a familiar E-H plane junction type of hybrid were used in 
ombination with waveguide filters built in a straight piece of waveguide, 
twenty-four E or H plane right angle waveguide bends would be required in 
the construction of each six-channel branching network. To avoid this extra 
complication and expense a hybrid configuration with ‘driven’ output arms 

parallel to the well matched ‘driving’ input arm was sought. 

The hybrid configuration settled upon was the one shown in Fig. 6. Here 
the electrical and geometrical recuirements discussed above are met simul- 
taneously by connecting the driving arm C to the symmetrically located 
driven arms A and B through a smoothly tapered E plane Y junction. The 
taper is approximately one half wavelength long in the center of the 3700-4200 








f, OUTPUT f,, OUTPUT 


HYBRID JUNCTION VW “DUMMY LOAD 
t¢ ti REFLECTION FILTER j QUARTER - WAVE 
L FOR FREQUENCY f WAVEGUIDE LINE 


Fig. 5—N Channel branching filter 


mc band. The driving arm D is connected to arms A and B through a coaxial 
line. This line is coupled to waveguide D by means of a conventional probe. 
The center conductor traverses the Y junction space in such a way that it is 
normal to the electric vectors of guides A, B, and C, and is effectively coupled 
to A and B but not to C by means of a probe P fastened through it normally 
(Fig. 6). 


THE BAND REFLECTION FILTERS 


The ideal reflection filter for the circuit of Figs. 3 and 5 would reflect per- 
fectly within a certain band and pass perfectly outside of this band. However, 
the ratio of bandwidth to band spacing in a given branching filter (20 me to 
80 mc) is such that a sufficiently good approximation to this ideal can be 
obtained in theory if each reflection filter employs only three resonances, 
Fig. 7(a). These could be effectively series resonant circuits placed at quarter 
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wave intervals along the guide, properly distributed in impedance leve 
all tuned close to the center frequency of the channel to be extracted, Fig 
The practical question was to tind how to obtain these resonances in an e 
constructed and adequately adjustable form. 

Early experiments indicated that a probe inserted in the broad side of (hy 
guide far enough so that its end formed an appreciable capacitance with ‘he 
opposite side could be made to resonate in a series resonant fashion. Im; 
ance levels available through this means of coupling were, however, far lower 


than required. Accordingly an alternate method in which the probe is 
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Fig. 6—Hybrid junction. 


inserted in the narrow side of the guide was selected (Fig. 8). The probe is 
made long enough to approach the opposite narrow wall of the guide and is 
tipped by a capacitative disk. The capacity of the disk and consequently the 
resonant frequency of the circuit is adjusted by means of a screw in the wal! 
just opposite the disk. Since this rod is inserted perpendicular to the narrow 
guide wall it is normally uncoupled to the principal mode in the guide. An 
adjustable coupling is achieved by inserting a screw in the broad side of th 
guide just above the probe. Insertion of this screw disturbs the symmetry oi 
the field and couples the rod to the guide. Increasing the screw insertio1 
increases the coupling and consequently varies the impedance level of the 
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equivalent series rcsonant circuit continucesly from infinity down to any value 
vithin the range required. 

It should be pointed out that impedance and frequency adjustments are not 
in practice completely independent, but do at least permit the realization of 
any required value. Furthermore, mutual coupling between the probes inter 
feres with the exact realization of the circuits of Fig. 7. This coupling does 
not interfere, however, with the realization of satisfactory filter character- 


istics. 
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Fig. 7— Band reflection filter 
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Fig. 8—Waveguide band reflection filter. 


When the reflection filter of Fig. 7 was embodied in waveguide form through 
use of the series resonant circuits just described, the configuration illustrated 
in Fig. 8 resulted. It was found that, with the exception of a slight change in 
disk to guide wall spacing, a single configuration could be tuned to any of the 
desired channels. 


ELECTRICAL PERFORMANCE 


Individual components as described above were constructed, adjusted and 
assembled to form a five-channel branching network (Fig. 9). The electrical] 
performance of this network was trimmed systematically, then was measured 
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Input standing wave ratio of hybrid branching filter 
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With all outputs ter- 


minated the standing wave ratio observed at the input line was under 0.6 db in 
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allchannels. This quantity is plotted in Fig. 10. The insertion loss m: 














































































between the input line and the various output lines varied from about |)5 spt 
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| channel dropping unit, it appears logical to assume that the progressive 4 
increase in loss is due to passage of the higher frequencies through the lower 
i frequency circuit components. 
| The measured discrimination against other channels and image responses | 
was 20 db or more. This modest amount of selectivity is sufficient since the | 
primary function of the present filter is branching, and hence the amount of 
ie 
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NON-REFLECTING BRANCHING FILTER 


scrimination needed is only that required to prevent a significant energy loss 
the channel bands. If crosstalk considerations indicate that more r-f dis 
$ crimination is required, this can be obtained by placing auxiliary filters in the 


branch arms. This can be done without complication since the branching 


* filter is a constant resistance device. 
[he delay distortion was examined but was found to be less than the errors 


of measurement, i.e. 2 millimicroseconds. 
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A Note on a Parallel-Tuned Transformer Design 
By V. C. RIDEOUT 


An analysis of the parallel-tuned transformer used in radio-frequency amp 
fiers has been made for a slightly over-coupled case. The resulting desig 
formulas are simple and practical. 

Two cases are discussed: (a) the so-called matched transformer, with resist 
ance loading on each side; (b) a transformer with loading on one side onl 
which has the same pass-band and phase characteristics as the matched trans 
former, but gives 3 db more gain when used as an interstage. 

A special arrangement of (a) where the matched transformer design is use: 
with one resistor removed, giving a transformer with a considerably doubl. 
humped pass-band characteristic and about 6 db more gain, is also discussed 


INTRODUCTION 
HE parallel-tuned transformer has been used in_ radio-frequency 
amplifiers for many years.' In an excellent paper,’ Christopher based 
design formulas on the principles of the broad band filter. In other publica- 
tions, simple circuit analyses have been used and design formulas have been 


based on the assumption of a small ratio of bandwidth to mid-frequency 


which often fails to be adequate when the wide bands required for modern 
television and multiplexing services are encountered. 

A transitionally flat transformer design may be obtained by setting the 
first three derivatives of the absolute value of the transfer impedance, wit! 
respect to frequency, equal to zero. The resulting design formulas ar¢ 
somewhat unwieldy. 

The transformer design to be described here is based upon two simple 
circuit conditions* applied to the fundamental case of a_parallel-tuned 
transformer with resistance loading on each side: 

I. Both sides of the transformer are tuned to the same frequency. 

II. The transmission loss is zero at the tune frequency. (This condition 
is responsible for the term ‘“‘matched transformer’? used to describe this 
case). 

The resulting transformer has a slightly double-humped characteristi 
with less than 0.005 db dip for a coupling coefficient of 0.5. Because of the 
slight overcoupling this transformer design gives a little more gain and 
bandwidth than the critically coupled case. Its main advantage lies in the 
fact that simpler design formulas can be used. 

'H. T. Friis and A. G. Jensen, High Frequency Amplifiers, B. S. T. J., Vol. II, Apri 
1924, pp. 181-205. 

2A. J. Christopher, Transformer Coupling Circuits for High-Frequency Amplifier 


B.S. T. J., Vol. XI, Oct. 1932, pp. 608-621. 
’ This analysis is based on an old unpublished report by H. T. Friis. 
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DERIVATION OF DESIGN FORMULAS 


or Symmetrically Loaded Transformer. 


Matched, 
In the circuit of Fig. 1: 

= 1(R: + 1/joC:) + I(— 1/400) 
Ti(— 1/420) + Io(jwls + 1/foCs) + 13(—joM) 
To(— jwM) + Is(jwle + 1/jwC2) + I4(— 1/jeC2) 


= I3(— 1/jwCe2) + I4(Re + 1/70C2) 
Ef If 
eC YO WAM. 


Ss R 
2 


= 2< 
M 
=KVL, Lo 
9 symmetrically loaded transformer 
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Cy 1 





E,(~v 








Fig. 1.—Parallel tuned, matched, « 


and I; among these equations gives: 


Eliminating J,, J» 
a 2 + <2 ; 2 4 s? 
Ly Re + Iai 2 | (1 7 hile) ck, 1 CR) | 


Eiji = - M M 


; ld >) Eee : Lo 
ae (| ( —_ 7 ) cic RR | + o| a (LiCi + L2C2) 
; (a Fs a) = aa 
M w M 


If the coefficient of mutual coupling is k, and the first and second induct 
ince and capacitance are resonant at w: and w, radians per second, respe: 
tively, we can substitute in (2) the expressions, 

/ p 2r 2;,. . 
M = ky L,L.,L:= 1 a3 1,L2= 1 wale, = F, R, 
This gives the general expression, 


.se R, VC, Co f @) l 
E,/E, = | tis 
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The application of circuit condition I, (w: = we = wo) to (4) gives, 


dat Ri wor/C, Co ( | ) 
E/E : e 
; k Ww C, R, ; a) 9 Rs 

— (1 — k°) t . 
WwW Wy ¢ 1 R, Wy C 2 Rs 
) i 9 ) | = k? Wi 

—ji(—) a —#) —-—(2+ 73 ) . 
. (2) | = ( wy C1 C2 Ry Re ’ w | 


Circuit condition II (zero transmission loss at w = ao) is satisti 


E./ Ey = 2VR,/R,. Substituting this condition in (5) and solvii 


k gives: 


waCiCoR Rs + 1 + (ani Ri — woC2 Ro) 
(wo CR)? t (wp Co Ry)? T (we C1 Co R; Re)? + | 


Since k must be real, 
wo Ri = ayCoRe 
and from (6) and (7) 
wo iRi = aploRo = Yl — RR & 


From (5) and (8) the transmission formula for this transformer is, 


7 fae , 2k . P 
A vi-t Gs} (1 : (5) it~ *) 
Es EF VGalvi-e - 


9 
((Je--fe+n2)] 
Wh Wy WwW 
: ky = FE ,¢ ; 
or, Es os R, ( , wnere, 
OS en ee mts a Fo € ) a ~ 1 
Eo k? Ci/ Ul — FP Wo : 
(10) 


+{(¢) (1 — k*) ——(2+#) +2 | 5 
. \@o Wo Ww ) 


js nite wy) (1 — k°) — (w/an)"(2 +R) +1 4 
= arcte : (11 
2k (w/wo)?(1 — k*) — w/a 
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e transmission loss, <, defined as the ratio of the maximum available 
er from the generator to the output power may be obtained from (10), 


| ae A 74] E 
= l (1 — k) 
iR, | Es k? [1 - 
a 


12 
“ 
L l oe . ; ° - WwW 2 n hb nl =i 
s Ww Ww 
the loss in decibels equals 10 log &. 
[he transmission delay, 6, may be obtained from (11 
w \! ; a Sua \2 
( ) (1 — ky)? - es) (1 — ky)? + 4k —14 =) 
do ae Wi WwW \W 
dw A (13) 
Vi- A 
2k ‘& 


seconds. 


The input impedance seen from the generator terminals is, 


Phe formulas (12), (13) and (14) give the important characteristics of this 


transformer. Table I gives expressions for these characteristics in the 


general case, and for k = 0.5, for several picked frequencies. Of these 

frequencies wo and wo/ WY 1 — k? are the frequences of zero transmission loss, 
4 a . . / 

wy Y1 — kis the geometric midband, and wo/Y/ 1 + k and wo/YV 1 — k 


are the cut-off frequencies of an infinite filter made of identical transformers 
of this type. 

The input impedance, transmission loss, and the transmission delay are 
plotted against w/w) for a matched transformer with & = 0.5 in Fig. 2. From 
these curves, and from Table I the following important characteristics of 
this transformer will be noted: 

(1) The pass-band is approximately symmetrical with frequency, rather 
than with the logarithm of frequency, as in circuits which have a low-pass 


analog. 
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Fig. 2.—Loss, delay, and input standing-wave ratio versus relative frequency for the 


matched transformer for k = 0.5. 


Loss and delay curves also apply to the 


case of the mis-matched transformer for k = 0.632 


(2) The delay curve is also quite symmetrical. 


} 


Circuits which can be 


derived from low-pass forms have more delay on the low-frequency side, 


where loss increases faster. The symmetry of tne delay curve makes phase 


compensation easier. 
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(3) The input standing-wave ratio passes through zero on either s 
mid-band, and although it has a 0.6 db dip at mid-band for the « 
k = ().5, it is under one db over a larger frequency range than in the t 
tionally coupled case. 

(4) The bandwidth between points one db down on the loss curve is ver 
nearly equal to the bandwidth between the cut-off frequencies. 
bandwidth is simply related to the tune frequency fo and the geometr 
midband frequency /,, as shown below. 


Af =fi/V1 —k—-fo/VI +k 


a 


a 


= fok/V/1 — F 
> fak/ V1 — Rk? 


Thus, given required values of /,,, Af, C: (or R:), and C2 (or Re), the matche 
transformer can be designed by the use of formulas (3), (8) and (15). Thes: 
formulas are summarized at the end of this paper. 

A transformer of different phase and loss characteristics results il 
loading resistance is removed from one side of a matched transformer 
Equation (5), if Rp = *, becomes, 


E Jt = § k + ‘ 
, = VEEB /Ef A (1—-(sfa—# ) 
Fy k? { C, V | -— Wy 


A 
0 
_ ((°.) Q<)~22 43 
wy Wo w 
At w = wo, 
E,/Ey |? = C2/ 17 
while from (10) the corresponding ratio for the matched transformer is, 
E,/ Ez |? = 4€C2/C, 18 


Thus this mis-matched transformer has 6 db more gain at the tune fre 


quency than the corresponding matched transformer, whether used as a: 


interstage, an output or an input transformer. 
The transmission loss referred to the loss at w = wo Is, 


5: ( - E, 2 | w 2 . 2 
x = rt = —<_ l = k*) 
¢ "9 Fs k : &, \ 


2 \3 2 
(*.) (i + #) 
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Fig. 3.—Loss and delay versus relative frequency for the mis-matched transformer ol 
tained by making R2 in Fig. 1 infinite, for k = 0.5 


The transmission delay 6’ in this case is, 


V1 — k 


6’ = 
Wy R3 
WwW : . : 
) 2 4 w@ ? - 
-ya—ry — ya -#) + Gh -1)+{(- 
Wy Wy, @ 
er seconds. 


(20) 
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The loss is equal at four points, 



























Ww = Wo, w wo V1 — k?, w = wo V1 — k, Ww = Wo V1 


The curves of Fig. 3 show the loss and delay for this mis-matched 
former for the case where k = 0.5. The double-humped band-pass c! ara 


ORRIN 


teristic of the mis-matched transformer can be compensated for in an ; 
fier by de-tuning slightly to equalize the humps and using a single-tuned | 
circuit in another interstage in which, 


% 
QO = wmCR & V1/k? — 3/4 — 3/32 n i 
The above combination of a mis-matched transformer and a single-tun 
transformer, in successive interstages, will give approximately 6 db more ey 


gain than if two matched transformers were used. The actual gain obtaine lies 





will depend upon the ratio of the input and output capacitances. This gai 
advantage carries with it the penalty of increased sensitivity to tube capa 


on eee 


ss mw! 
' ’ 


=KVLUYL 





ERS IE EYES 
— 


4 
a Se 
Pig 








Pe ie ee ae aca, 


£ san 


n> 


Fig. 4.—Mismatched or unsymmetrically loaded transformer. 


itance variations. However, the mis-matched transformer may be used 
to advantage in the first interstage of an amplifier to minimize the effect pin 
of the second tube on the signal-to-noise ratio. It may also be used in the ® bel 
last interstage to offset any lack of power-handling capacity in the second- 
last tube. 3 


| (b) Mismatched or Unsymmetrically Loaded Transformer 

A transformer with only one loading resistor, as in Fig. 4, can always bi 
designed to have the same bandpass and phase characteristics as the matched 
transformer‘ shown in Fig. 1. No power transfer will occur, but the rati 


3 

- - . . g 
of output voltage to the square root of the available input power can be 4 
determined, and set equal to this ratio for the matched transformer, excep! | 7 


‘ This important principle has been described by S. Darlington in a paper: Synthesis 0! 
Reactance 4-Poles, Journal of Mathematics and Physics, Vol. XVIII, Sept. 1939, pp. 257- | 7 
353. Independently, this principle was demonstrated to the author by W. J. Albershein 
in an unpublished memorandum dealing with the transitionally flat transformer. 
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For the matched case this ratio is, from equation 


(: 1 
Wo 


1 


Los) pemele 1 — k C: 
—j\— - — 
Wo 2k? | CLR, 
. (¢) (2+ kh)V1i — kB Ce 
in 2k Ci 


“i 


Vi-k#./C 


2 (9 
R, 


Vc 


2k? 


For the mismatched transformer the corresponding ratio can be obtained 


from equation (4) by letting R: approach infinity, giving, 


Ev/4R; _|1 2 ,/ C2 w 
| te? 2k’ wi ) CiRi wie 
1 — k” C2 w 1—k” / 
, 7 | To —j —WVRiC iC; 
2k | Cy Ri J wiws 2k ens 
, , /pir 7 be, [nipping 
; 12 ) /RiCiC2 . Wy We 2 
+jo in pon V RiCiC2 _ - Me meets (23) 
W1 wo 2k’ w 2k 


The two transformers are to have the same phase characteristics, and the 


same band-pass characteristics except for a constant factor N where, 


N= V 


EY 


/4R; 


Vy EP 


(24) 


The terms involving voltages in (22) and (23) can be eliminated by com- 


bining with (24). 
below. 


1— & 
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2+k° 
> VRiCiCe = 


Equating the coefficients of like powers of w gives (25), 
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The last three equations were simplified by the use of equation (8) 






Equations (25) yield the tive relations below: 






RiCi = R:C,/2 
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SUMMARY OF DESIGN FORMULAS 


ah; 


A transformer must usually meet certain requirements as to bandwidt 
and mid-frequency. Loss curves for various values of k, plotted agains: 
relative frequency may be used. To a very close approximation the g 
metric mean frequency fm = wm/27, may be used for the mid-band fr 
quency, and Af, the bandwidth between cut-off points, may be used for 
bandwidth between one db points. Then, 


4 9 
Wm = wo/YV1 — k? é 


And, from (15), 





bins os i Fs 
~ V1 + (Af/fo)? & 
| a 
The other necessary relations for the matched transformer and for the 
mismatched transformer with the same transmission characteristics, 
derived from equations (3), (8) and (26) are given below. % 
(a) Matched Transformer. 4 
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Mismatched Transformer: 
1 
Wh; 

1 Wo 
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re Fig. 5.—Mismatched transformer (a) fed by constant current generator, (b) with input 
and output reversed 
or | The ratio of the output voltages in the two cases is given by: 
=n / 2 2 ,/ R . 19 
Es y? FE 4 4 1 ( 2 E, 3] 
. ale one = ese - (5 
E, E, - 4R, R,¢ 9 E, 
The relative gains of the two transformers, as given by equation (31 
will be examined for three important cases: 
= (1) Input Transformer—In this case Ri = R,’, and C2 C,", so that 
© from (31), the mismatched transformer gives four times, or 6 db more 
gain. 





= (2) Output Transformer—The voltage source in Fig. 4 may be replaced 
@ by the constant current source /,' = E,’'/R,’, by Norton’s Theorem, as 





shown in Fig. 5(a). By the Reciprocity Theorem we may exchange the 


positions of the current source and the output voltage as shown in Fig. 5(b), 





¢ without changing the value of the. latter. This may also be done for the 
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matched transformer. The resultant circuit is that for a pentode « 
transformer in each case, where the plate resistance of the tube is s: 
compared to circuit impedance that we can consider it to be a co! 
current generator of g,,/, amperes, where g,, is the grid-plate transcor 
ance and £, the grid signal voltage. Thus in this case the mis-mat 
transformer also gives 6 db more gain than the matched transformer. 
(3) Interstage Transformer—If the voltage generators of Figs. 1 a 


are replaced by their equivalent current generators g,£, = E./R 
o a8 it , a“ 
QmE, = E,/R,, then (31) becomes, 


E:? _ , RiC2| gE, 
x RC. | gn Ey 


The first relation in (26) may be substituted in (32) giving, 
/ }2 v7! pl 8 

Es ey C1C 22m E, 

Ey 


 16CsQu Ey 
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t- Stati: 


ed 
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to 


number | 
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paper m 
“Mather 
properti 
33 The ct 


Thus, in the interstage case, this mismatched transformer has only 3 


db more gain than the matched transformer. 


Here, as in the case of the first mismatched transformer discussed in ¢] 
paper, there is an increased sensitivity to capacity variations. In practi 
it appears that this is not always serious in wide-band intermediate-f 


quency amplifiers. Some improvement in noise figure and power hand 
may be obtained, as before, without encountering difficulties with doul 
peaked pass-band characteristics. 

The design formulas given in this paper have been used successfully in 
design of experimental intermediate frequency amplifiers in the 65 mc reg 
having band widths of the order of 10 to 20 mc. 
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Statistical Properties of a Sine Wave Plus Random Noise 
By S. O. RICE 
INTRODUCTION 


N SOME technical problems we are concerned with a current which 
consists of a sinusoidal component plus a random noise component. A 
number of statistical properties of such a current are given here. The present 
paper may be regarded as an extension of Section 3.10 of an earlier paper, 
‘Mathematical Analysis of Random Noise’, where some of the simpler 
properties of a sine wave plus random noise are discussed. 
The current in which we are interested may be written as 


I = Qcos gt + Iy 
(3.4) 
= Rcos (qt + 9) 


where ( and gq are constants, ¢ is time, and /y is a random (in the sense of 
Section 2.8 of Reference A) noise current. When the second expression in- 
volving the envelope R and the phase angle 6 is used, the power spectrum of 
Ty is assumed to be confined to a relatively narrow band in the neighborhood 
of the sine wave frequency f, = q/ (27). This makes R and 6 relatively 
slowly (usually) varying functions of time. 

In Section 1, the probability density and cumulative distribution of / are 
discussed. In Section 2, the upward ‘‘crossings” of J (i.e., the expected 
number of times, per second, / increases through a given value /;), are 
examined. 

The probability density and the cumulative distribution of R are given in 
Section 3.10 of Reference A. The crossings of R are examined in Section 4 
of the present paper. 

The statistical properties of 6’, the time derivative of the phase angle @, 
are of interest because the instantaneous frequency of J may be defined to 
be fg + @/(2r). The probability density of @ is investigated in Section 5 
and the crossings of 6’ in Section 6. 6’ is a function of time which behaves 
somewhat like a noise current and may accordingly be considered to consist 
of an infinite number of sinusoidal components. The problem of determin- 
ing the ‘power spectrum” W(f) of 6, i.e., the distribution of the mean 
square value of the components as a function of frequency, is attacked in 

1 B.S.T.J. 23 (1944), 282-332 and 24 (1945), 46-156. This paper will be’ called 
‘Reference A’’, 
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Sections 7 and 8. The correlation function of 6’ is expressed in ter »s 
exponential integrals in Section 7, the power spectrum of /y being as 
symmetrical and centered on f,. In Section 8, values of W(f) are ob 
for the special case in which the power spectrum of /y is centered on / 

of the normal-law type. 

It is believed that some of the material presented here may find a 
the study of the effect of noise in frequency modulation systems. | 
example, the curves in Section 8 yield information regarding the noise | 
spectrum in the output of a primitive type of system. Also, the proc: 
employed to obtain the expression (5.7) for 6’ may be used to show t 


Qeos|(A, wo) Cos wot + gt] + Ivy = Reos (qt + 4) 
the sinusoidal component of d@/dt is* 
—A (1 — e *) sin wo 


where p is the ratio (”, (2 7). This illustrates the ‘crowding effect” of t/ 
noise. The statistical analysis associated with R and 6 of equations (3.4 
(when the sine wave is absent) is similar to that used in the examinatio 
the current returned to the sending end of a transmission line by reflect ions 
from many small irregularities distributed along the line. This suggests 
another application of the results. 
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1. PROBABILITY DISTRIBUTION OF A SINE WAVE PLUS RANDOM Nos! 
A current consisting of a sine wave plus random noise may be represented 

as 
I = Qeos gt + Ty (1.1 


where () and g are constants, ¢ is the time, and /y is a random noise current 
The frequency, in cycles per second, of the sine wave is fg = g/(2m). In all! 


2 The first person to obtain this result was, I believe, W. R. Young who gave it i 
unpublished memorandum written early in 1945. He took the output of a frequency 
modulation limiter and discriminator to be proportional to either the signal frequency 
to the instantaneous frequency (assumed to be distributed uniformly over the input band 
of /y according to whether () is greater or less than the envelope of /y . His memorandu 
also contains results which agree well with severa| obtained in this paper. 
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- work we denote the power spectrum of /y by w(/) and its Correlation 


| 


/ 


nection by ¥(r). The mean square value of /y is denoted by y 
rhe study of the probability distribution of J is essentially a study of the 


tegral® 


1 pe. Fy O cos 6] 
pl) | al | dé (1.2) 
TV Wo ~0 Vv 


where 


o(x) é ws ey 

V 29 
and p(/) is the probability density of /, ie. p(/)d/ is the probability that a 
value of current selected at random will lie in the interval /, / + d/. An- 
other expression for p(J) is given by equation (3.10-6) of Reference A, 


namely 


pil) | e ver” Jo(Qz) dz 1.4) 


26 J—« 


where Jo(Oz) denotes the Bessel function of order zero. 


The substitutions 


] QO 
) a 1) 
V Yo V Wo 
enable us to write (1.2) as 
1 r 
Ply) = WV Yo PU) = | gly — a cos 6) dé, (1.6) 
T 0 


where p(y) denotes the probability density of y. This is the expression 
actually studied. Curves showing ~;(y) and the cumulative distribution 
function 


I oy 
pU,)dI, = | pilvi)dy, 
v-- 2 “ x ° 3 ~ 
1.7) 
1 aT 
= (ta) iV — adcos @) dé, 
T 0 
where 
ax 
a) (0) = o(x,)dx, = , + 5 erf (x vV/ 2) (1.8) 
J— 2 


7W.R. Bennett, “Response of a Linear Rectifier to Signal and Noise,” Jour. Ay 


Soc. Amer. Vol. 15 (1944), 164-172, and B.S.7.J. Vol. 23 (1944), 97-113 
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are shown in Figs. 1 and 2. The curves fora = 10 anda = +/10 were TI 
puted by Simpson’s rule from (1.6) and (1.7), and the curves for a = | ® :y th 
were computed from the series (1.10) given below. Since both ¢(x) and 3 Se 
q orc 
LS } f i if T = x : 
| | | | | | | (in | 
a | fk | | pe 
1.0 | | | d 
o7s)}—_+—__1— —+———+ } 1 ——-f-----+- | ; 
r l a=00 
i 
0.5 ——J t+———— 1 cele eae 1 
0.4 fret ; =, 9¢ 
>_> | 
ee te ie apie 1 
| Sy = | y 9s 
] > } eee 
Pes es MO A a. +> 
| | | x 
01 Pt we x ] a6 
— | oO 
>» | J 
~— < 
a ou "7 ° 
—|lu | | 4 = 
41, 0.075 }—7+-—F | % 
© | | | - 
+ 7 
wy 0.05 -— T +—— es 
| 5 
_ “— [| a 
| < 
0.03 + ~ 4 —— —— + 3 
| j "4 
| a 
| x 
0.02 }— —+——+—--—- a? . 
| | - 
0.015 —+— T SY 
| 
0.0! | 
0.0075 | 
0.005 
6 062 CA O08 08 10 12 t@ 6 18 20 22 2A 28 
I1/ (RmSI) ‘ we 
Fig. 1—Probability density of sine wave plus noise. : 
4 
/ / — ei 
I = Qcos gt+ In, a = QO/V fo, vy = 1/V Yo, V vo = rms Ty : 
| bi(y) dI/V Yo = probability total current lies between 7 and J + d/ Ae 
| v(1 + a?/2)-'? = J/(rms J). Curves are symmetrical about y = 0. 5 
} 3 
¢—i(x%) are tabulated‘ functions the integrals in (1.6) and (1.7) are well i 
suited to numerical evaluation. , Tl 
a : re ‘ = in 
i 4 y(x) is given directly and g-,(x) may be readily obtained from W.P.A., “Tables 
Probability Functions,”’ Vol. II, New York (1942). The functions g‘”?(y) are tabulated 
in Table V of “Probability and its Engineering Uses” by T. C. Fry (D. Van Nostrand Co., 
1928). 
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The form assumed by /;(y) as the parameter a becomes large is examined 
n the latter portion (from equation (1.12) onwards) of the section. 

Series which converge for all values of a but which are especially suited 
for calculation when a < 1 may be obtained by inserting the Taylor’s series 
(in powers of x) for g(y + x) and g_i(y + x), x = —a cos @, in (1.6) and 
(1.7) and integrating termwise. When we introduce the notation‘ 














m m 
(m d 1 d y2/2 ‘ 
¢ (y) = —¢y) = > ps2 (1.9) 
dy™ V 2x dy™ 
99.99 
* 99.98 -—-—_ +>} ot 
V 
~ 99.9-—+—T a aia otc! 
= 
w 99.8-—+——} T ; a ah 
z } 
s 
| en ee + t t 
3 99;— —+---— +-- co - T - + 
e | 
© oon | } 4 4 4 
< 
- 95+-——+ +— 4 —__} + 
re } 
= eet see a —— } 4 ah } } } 
3 90 
a 
< | s 
¢ 80K + + —+ + + + + + + + 
ac | - ] 
a | "said | | 
as 70 -————_—+—- _. << — a + + + +- + se S| 
2° i 
° 60 4+-——_——_--+ + + + + T - + + + 
50 1 











) 02 04 06 O8 10 12 +64 #+%1'16 %18 20 22 24 26 
I1/(RMS TOTAL CURRENT ) 


~ 


Fig. 2—Cumulative distribution of sine wave plus noise 


y 
Ordinate = 100 pi(yi) dy, . See Fig. 1 for notation 


x 


we obtain 


x 1 2n ; 
nao min! \2 (1.10) 


J 


ie = 1 : een 
pilyvdvy = galy) + ps <a gy) 
x n=i Min! \Z 


The second equation of (1.10) may be shown to be valid by breaking the 
interval (— «, y) into (— «, 0) and (0, y). In the first part, 





| pi(vi) dy; = o-1(0) 
ce) 
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since both sides have the value 1/2. In the second, term by term int 
tion is valid since the series integrated are uniformly convergent as m 


1/2 ) 1/4 D 
) - (1 + O(n 1) + Oly 2 I], 
£m TH 


in which we suppose that y remains finite as—> x. This may be obtai 


seen from the inequality 


by using the known behavior of Hermite polynomials of large order.° 
When 0 > rms I, so that a is very large the distribution approaches t 
of a sine wave, namely 


ply) ~ * 


9 


a= ¥y) 


1/2, 


r 1 1 ea 
pi(m)dya ~ 5 + —arcsin<, y| <a 
a) “ T a 
In order to study the manner in which the limiting expressions (1.12) 


approached it is convenient to make the change of variable 
x= y— @ cos ¢, d@ = {a? Si x)?] Ll2 dx 
Z2=2x-yta 


in (1.6). We obtain 


ny ta 
pily) = ; | g(x) [a — (y — x)? dx 
T ¥y-a 
| ¢g(s + y — a)[s(2a — z)] 
T 0 

An asymptotic (as a becomes large) expression for ~;(y) suitable for the 
middle portion of the distribution where a — | y > 1 may be obtained from 
the first integral in (1.13). Since the principal contributions to the value of 
the integral come from the region around x = O we are led to expand the 
radical in powers of x and integrate termwise. Legendre polynomials enter 
naturally since they are sometimes defined as the coefficients in such an 
expansion. Replacing the limits of integration y + a and y — a by +~ 
and — «, respectively and integrating termwise gives 


2 2 1/2 2% : 
(a — vy) « 1.5.9. (2 =f) 
pity) ~ : , E + > (-) r > I 
‘ 7 nel er == oy)” 


2(a* — ¥°*) ae — 7} 


a 3 ewe 3t 3(350 + 301 + 3 
(a 1+ gett, + Pet 


5A suitable asymptotic formula is given in Orthogonal Polynomials, by G. Szegi, 
Am. Math. Soc. Colloquium, Pub., Vol. 23, (1939), p. 195. 
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ere? = y*/(a? — y’) and P2,( ) denotes the Legendre polynomial of 
er 2n. We have written this as an asymptotic expansion because it 


viously is one when y, and hence /, is zero in which case 
1.3.5---(2n — 1) 


Po, (0) = (-—)’ 
2.4.---2n 
When vy is near a or is greater than a, a suitable asymptotic expansion may 


be obtained from the second integral in (1.13) by expanding (2a 


in powers of z/ (2a) and integrating termwise. The upper limit of integra- 


tion, 2a, may be replaced by ~ since g(s + vy — a) may be assumed to be 


negligibly small when z exceeds 2a. We thus obtain 


where we have used the notation 
(a@)o 1, (a) ala+i1)--- ( nl 


The integrals occurring in (1.15) are related to the parabolic cylinder 
function® D,,(x). Their properties may be obtained from the known 
properties of these functions or may be obtained by working directly with 
the integrals. 

Suppose now that a is very large so that only the leading term in the series 
(1.15) for pi(y) need be retained. 

Then 


where 


'? dz (1.17) 


By writing out ¢(s + s), expanding exp (—zs) in a power series, and inte- 
grating termwise we see that 


2 
g(s)2 ae ( 


T (=0 


Fis) = 


' me 
= (2r) 's' “ols V 2) Ky(s°/4) 
where AK denotes a modified Bessel function. The relation (1.18) may also 
6 Whittaker and Watson, ‘“‘Modern Analysis,” 4th ed. (1927), 347-351 


7A table of A4fx) is given by H. Carsten and N. McKerrow, Phil. Mag. $7, Vol 
(1944), 812-818 
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be obtained from pair 923.1 of ‘Fourier Integrals for Practical A 
tions,” by G. A. Campbell and R. M. Foster.® 

A curve showing F(y — a) plotted as a function of y — @is given in | 
It was obtained from the relation 


| 


F(s) = 2 4a-32x(— s/4/2) 




















where 
@ 
z—w?)2 
x(x) = [ e "au 
“0 
0.24 
0.22 ee a 
0.20} 4 } 4 } 
0.18 + + . + + 
| 
0.16 4 a 
} 
ax 0.14 + + t 
© | | 
' } 
> 0.12 4 4 t 
~ | | 
uw | 
0.10 + } + + + = 4 4 4 + 4 
| | 
| } 
0.08 + + ' + + 
| | 
0.06 }+——+ + ee | + + + + t + 
0.04 F | 4 
0.02 +—+ ' + - +-\b-—- 
r¢) 1 J | 
-4 -3 2 -| 0 \ 2 3 
Y-a 
Fig. 3—Probability density of sine wave plus noise. 
When rms Iy << Qand / is near Q, p(y) ~ a72F (y a),y—-a=(I Q)/(rms Ixy 


See Fig. 1 for notation. 


This function has been tabulated by Hartree and Johnston.’ 
The probability that J exceeds (), or that y exceeds a, is, integrating the 
second of expressions (1.13), 


oO 4% dz + 
Oy» , steiae 
i ily) 4 rd V/2(2a — 2) J, * na 


An asymptotic expansion may be obtained by expanding (2a — z)~!/? as in 
the derivation of (1.15) but we shall be content with the leading term 


* Bell Telephone System Monograph B-58A. 
* Manchester Lit. and Phil. Soc. Memoirs, v. 83, 183-188, Aug., 1939. 
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sing 
oS 
wo a a® pr 
/2 1/2 1/4 —1/2 yy 
i has dz | g(x) dz = | g(x) dx | 2 dz 24 r(3) 
0 Jz “0 “0 
we obtain 
a 
\ 1/4 3/2 ww 2 1/2 = 1/2 
| ply) dy ~ 2 te PP Q)at? = 0.185- +a (1.19) 


For use in computations we list the following values 
r(4) = 3.62561, r(?) = 1.22542, r(2) = 0.90640 
2. EXPECTED NUMBER OF CROSSINGS OF I PER SECOND 


In this section, we shall study two questions. First, what is the proba- 
bility PU, t:)dt of J increasing through the value /, (i.e. of J passing 
through the value /; with positive slope) during the infinitesimal interval 
t1,t; + dt? Second, what is the expected number .V(/,) of times per second 
] increases through the value 7;. When /; is zero, 2N(0) is the expected 
number of zeros per second, and when /, is large \V(/,) is approximately 
equal to the expected number of maxima lying above the value /, in an 
interval one second long. 

We start on the first question by considering the random function 


z= F(a,, do, +++ ay: t) 


where the a’s are random variables. The probability that the random curve 
obtained by plotting z as a function of é increases through the value s = 2; 
° . - 1 

in the interval /;, ¢; + dtis 


ax 


at | np(zi, 7; 4) dn (2.1) 
0 


where p(é, 7; ¢:) denotes the probability density of the random variables 


& = F(a, de, +++ , dy; ty) 
OF 
»=I/- 
Ot jt=t,. 
In our case z becomes the current J detined by equation (1.1). The 


method used to obtain equation (3.3-9) of Reference A may also be used to 
show that the quantity p(/,, , ¢:) (which now appears in (2.1)) is given by 


WwNo e 
PU, n, 4) = mely — acos gt)g(x + 6 sin qt;) (2.2) 


0 


10 This result is a straightforward generalization of expression (3.3-5) in Section 3.3 of 
Reference A where references to related results by M. Kac are given. A formula equiva 
lent to (2.1) has also been given by Mr. H. Bondi in an unpublished paper written in 1944 
He applies his formula to the problem studied in Section 4. 
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where ¢( ) denotes the normal law function defined by equation (1.4 


a 


” 2 ~ 22 . . 1 —ypo I; 
—-W = 4° | wif )f df, N o> ’ y= oe ; 
¥ Jo te ; T { Yo j Vw 
' dv . n Qq 2af, : 
( = e 4 = 9 = 5 i >= = . 
.” Yo Vv —ypo Vv —ypo N 0 


? 


Equation (3.3-11) of Reference A shows that .Vo is the expected numb 
zeros per second which /y would have if it were to flow alone. 

Let P(/,, t;)dt be the probability that J will increase through the value / 
during the interval ¢;, 4; + dt. Then (2.1) and (2.2) give 


PU, t,) = | np, n, ti) dn 
“0 
is 24 


= tNog(yv — a cos gt) i xg(x + b sin qt) dx. 
0 


The integral in (2.4) is of the form 


sa 2 
| xgo(x + v) dx = ¢g(v) — 2 [ g(x) dx 
0 


vv 


-" + o(v) +9 [ g(x) dx 


“0 


= —v+ o(v) + v9_;(v) 


a ae tal 14 
= —_ + (r) ($553 -*) 


where v replaces 6 sin gf; and ,F, denotes a confluent hypergeometric func- 


tion. 


The distribution of the crossings at various portions of the cycle (of the 


sine wave) may be obtained by giving special values to gf, in (2.4). 


The expected number of times J increases through the value /,; ia one 


second is 


ao © i a ee 
N(I) = Limit | Puls. t;) dl, 
t+0 1 Jo 


(2.6 


Ls 


bsin9 
= No] ¢ly — acos 8) E sin 0) + dbsin@ [ g(x) ax |ao 
Jo 


0 

where we have used (2.4) and the second equation of (2.5). The integrand 
in (2.6) is composed of tabulated functions and is of a form suited to nu 
merical integration. Expanding ¢(y — a cos 6) in (2.6) as in the derivation 
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1.10), replacing the quantity within the brackets by the series shown in 
e last equation of (2.5) ,and integrating termwise leads to 


V(Ih) = No(w/2)'? yo? (5) Fi (-3: +4; -5) (2 


n=o nin! 2 


he series (2.7) converges for all values of a, y,and 6. This follows from the 
inequality (1.11) which may be applied to ¢ ’(y), and from the fact that 
the ,/; is less than exp (6/2) as may be seen by comparing corresponding 
terms in their expansions. 

The expected number of zeros, per second, of / is 2.V(O) where we have set 
],, and hence y, equal to zero. In this case the integral in (2.6) may be 
simplified somewhat and we obtain 


2N(0O) = No « “To(B) + : le (°. a) (2.8) 
4Q@ a 


where / (8) is the Eessel function of order zero and imaginary argument and 


a+b a —b 
* : 4 


Te(k, x) = e “To(ku) du. 
Jo 
The integral /e(k, x) is tabulated in Appendix I. 

I have been unable to obtain a simple derivation of (2.8). It was orig- 
inally obtained from the following integra] 


3 2 


, N, 
NU) = - | d0 o(y — a cos @) | xg(x + bsin 6) dx = (2.9) 


“0 


which may be derived from the second equation of (2.4) and the first of 
(2.6). Setting /, and y equal to zero and writing out the y’s gives 


p@® 


2N(0) = = | do | dx 
Tas —* 


~~ 
x exp [—3(x° + 2bx sin 6 + a’ cos 6 + B sin’ 6)]. 
Equation (2.8) was obtained by applying the method of Appendix III to 
this expression. 
3. DEFINITIONS AND SIMPLE PROPERTIES OF R AND O 


The remaining portion of this paper is concerned with the envelope R and 
the corresponding phase angle 6. These quantities are introduced and some 
of their simpler properties discussed in this section. 

Suppose that the frequency band associated with Jy is relatively narrow 
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and contains the sine wave frequency f,. The noise current may 
solved into two components, one “‘in phase” and the other “in quadra: ur 
with Q cos gt. Using the representation (2.8-6) of reference A and provee 
ing as in Section 3.7 of that paper: 

M 


Iy = > Cn COS (Wnt — gn) 3] 


n=l 


M 
= 2. Cn COS [(wn — g)t — gn + ql] 
4 == 1 
= I, cos gt — I, sin gt 3.2 
where 


M 
i, = =. Cn COS [(Wn — g)t — ¢nl 


n=1 


M 
I, = =: ce sin [(w, — g)t — ¢nl 


n=l 
a, = arf ny Te _ nf, és ans 2w(fnJAf 


w(f) denotes the power spectrum of Jy and the ¢,’s are random variables 
distributed uniformly over the interval (0, 27). 
The total current J may be written as 
I = Qcosqt+ Iv 
= (0 + I.) cos gt — I, sin gt 
= R cos 6 cos gt — R sin @ sin gt 
= Ros (gt + 8) 
where we have introduced the envelope function R and the phase angle @ 
by means of 
Reos#=QO+ IT, 
Rsin@d =I, 

Since J, and J, are functions of ¢ whose variations are relatively slow in 
comparison with those of cos gt, the same is true of R and (usually) @. 

A graphical illustration of equations (3.4) and (3.5) which is often used is 
shown in Fig. 4. 

In accordance with the usual convention used in alternating current 
theory, the vector OQ is supposed to be rotating about the origin O with 
angular velocity g. If Jy happened to have the frequency g/2z, its vector 
representation OT would be fixed relative to OQ. In general, however, the 
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enuth and inclination of OT will change due to the random fluctuations of 
.. Thus the point 7 will wander around on the plane of the figure. If 


rms Jy is much less than (), 7 will be close to the point ( most of the time 


In this case 


6 = tan” oe 
t a \ 
do di, _ iI, 
dt dtQ Vv 
Le 


and a number of statistical properties of R and @ may be obtained from th‘ 
' and 'f 


corresponding properties of noise alone when we note that /,, / 
behave like noise currents whose power spectra are concentrated in the 


lower portion of the frequency spectrum. 





OS=Q+tIc=RCOS9 T 
QT=In 





e) —» 


Fig. 4—Graphical representation of J = Qcos gt + In. 

By squaring both sides of equations (3.1) and (3.3) and then averaging 
with respect to ¢ and the ¢,’s we may show that J,, 7,, and /y all have the 
same rms value, namely yo”. 

It may be seen from (3.3) that the power spectra of 7. and J, are both 
given by 

7) 


wfa tf) + wlfe—f) 3. 


where it is assumed that 0 < f <f,. Likewise the power spectrum of the 
time derivative J; of I, is 
4’ f*[w( fg t+ f) + whe — fl (3.8) 

This follows from the representation of /{ obtained by differentiating the 
expression (3.3) for 7, with respect to ¢, the procedure being the same as in 
the derivation of equation (7.2) in Section 7. The power spectra shown in 
Table 1 were computed from equations (3.7) and (3.8). 

The correlation function for 7, and hence also for /,, is, from equations 
(A2-1) and (A2-3) of Appendix I, 
= 
| w(f) cos 2r(f - ‘te \r df 
“U 


I(t +7) = g = 
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where the bar denotes an average with respect to ¢ and g is a functic 
From (A2-3) the correlation function for // is —g” where the doubk 
on g denotes the second derivative with respect to 7. 

Attention is sometimes fixed upon the variation in distance betwee su 
cessive zerosof/. The time between two successive zeros of J at, say, fy .nd 
is the time taken for gf + 6, as appearing in R cos (qt + 4), to increase )y + 
This assumes that the envelope R does not vanish in the interval. Ff 
moment we write @ as @(¢) in order to indicate its dependence on the time 
Then ¢ and ¢; must satisfy the relation 


qty + O(¢1) ‘cna gto _ A(to) = 49 
Since 6(1) is a relatively slowly varying function we write 


O(t;) — O(to) = (ty — to) O’(to) + (ty — to)?O” (to) /2 + 


TABLE 1 
POWER SPECTRA OF ly ,J]¢-, 1/5, AND 1; 


In Te and T/s I; 
w(f) = wo = Yo/B for f, — B/2 <f | 2wofor0 <f < B/2 \8x*f?w forO < f < B/2 
< fa t+B/2 
w(f) = 0 elsewhere 0 elsewhere 0 elsewhere 
fq = mid-band frequency 
w(f) = wo = Yo/Bforf, —B<f<f,| wofor0<f <p 4 fw for0<f <8 
w(f) = 0 elsewhere 0 elsewhere 0 elsewhere 
f, = top frequency 
Yo f—f,)2/ (2e2 2Yo 3 f2/ (292 Sr2f2Wo _ 52/292 
w(f) = e ’ 7—e* 7 6 
oV 2n oV ln oV 2x 


where the primes denote differentiation with respect to ¢. When this is 


placed in (3.9) and terms of order (¢; — ¢o)” neglected, we obtain 


1 ri 1 
= 1 + — @'(t) (3.10) 
Q(t; ei to) dr Qn 
which relates the interval between successive zeros to 0’. 
The expression on the right hand side of (3.10) may be defined as the in- 
stantaneous frequency: 


1 dé 


. (3.11) 
2m dt 


Instantaneous frequency = f, + 


This definition is suggested when cos 2xft is compared with cos (qt + 4 
and also by (3.10) when we note that the period of the instantaneous fre- 
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uency 1s approximately equal to twice the distance between two suc¢ essive 
eros which is 2(¢; — fo). 
lhe probability density of R is" 
R R+0 
exp i — = | Io RO/ Wo) (3.12) 
v0 - 0 
v ‘ by 
where Jp (RO/Wo) denotes the Bessel function of order zero with imaginary 
argument. In Section 3.10 of Reference A, it is shown that the average 


value of R" is* 


R = (2y)"'r r+ 1) (—%, 1; -). (3.13) 


where p = (7, (2Wo), of which special cases are 


R=e “(ao 2)" 11 + p)lo(p/2) + pli(p/ 2)| 


R? _ ( )? + 2Wo 


(3.14) 
Curves showing the distribution of R are also given there. 


4. EXPECTED NUMBER OF CROSSINGS OF R PER SECOND 


Here we shall obtain expressions for the expected number .Vz of times, 
per second, the envelope passes through the value R with positive slope. 
When R is large, .Vx is approximately equal to the expected number of 
maxima of the envelope per second exceeding R and when R is small .Vx is 
approximately equal to the expected number of minima less than R. For 
the special case in which the noise band is symmetrical and is centered on 
the sine wave frequency f, .Vx is given by the relatively simple expression 
(4.8). 

The probability that the envelope passes through the value R during the 
interval ¢, ¢ + dt with positive slope is, from (2.1), 

p@ 


dt | R'p(R, R’, 8) aR’ (4.1) 
aT) 

where p(R, R’, t) denotes the probability density of R and its time derivative 
k’, t being regarded as a parameter. 

An expression for p(R, R’, t) may be obtained from the probability density 
of T.,1,,1,, 1. From our representation of a noise current and the central 
limit theorem it may be shown (as is done for similar cases in Part ILI of 
Reference A) that the probability distribution of these four variables is 


4 In equation (60-A) of an unpublished appendix to his paper appearing in the B.S.7.J. 
Vol. 12 (1933), 35-75, Ray S. Hoyt gives an integral, obtained by integrating (3.12) with 
respect to R, for the cumulative distribution of R. 

*The correlation function for the envelope of a signal plus noise, together with associated 
probability densities of the envelope and phase, is given by D. Middleton in a paper 
appearing soon in the Quart. J]. of Appl. Math. 
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normal in four dimensions. If the variables be taken in the order «jy, 


above the moment matrix is, from equations (A2--2) of Appendix II, 
bo 860 0 db | 
0 by —b, O + 


M = 


0 —b b, 0 








b, 0) 0 b 


where the 6’s are detined by the integrals in equations (A2-1). The i: 
matrix is 
if & 0 0 —b] 


M" = , B= bobo — bi $3 
B 0 bb b O 


0 ob bo 0 








g —b, 0 0 bo | 


which may be readily verified by matrix multiplication, and the determinant 
M | is B*. The normal distribution may be written down at once wher 


use is made of the formulas given in Section 2.9 of Reference A. The sub 



















stitutions 
I. = Reos 6 — Q, I, = R’cos 6 — Rsin 0 0’ 
I, = Rsin 8, I, = R’sin 0 + Reos 6 6’ ™ 
dI di dI-dI, = R*dRdR'd6d6’ 
enable us to write 
bol? + 12) + bo? + 12) 
—2b,(1 J, — 1412) = bo(R? — 20Rcos 6 + 02) 
+bo(R”? + R60’) 
— 26, R°0’ + 26,0(R’sin 6 + RO cos 6). 
Consequently the probability density of R, R’, 0, @’ is 
p(R, R’, 0, 6’) = nd exp { ie [bo(R° — 20Rcos 6 + Q°) 
4n°B | =a or 


\ 


+ b(R” + R’0") — 26, R60’ + 26,Q0(R’sin 6 + RO'cos 8) Y 


In this expression R ranges from 0 to x, 6 from —z to m, and R’ and @ 
from —* to +2, The probability density for R and R’ is obtained by 






ee 
‘Sirichi Ss eecasitacak see tlk Fomeee 


pe een ne 


th 
in 











1 ‘ : — 
4 exp 4 — [B(R’ — 2RQcos 6 + Q°) + (b)R’ + b,Qsin 6)*]> 
\ 


; bo\"? R RO | R+ Q 
Nr = (*) b Io ( bo ) exp ( 2b, -) 





VE 
4) 
$3 
inant 
whe! 
sub 
4.4 
k 

15 
d 6 
by 
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ntegrating (4.5) with respect to @ and @& over their respective ranges. The 
integration with respect to @’ may be performed at once giving 


} _-«-RQey** ¢* 

® (RR, t) = | d0 

Ej V Bh ~* : 
4 (4.6) 
m expi- [B(R? — 2ROQcos @ + (0) + (bo R’ + b,Qsin 6)°| $ 

4 j 2Bbo ‘ : 


as 


a From (4.1) and (4.6) it follows that the expected number Vz of times per 
Fl. ‘ aoa . 
® second the envelope passes through R with positive slope is 
_ R(2x)** [* te ; 
Ne = d0 R’dR 


Bbo i “0 
i (4.7) 


\ 


2Bby 


When the power spectrum w(/) of the noise current Jy issymmetrical about 
the sine wave frequency f., 6; is zero and B is equal to bobs. In this case the 
integrations in (4.7) may be performed. We obtain 


7 (4.8) 
1 a Probability density of 
= 4 ‘3 
Qn envelope at the value R 
where the second line is obtained from expression (3.12). As will be seen 
from its definition (A2-1), 69 is equal to the mean square value Wo of Jy 


(and also of J, and J,). 
Introducing the notation 


1/2 
v = Rbo ~ = R/rms Iy 
st (4.9) 
a= Ab ~ = QO/rms Iy, 
which is the same as that of equations (3.10-15) of Reference A except that 
there P denotes the amplitude of the sine wave and plays the same role as 
( does here, enables us to write (4.8) as 


; be a 2442)/2 be ie 
Nr = | | v]y(avye = | | Pie). (4.10) 


The function p(v) is plotted as a function of v for various values of a in Fig. 
6, Section 3.10, of Reference A. 
It is interesting to note that 


(b2/bo)!/*/a = Expected number of zeros per second of J, (or of J,) (4.11) 
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This relation, which is true even if the noise band is not symmetrical a 
f,, follows from equation (3.311) of Reference A. 

When ( > rms /y and f, is not at the center of the noise band it is « 
to obtain the asymptotic form of Vx from the approximation (3.6), 


R~O+ I, 


Part of 1 
pe yrtion 


instead of the double integral (4.7). When O>> rms /y and R is i: 
neighborhood of ( (as it is most of the time in this case), .V x is approxima 


a second | 
equal to the expected number of times /,, increases through the value / 


R — (in one second. Thus, regarding /, as a random noise current 

have from expression (3.3-14) of Reference A 

12, / (2b 
1 


Ve~we X [1,2 the expected number of zeros per second oj 


and when we use equation (4.11) we obtain 
1 1/2 
N Rr ~~ (bo bo) “€ 


2r 


TABLE 2 


) = by /B OVER A BAND OF WipTH B 


The | 
Band extends from f, 8/2 to (3 /6)!/2 = ().7248; equatio 
B/2 


/@q 


Same as 1 and in addition Q = 0 
Same as 1 and in addition 
QO>>rms/, 


As w 
Band extends from f, to f, + 8 


@ assoc 
and Q >> rms/, 


quency 
modifie 


- ' . given v 
lable 2 lists the forms assumed by (4.10) and (4.12) when the power spe: 


: band o' 
trum w(/) of the noise current /y is constant over a frequency band of width 3. 


Ke ; ; : 7 : those fe 
Phe quantity 69 in the expressions for b, represents the mean square value 


of Ty. 

In the general case where the band of noise is not centered on /, and 
where R is not large enough to make (4.12) valid we are obliged to return to 
the double integral (4.7). Some simplification is possible, but not as mu 


relative 
The 
sion (4. 
respect 
for R, ¢ 
as could be desired. Introducing the notation 
= RO/bo, y = b,0(Bby)>”? PLR, 8, 


(boR’ + 6,Osin 0)(Bbo) 
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ibles us to write (4.7) as 
3/2 . 12 R2+Q2 
(27) © “(R/bo)(B/bo) “e 


7 


dé (x — y sin Oe"? =? dy 


v--4¥ ~“ y sin 9 


Part of the integrand may be integrated with respect to x and the remaining 
portion integrated by parts with respect to 6. The double integral in the 
second line of (4.13) then becomes 


b/2 ? l 
dx ajya e 


1 x cos 6 6) 2/2 
cos Aye‘ dé 


7 


= ya e ia ‘7 (y cos 6 + a/y)e 


r 


= (3) y . 
= 27 b — (- : [T,(@) + 
n= M1: a 


The series is obtained by expanding exp |—(y sin @)?, 2] 


yY cos 8 


equation in powers of sin @ and integrating termwise. 


” dé 
5. PROBABILITY DENSITY Of di 
As was pointed out in Section 3 the time derivative @ of the phase angle 
@ associated with the envelope is closely related to the instantaneous fre- 
quency. The probability density p(6’) of 6’ may be expressed in terms of 
modified Bessel functions as shown by equation (5.4). Curves for p(6’) are 
given when the sine wave frequency /, lies at the middle of a symmetrical 
band of noise. Although the expressions for p(6’) are rather complicated, 
those for the averages 6’ and | @’| given by equations (5.7) and (5.16) are 
relatively simple. 
The probability density p(@’) may be obtained by integrating the expres- 
sion (4.5) for p(R, R’, 0, 6’) with respect to R, R’, 6. The integration with 


respect to R’, the limits being — x and + ~, gives the probability density 
for R, 6, 6: 


St fa \1/9 : 
P(R, 6, 6’) = z (, ) exp | —aR° + 2bRcos 6+ csin @ 
r: 


i) 


—b.0"/(2B)| 
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where 
B = bobe — bj b = O(b. — b,0’)/(2B) 
a = (bp — 26,6’ + bo8’)/(2B) — ¢ = 0°1/(2Bbo) = bip/B 


p= 0? (2bo) YY = b a : 
The int 
and b¢, 1, b: are given in Appendix IT. the inte 


Integrating with respect to R gives the probability density for 4, / detiniti 
Expanding exp (2bR cos @) in powers of R and integrating termwise, 


I 2 si c sin20—bobop/ B 
6, 6’ = je sin robop 
P\ ; 167a is) ; 


> n+1 ‘ cos ey 
ned of 08 oe 
I (; > ') We ma 
P appear: 
When we integrate @ from —7 to r to obtain p(6’) the terms for which » 
is odd disappear and we have to deal with the series, writing y for 6°/a, 
oO Whe 
2m a 1 9 9 9 
> ae (y cos 6)" = (2y cos 6 + 1) exp (y cos’ @) 
m=( H 


Thus, the probability density of @’ is 


1 2 we gr 9 c 8in20+¥ cos? 6—b»bop/ B 
ie 16ra \ab)B [ (2y cos 6 + 1) &' ae eee ae 
( 0 Jr 
In tl 


1 a gs y—-c y-c Ba 
a. (y+ 1)/ - (5. | noise ¢ 
8a (5) | . , ( 2 ) - ( 2 ) case by 


exp k 7 —_ bs > 
= B 


_ bo <— 2b,0' 
P bs — 26,0’ + B06” 


From (5. 


y= < 


c+y ea bobop —e be — 26,0’ + 2bp” where 
: B 2 b: — 2b,’ + 6/6” 


It will be noted that for large values of | 6’ | the probability density of 4’ 
varies as | &’ |-*. Although this makes the mean square value of 6’ infinite, 
the average values @’ and | 6’ | of 6 and | 6’ | still exist. In order to obtain 
@’ it is convenient to return to (4.5) and write 

Pa re 2 +o x } 
= J dé [ dR i] dR’ | dé’ 6’ p(R, R’, 6, 6’) (5.6 \ 2 T} 
-¥ “0 — © a) 


Bennet 
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e integration with respect to 6’ may be performed by setting R@’ equal to 
nd using 
” ae (3/a)(m/a 
in 


Phe integral in R’ reduces to a similar integral except that the factor 


the integrand is absent. Performing these two integrations and using the 


definition of B leads to 
aT = 
dé | dR(R — V cos @) 


aT) 


1 

exp x (R° — 20R cos 6 + O 

2b é . 
— 

When this is done cos @ 


We may integrate at once with respect to R. 
rhus 


appears and the integration with respect to 6 becomes easy. 


6’ = (b1/bo) exp [—Q?/ (2b0)| = (b1/ bo) 


When the noise power spectrum is equal to wo in the band extending f 


ar 


— B/2 to fo + 8/2 and is zero outside the band, b, 


Hence, from (3.11), 


ave. instantaneous frequency = f, + 6’/(2 


(27) 


= fot Ue — fo)(l — 


qd 


In the remainder of this section we assume the power spectrum of 
noise current to be symmetrical about the sine wave frequency / In 
case 6, and ¢ are zero, B is equal to bob2 and (5.4) becomes 


p(8") = 4(bo/be)'2(1 + 2?)-*%e 


[iy + 1)Jo(y/2) + ylily 
5 (bo/be)'/2(1 + 3?) 


where ,/'; denotes a confluent hypergeometric function” and 
2? = bo6'*/be, > = What = 2 5.10 


When the noise power spectrum is constant in the band extending from 


f, — B/2 tof, + B/2 (see Table 2, Section 4) 


(bo/bo)'/? = 3B, 


® The relation used above follows from equation (66 
Bennett’s paper cited in connection with equation (1.2 
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Fig. 6—Cumulative distribution of time derivative of phase angle. 


Notation explained in Fig. 5. 
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rhe probability density p(6’) of # and its cumulative distribution, ob 


ed by numberical integration, are shown in Figs. 5 and 6. 


The probability that @’ exceeds a eiven 4; is equal to the prob ibility 


exceeds 31, where 2; denotes (by) 6 #,,and both probabilities are equal to 
—— . 3 ie : 
1+s) iF; -1; p(1 + lz oa 
) Sx. ? 
[he probability that 4° > 6; becomes 21) as {1 — x 


When () > rms /y the leading term in the asymptotic expansion of the 
Fy in (5.9) gives 
bia’) ~ is ; ; h,/Q) 5.13 
i aN in y t 


when it is assumed that 3° < 1. This expression holds only for t 


portion of the curve for p(@’). Far out on the curve, (6) still varies as 
0%. Equation (5.13) may be obtained directly by using the approximation 
3.6) that 6’ is nearly equal to /; © and noticing that 6, is the mean square 
value of I. 

If the sine wave is absent, p is zero and 


b( 0’) lih, hb 1 7 


which is consistent with the results given between equations (3.4-10) and 


3.4-11) of Reference A. In this case (5.12) becomes 


~ 
~ 


\lthough the standard deviation of 4’ is infinite an ide: 
the distribution may be obtained from the average value of 6’ Seiting 
6; equal to zero in (4.5) in order to obtain the case in which the noise band is 


symmetrical about the sine-wave frequency leads to 


) —F _T . x 4 
4’ a | dR | dé | dR’ | le’ a’ R 
» J | a 4 - 


exp 3 |- R° - 20R cos 0 O°) /b R™ + RO) /bo) 


he integrals in R’, 6’ cause no difficulty and the iniegral in @ is proportiona 
to the Bessel function /)(OR bo). When the resulting integral in R is 


evaluated we obtain 


where p = (?/(2bo). 


18 See, for example, G. N. Watson, ‘'T ry of Bessel Fur ms,’ Cambr 1044 


394, equation (5). 
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When p is zero equation (5.16) agrees with a result given in Section 
reference A, namely, for an ideal band pass filter 


t = %% he a Ta 


T1 ; V 3(fo + fa) 





ave 


where 7 is the interval between two successive zeros and 7, is its ave 
value. 7 is equal to ¢; — fo of our equation (3.10) from which it follows + 


(7 —71)/n = — 0 g 17 


6. EXPECTED NUMBER OF CROSSINGS OF @ AND S PER SECOND 
( 

\fter a brief study of the expected number of times per second the plas 
angle @ increases through 0 and through mr (where it is assumed that —7z 
@ < m) expressions are obtained for the expected number .Vy- of times pe 
second the time derivative of @ increases through the value 6’. 

The point 7 shown in Fig. 4 of Section 3 wanders around, as time goes by 
in the plane of the figure. How many times may we expect it to cross some 
preassigned section of the line OQ in one second? To answer this problem 
we note that, from expression (2.1), the probability that @ increases throug! 
zero during the interval ¢, ¢ + dt with the envelope lying between R and 
R + dR is 


2 

dt aR | 6’ p(R, 0, 6”) do’ (6.1 
“0 

where the probability density in the integrand is obtained by setting @ equa 

to zeroin equation (5.1). The expected number of such crossings per second 


IS 


’ 


3/2 ee by (R--Q)2/ (2B) 
dr) (by B) rine” ‘ 


( 


« 


| do’ 6’ exp [—by R°0”/(2B) + 6, R(R — Q)0'/B] 
which may be evaluated in terms of error functions or the function ¢ 
detined by equation (1.8). For the special case in which the power spe 
trum of the noise current Jy is symmetrical about the sine wave frequency, 
b, is zero and (6.2) vields 
$/27—1,1/2 —(R—Q)2/ (2ho) , 
(29) "bo bo “e : dR (6.3 

From equation (6.1) onwards we have tacitly assumed that the range of 4 
is given by —a < 6 < m because setting 6 equal to any multiple of 27 in our 
equations leads to the same result as setting @ equal to zero. This is due to 
6 occurring only in cos @ and sin 6. When @ increases through the value 7, 
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it does when the point 7 crosses, in the downward direction, the extension 
of the line OQ lying to the left of the point O in Fig. 4, we imagine the value 
of 6 to change discontinuously to the value — 7. 

lhe expected number of times per second @ increases through may be 
obtained from (6.2) and, in the symmetrical case, from (6.3) by changing the 
sign of Q since this produces the same effect as changing @ from 0 to 7 in 
p(R, 0, 6’). 

The expected number of crossings per second when R lies between two 
assigned values may be obtained by integrating the above equations. For 
example, the number of times per second @ increases through zero with R 
between Q and R;, is, from (6.3) for the symmetrical case, 


(47) l(b, bo)! 2 erf [(2do) V2) Ry — () } (6.4) 


5 


where we have used the absolute value sign to indicate that R; may be either 
less than or greater than Q and 


1/2 t? ™ 
: erf x = 2r e dt 6.5) 


Expressions for b) and b, are given by equations (A2-1) of Appendix IT, 
The mean square value of Jy is bo, and when the power spectrum of /,y is 
constant over a band of width 8, by = 7°67b9/3. 

In much the same way it may be shown that the expected number of times 
per second @ increases through 7 with R between 0 and R;, is 


(4or)"(bo/bo)'/? ferf [(2b9)—!/2( Ri + O)| — erf [(2b5)70}| 6.6) 


A check on these equations may be obtained by noting that the expected 
number of zeros per second of J,, given by equation (4.11), is equal to twice 
the number of times 6 increases through zero plus twice the number of times 
9 increases through 7. Setting R, equal to zero in (6.4), to infinity in both 
| (64) and (6.6), and adding the three quantities obtained gives half of (4.11), 
\ asit should. 

: Now we shall consider the crossings of 6’. The equations in the first part 
+ of the analysis are quite similar to those encountered in Section 3.8 of 
Reference A where the maxima of R, for noise alone, are discussed. We 


} start by introducing the variables x, x2, - -- 26 where 
; x, = I, = Reos 6 — Q, x= I, = Rsin@ (6.7) 
r } and the remaining x’s are defined in terms of the derivatives of J, and J, and 





) are given by the equations just below (3.8-4) of Reference A. 
Here we shall consider the noise band to be symmetrical about the sine 
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wave frequency /, so that 6; and 6; are zero. Then from equations 


and (3.84) of Reference A the probability density of x1, 2, +++ x¢is 
; aL load + 22) + 2be | ) 
exp -_ {Og\X, Tt Uy Tt L£02\X1X3 - aXe 
828 bo DP 2D his 
+ (D/be)(x3 + x5) + bo (x3 + +)1] 
where D = bob; — b: and the b,’s are given by equations (A2-1). R 
ing the x’s by their expressions in terms of R and @, similar to thos 
above equation (3.8-5) of Reference A, shows that the probability densit 
for KR. k’. R', 6. 0, 0” is 
P(R, R’, R"”, 0, 0, 0") = exp {| —~—~ [b,(R° — 2RO cos 6 + O 
82° be D 2D 
(D/ bs) (R® + R*0’*) + 2bo(RR” — R?6’*) 6.9 
bo(R”? — 2RR’0? + 4R°O? + 4RR’0'0” + R20" + R09”) 
— 26.0(R” cos 6 — R6” cos 6 — 2R’6’ sin 6 — RO” sin 0) 
It must be remembered that (6.9) applies only to the symmetrical case i 
which 6; and 6, are zero. 
Integrating R’ and R” in (6.9) from —* to x gives the probabilits 


} 


density of R, 0, 6’, 0’. The integration with respect to R” is simplified by 


changing to the variable R” — R@’. The result is 


p(R, 0, 6’, 0”) = R*(2x)-2(bobeD)“ 21 + u)-? 


exp(- i & — 2ROQcos @ + O + bo R’ 0” ‘be 


(6.10 
(Obe sin 6 + dbo Re") 
(1+ u)D 
where thob)0’” D. The expected number of times per second the tim 


derivative of # increases through the value 6’ is 


No = ao | ar [ d6’’ 6"’ p(R, 0, 0’, 0’’) 


“0 0 


x 
2 e 1/2 
mw (b06/bo) ao [ rdr 
A /0 


r 


“6 611 
[ xdx ! 
0 


exp [—yr? + 2racos 6 — a® — 6(x + asin 6)?| 


. 








whe 


\ 


wh 


cor 


bo 
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ere we have set 


Inna a4 Stone tie Cas 


ry = R(2b,)“ x rb9” | bo 


a= QO(2h) “=p ¥ 1 + 0) 8" /b» rr 6.12 


6= — 
(1+ “)D 
r being regarded as a constant when the variable of integration is change: 
from @” to x. 

he double integral in @ and x occurring in (6.11) is of the same form as 
(A3-1) of Appendix III and hence may be transformed into (A3-3 Here 
1=ra,c = —6a*,c+6=0. The diameter of the path of integration ¢ 
may be chosen so large that the order of integration may be interchanged 
and the integration with respect to r performed. The result is again an 
integral of the form (A3-3) in which a? = 0. When this is reduced to (A3-6 


it becomes 


Ver = € °° (2wy) be “(b9d) 1? [e-"? 277 (6/2) 
(6.13 
+ (1 + yd)(1 + yd/2) "Te \yd(2 + v5), py +6p 2} 
where we have used Je(—k, x) = Je(k, x) which follows from the definition 
A1-1) given in Appendix I. 
When there is no sine wave present, p is zero and (6.13) becomes 
bs bo , 2 
4 1 be - { bo " b + 46 
No - 6.14 


2ry \ bod 7 bs 8 
2r {1 Ty 6 


his gives a partial check on some of the above analysis since (6.14) may be 
obtained immediately by setting a equal to zero in (6.11). Another check 
may be obtained by letting p — =~ and using Je(k, =~) 1 k 
6.13) becomes 


Vor ~ (2r)~"(b4/ be)! 20 F (6.15 


e é . ° ° , , 
which agrees with the result obtained from & ~ J,/(Q. 
For the case in which the power spectrum w(/) of the noise is equal to the 


constant value wo over the frequency band extending from / 8/2 to 
f . =e ? 
fo + 6/2, 
bo = Bwo, be = r-B°w' 9 3 = 7°8"bo 3, bs = 718°wW/5 7*3*b,)/ 5 (6.16) 
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These lead to 


3 = (bo/be)!/20’ = 31/26’ / (xB) D/bs = bybo/bs — 1 = 9/5 — 1 


’ 9 5 


u =4 b32" ia Se 6 = — 
4(1 + 52°) 


y=1+2 


and the coefficient in (6.13) may be simplified by means of 


1 ey i ‘ 4 
2wy \bod 1+2\ 15 ; 


From (6.14) we see that (6.18) is equal to Vg, when noise alone is present 
(and is of constant strength in the band of width 8). The curves of \ 
versus 3 shown in Fig. 7 were obtained by setting (6.17) and (6.18) in (6.13 
Ny /8 approaches e °/(z 4/3) asz > ~. 

When the wandering point T of Fig. 4 passes close to the point 0, 4 
changes rapidly by approximately + and produces a pulse in 6’. In dis 
cussions of frequency modulation @’ is sometimes regarded as a noise voltage 
which is applied to a low pass filter. Although the closer T comes to 0 
the higher the pulse, the area under the pulse will be of the order of x and 
the response of the low pass filter may be calculated approximately. 

That the pulses in 6’ arise in the manner assumed above may be checked 
as follows. We choose a point relatively far out on the curve for p = 5 in 


Fig. 7, say 2 = v/ 30’/(8r) = 1.6 or 0’ = 2.98. The number of pulses per 
second having peaks higher than 2.98 is roughly Ng, = .0098, and half of 


these have peaks greater than 6’ = 3.88 which is obtained from Fig. 7 ior 
Ng = .00458. From Fig. 6 we see that 6’ exceeds 2.98 about .0018 of th 
time. Thus the average width at the height 2.96 of the class of pulses 
whose peaks exceed this value is .0018/(.0098) = .2/8 seconds. This figure 
is to be checked by the width obtained from the assumption that the typical 
pulse arises when T moves along a straight line with speed v and passes 
within a distance 6 of O. We take tan 06 = ot/b = at so that 
0° = a/(1 + a’). From this expression for 6’ it follows that a pulse of 
peak height 3.86 (the median height) has a width of .3/8 seconds at @ 

2.98. This agreement seems to be fairly good in view of the roughness of 
our work. A similar comoarison may be made for p = 0 by using the 
limiting forms of (5.15) and (6.18). Here it is possible to compute the 
average width instead of estimating it from the median peak value. Exact 
agreement is obtained, both methods leading to an average width of 2, (44 
seconds at height 6’. 
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Fig. 7—Crossings of time derivative of phase angle. 


= expected number of times per second d@/dt increases through the value @’. p 


and the power spectrum of Jy are the same as in Fig. 5 
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dé 
7. CORRELATION FUNCTION FOR 
( 

In this section we shall compute the correlation function (7) of 
We are primarily interested in Q(r) because it is, according to a fundamy 
result due to Wiener, the Fourier transform" of the power spectrum 
of 6’(t). 

We shall first consider the case in which the sine wave power is very 
compared with the noise power so that, from (3.6), @ is approxima o—_ 


; , , a ‘ape equati 
I, (and @ approximately 7,/Q. Then using (A2~-3) and (A2-1) 4: 


conver 
We: 


a (A (t+ 7) = OPI t + 7) Re: 
of 4’(2) 


nw 


—g"0° = 49r O° | w(f)(f — fy)’ cos 2x(f — fy)r df 
Ay) ‘ : functi 


it wou 


oe ee ’ difficul 
When w(f) is effectively zero outside a relatively narrow band in the neig aii 
or a A . ; 7 a enter: 
borhood of f,, as it is in the cases with which we shall deal, (7.1) leads to t! Siesmes 
“ Cau! 


relation (divide the interval (0, ~) into (0, f,) and (f,, ~), introduce new 


a deal w 

variables of integration f/f, fo — f, fe = f — fq in the respective interva emai 
Ai d ' 

replace the upper limit /, of the first integral by ~ , combine the integra 


and compare with (2.1-6) of Reference A) 
Power spectrum of 6’(¢) Wf) 


sideret 


We 
= 4? PO[w(fa + f) + wife — Sf] 
This form is closely related to results customarily used in frequency modu 


tion studies. It should be remembered that in (7.2) it is assumed t! 
O0<f<f,andrmsly <0. 


Additional terms in the approximation for Q(r) may be obtained by and tl 


expanding 


8 


0+ 1. 


in descending powers of (), multiplying two such series (one for time ¢ and 


6 = arc tan 


the other for time ¢ + 7) together, and averaging over t. If [.1, 74; and where 
I 2, 1.2 denote the values of /,, 7, at times ¢ and ¢ + 7 respectively, th As be 
average values of the products of the 7’s may be obtained by expanding t!. Phe r 
characteristic function (obtainable from equation (7.5) given below | F sion | 
setting 3; = s5 = s7 = ss = 0) of the four random variables J .1, J s1, [¢2, / this n 
This method is explained in Section 4.10 of Reference A. When w(/ 

symmetrical about f, it is found that 


; : ae : i tele proce 
18 The form which we shall use is given by equation (2.1-5) of Reference A. 





PROPERTIES OF SINE WAVE PLUS NOISI 


From the exact expression for Q(7) obtained below it is seen that the last 
equation in (7.3) is really asymptotic in character and the series does not 

mnverge. We infer that this is also true for the first equation of (7.3 

We shall now obtain the exact expression for the correlation function Q(+ 
of 6’(t) when f, isat the center of a symmetrical band of noise. At first sight 
it would appear that the easiest procedure is to calculate the correlation 
function for 6(¢) and then obtain Q(r) by differentiating twice. However, 
difficulties present themselves in getting @ outside the range —7, 7 
enters the expressions only as the argument of trigonometrical functions 
Because I could not see any way to overcome this difficulty I was forced to 
deal with 6’ directly. Unfortunately this increases the complexity since 
now the distribution of the time derivatives of 7, and J, also must be con 
sidered. 

We have 


tan @ 


(O + ] ip 
70(0 + I,): O+ 7 


and the value of 6’(t)@’(¢ + 7) is the eight-fold integral 


+20 o +x 


[ata fo dtepita, 


x . x 


, , 


(0 T I )1 a a ae 
(0+ Ia? +h 


where p(I.1, +++ , I,2) is an eight-dimensional normal probability density 


\s before, the subscripts 1 and 2 refer to times ¢ and ¢ + 7, respectively. 
lhe most direct way of evaluating the integral (7.4) is to insert the expres- 
sion for p(Te1,--* , I’») and then proceed with the integration. Indeed, 
this method was used the first time the integral (7.4) was evaluated. Later 
it was found that the algebra could be simplified by representing p(/ 

I's) as the Fourier transform of its characteristic function. The second 
procedure will be followed here. 
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1 ° ° e ° * , , , i 
Ihe characteristic function for J¢1, Do, 201, Dea) Let) Lea, La, Tee is, | 


(A2-2) and (A2-3) of Appendix II and Section 2.9 of Reference A, 


of , , ; , : 7 
ave. exp Uzi oe t sol 9 Sal +1 a sal oo . S5l e1 a. 6l 2 + S71 51 T Zsl 


exp (= +) [bo(zi - 2 4 23 + 34) r bo(z5 r 26 


)\ol 


+ 2b, (5137 


3538 —— 2627 ) a 


Since we have included },, h, h’, hk” this holds when f, is not necessaril\ 
the center of the noise band. However, henceforth we return to our assu 
tion that f, is placed at the center of a symmetrical noise band and 
by, h, h’, h” to be zero. 

The probability density of 7.1, --- I. which is to be placed in (7.4 


the eight-fold integral 


pla, -*° 50) Tod ee fim 


exp [—itiTe — +++ — isgli2] X [ch.f.] 


where ‘‘ch.f.’’ denotes the characteristic function obtained by setting / 
h, h’, h” equal to zero on the right hand side of (7.5). 
The integral (7.4) for (7) may be written as 


Q(7r) = J,-—-Jeo—J3+ J 


where J, is the 16-fold integral 


a [ata oh | ” dT's (2m) . | dz, ++: 


0 x v— 


, 


exp [—1237¢ se I i2g I 2] 
(O + Ta)(Q + Teo) Tin Tee 
(0+ Ta)? + TAQ + Tea)? + Ti 
and Jo, J3, J, are obtained from J; by replacing the product (QO + / 
, ; / y/ 
(O + Te)Lidia by (O + Teileel ee, Teada(Q + Teale, Ted al alc 
respectively. 
The integration with respect to 7.; and J,; in (7.8) may be performed at 
once. Wereplace()+ /.,and/,,by xand y, respectively, and use 


X [ch.f.] 


a +00 . 


a 9 — 2riz 
dx dy é oF a 
- = 


x? + y? 2 + ¢ 


J— 
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Phe integration with respect to 7,2 and 7,2. may be performed in a similar 
nanner. In this way we obtain a 12-fold integral. 
The integrations with respect to the /’’s may be performed by using 


+ a +2 





a ae .; : 
| dl e '*! f(s) dz f(Q) 
dn J % vx 
- 7.10 
1 . x » +s if >) 
} 1 , ‘ ( \ 
I dl C "#(s) dz = —}] : 
2m Jc A dz _jamo 
The result is the four-fold integral 
- x +2 @ loa g > ) 
2 sad ics ~ a4 
J, (27) dz, - Gites 
00 J— x (21 -- 23)(Z2 + 24) (7.11 
exp [— (00/2)(1 + 22 + 23 + 33) 9(31S2 + 3334) + 10(21 + 3e)]. 
i In the same way J, J3, /4 may be reduced to the integrals obtained from 
| (7.11) by replacing 2122(g” — g’°ssz4) by —g “2124, —g 3923 and 23z4(g” 
+) 1S to . ’ ’ . ‘ ai 
—g’*s 132), respectively. When the J’s are combined in accordance with 
(7.7) we obtain an integral which may be obtained from (7.11) by replacing 
” "9 ‘ 
S\e\g — g °Z324) by 
1.6 g” (S152 + 3334) + g’2(s1%4 — 2223)? (7.12) 
The terms 21 + 23 and s2 + 23 in the denominator may be represented as 
bh, € infinite integrals. Interchanging the order of integration and expressing 
(7.12) in terms of partial derivatives of an exponential function leads to the 
six-fold integral 
dad x 12 fs) af eto ~ +00 
id 2 2 € ; 
Q(r) = (4) [ du | dv| —g’-+g° — dz, ef dz, 
’ “0 “0 Og Ja‘ a=) /—x » 
2 2 2 2 (7.13) 
exp [— (bo + u)(si + 33)/2 — (bo + v)(s2 + 24)/2 
- 2( 2152 T- 2324) — (2124 — ZeS3) 10 (24 T 22)} 
8 where the subscript a = 0 indicates that a is to be set equal to zero after 
the differentiations are performed. 
When the four-fold integral in the z’s is evaluated (7.13) becomes 
aw ae 4 at 
/ € yo € 
Q(7r) = | du | dv| —g’ —+2°- , 
/ 0 /0 Og Oa" ja=0 
exp [—(Q°(2b) — 2g + u + 0)/(2D)} (7.14) 
4D } t 5 


a® 4) 


9 = i du | do|(g” — gg”’)(2 — 2F + 07/g) — gO" /gle”/(4D35) 
0 0 
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where 


D = (bo + 4)(0o + 2) g — a, F = O'(2by) — 2g + ut 2)/(21 

and Dy denotes the value of D obtained by setting a = 0. When differ 

tiating with respect to @ it is helpful to note that 

a f(D) vr oD ’ ee aD 
—> {'"(D)\ - + f(D) — 
Oa 0a : Oa 

and that only f(D) = df dD need be obtained since 0D da vanishes whe 


a= Q. 
In order to reduce the double integral to a single integral we make t 


change of variables 
O'(by + u — g) 
2[(bo + 2)(bo + 0) — g?| 


§ = O"(bo — : = g) (2D ), F = 7 + Ss 


r O" (bo + ae 2) (2Dpo) = 


(7.15 


O(r, s)/A(u, v) = —rs/Do, 4srDy = O'[(0" — 2g(r + 5)] 


The limits of integration for r and s are obtained by noting that the points 
(0, 0), (~, O), (x, ~), (O, ~) in the (u, v) plane go into (Q?/(2b9 + 2¢ 
0?/(2bo + 2g)), (07/(2b9), 0), (0, 0) (0, 0?/(2bo)), respectively, in the (r, s 
plane. It may be verified that the region of integration in the (r, s) plane 
is the interior of the quadrilateral obtained by joining the above points by 
straight lines. Equation (7.14) may now be written as 


9 ” 9 9 9 \ 
fi(g° — eg )\(2—2r7 — 2s + 0'/2) — g°O'/g| _-- 
Q(7) [| 4\6 : ulb 6 b\e *drds 


‘ QO7(Q? — 2g(r + s)] : 
12 ” ‘ (7.16 
ie: | a eZ 
a ro2 v1 ~ 572 92 
2g 2¢° 
where y; and y2 are the dimensionless quantities 
PO () — Dy —= Ie 2/9 
22 (2 2r 2s + O'/g) _ 
V1 | | sera — Be dr ds 
: O7[|0? — 2e(r + s)] 
/ / 2ee." “* dr ds 
‘ J QO? — 2g(7 +s) 
It is seen that 
a, \ 
V1 2g0 2 ve + | e” °dr ds}. (7.17 


Since the integrands are functions of r + s alone we are led to apply the 
transformation 


ao 


- 8 9 
‘ Ls = BW) x 
| f(r + s) drds = uf(u) du + [ f(u) du (7.18 
0 %°a@ 


28 — a 


vd a 


A 








an 








i 
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vere A is the area enclosed by the quadrilateral whose vertices are at the 


points (r, s) given by (0, 0), (O, a), (3, 8), (a, 0) and it is assumed that 8 and 
are positive. «is a new variable and is not the one introduced in (7.13 
Setting a = 0", (269) and 8 = (” (269 + 2g), using (7.18), and introduc- 


ng the notation 


p = O° /(2bp), k = g/b 

Or do /.19) 
3 O*/(2¢) = p k. Xd = c é 
s x § | a g 1 + } 


permits us to write 


ef p a 
alee J : oe 
| | e” “drds = / ue du + | edu 
J i >» AS 
A 


0 . p 
7.20 
1— he* | pe ; 
Lup’ beg 
and (7.17) yields 
y 1+k -» 
p Z , = 
piece a: 1 =_— é / 21 
Vo k Vi ot i - p i = j 


where we have expressed J in terms of p and k, 

The double integral detining ve may be treated in the same way as (7.20 
of r—s ep u aX / u 7 
e€ dr ds ue du, ff piA— we du 


JJ €é—r—s Jo €—u J, (XN— phE—U4 


— n)in the two numera- 


See 


Writing w= — — (E— uw) andA— uw =A— ESF | 
tors leads to 


pte du ‘oe 
Vo =§ - — edu 
‘ 09 &— 4 Jo 


y a 
bee 
ad 4 " 
7 é du p } 
cometay . | é au 
Jp €—U A.— 2p 
where we have used p(A — £)/(A — p) = —é to simplify the coefficient of 


the third integral. When the second and fourth integrals are evaluated, 
their contribution to ye is found to be equal to the terms independent of y, 


on the right of (7.21). Hence, comparison of equations (7.21) and (7.22) 
shows that 
af - aX ? 
ee du e du (= 49% 
y= — CS) 


y &E— UU Jp E— th 
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The integrals in (7.23) may be evaluated in terms of the exponential i 


gral £2(x) defined by, for x real, 


0 n 
x 


Fi(x) e'dt/t=C+ * log. a+ 
~e > n!/x"*! 


where C = .577 -+- is Euler’s constant and Cauchy’s principal value of the 
integral is to be taken when x > 0. Weset ¢= & — wand obtain 


nml 2!” 


ok - Be vy fala R) | 
n= \ Eile k] — 2Eilp(1 — k)/k] + Bile + Al 


where we have again expressed £ and \ in terms of p and k. 

A power series for y, which converges when —1/3 < k < 1 may be ob- 
tained by expanding the denominators of the integrands in (7.23) in powers 
of u/é and integrating termwise: 
yi = El — 2e? + | 

+ 171 — Wi + p/i)e? + (1 + A/1Ne™| ( 
+ QE — 201 + p/1!+ p2/2e? + (A+ A/1E + A2/2WNe 


+ 


~! 
ei 


The following special values may be obtained from the equation given 
above. When p = 0 
yi = —log, (1 — ?) 
(7.25) 
Ve == § 
rhis result may also be obtained by evaluating the integral obtained when 
we set O 0, 2, r, COS 01, 23 = 71 SiN 0), 32 = 12 COS Oo, 34 = re SiN be 
in (7.11) and (7.12). 


Near k = 1, 
yi = e “[Ei(p) — C — log, p(i — ?)] 


(7.26) 
ye = pyi- 1+ (1+ p)e” 
Near k = 0, 
yi = RL — &*)*/p, yo= VN (7.27) 
except when p = 0 in which case ¥y, is approximately k?, 
When p is large 
k, 1k’ , 21k, 31k 
mr itt ttt. 
p p° p p " 
| (7.28) 
p 1!k 2!k° 
Yq ™ —1 a yy Ow + my + 
k p p° 
except near k = 1 where both y; and ye have logarithmic infinities. The 


asymptotic expansion (7.3) for Q(r), which was obtained by the first method 








al 
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bi 
0! 
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| 26) 
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r Q(r) in terms of y; and yp. 
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this section, may be checked by inserting (7.28) in the expression 





7.16 


Values of y; and y2 tabulated as functions of & for various values of p are 


viven in Table 3. 


VALUES OF ¥; AND yo USED I 


90 


QO? 
Q4 
96 


Q7 


98 
99 
995 


9907 


2 


0 
.01005 
04082 

09431 


.1744 
.2877 
4403 


6733 


.0216 
.2228 
.4890 
.6607 


1.8734 
2.1507 
2.5459 
2.8285 


3.2289 
3.9170 
4.6072 
J.4ht0 


TABLE 3 


tT) = O(a (t - 7) = lg*(y — 9 gg 22° 
'd 3 
g(r) = (f cos 2r(f f rt df k = g(r)/g(0 
“U 
Values of 1 Value 
1 
0 0 0 0 0 0 

03526 .04224 .03854 .02000 .03171 .04147 
OS043 09003 .07979 .04105 06550  .O8654 
1379 1452 1246 06292 .1022 1363 
2110 2102 .1740 O8586 .1432 1926 
3056 2886 .2296 1101 1914 2579 
4278 3860 2942 1258 2481 3368 
5953 5129 3721 1636 3220 4379 
8416 6914 $729 1941 275 5803 
9798 7888 5242 2075 4866 6593 
1.1590 .9127 5866 2219 5641 7619 
1.2742 .9898 6241 2296 6138 8260 
1.4144 1.0834 6686 2378 6753 9058 
1.5948 1.2024 .7217 2466 7550 1.0093 
1.8486 1.3668 7939 2566 8711 1.1558 

2.0251 | 1.4815 | .8414 2623 9474 1.2605 * 
2.2762 | 1.6405 | .9073 2690 1.0704 1.4081 
2.7080 1.9066 1.0127 2778 1.2773 1.6610 
3.1341 | 2.1721 |1.1125 | .2846 |1.4838 (1.9175 
3.4445 2.3622 1.1866 2889 1.6367 2.1048 


03936 


O827 
1315 


1870 
2515 
3289 
$269 


5602 
6318 
7226 


&9? 


S486 
9333 
0546 
13606 


5 


Negative values of k have not been considered since they 
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0 
02051 
O4283 
06702 


QO384 
1238 
1576 
1975 


2461 


2693 


are not required for the case in which Jy has a normal law power spectrum, 
the case discussed in the next section. 


? dé... ‘ 
8. POWER SPECTRUM OF it WHEN Jy HAS NORMAL LAw PowWER SPECTRUM 
at 


The problem of computing the power spectrum |W(f) of 6’(t) appears to 


be a difficult one.* 


of work we have had to do two things which are rather restrictive. 


In order to obtain an answer without an excessive amount 
First, 


we confine our attention to the case in which the power spectrum w(/) of 


*Since the above was written the general f. m. problem has been studied by D. Middle- 


ton. He generalizes our (7.11) and (7.12), 
integrand in powers of g, and integrates termwise. 
a.m. theory. 


W 


introduces polar coordinates, 
f) then follows somewhat 


expands the 


as In 
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Ty is of the normal law type (our method could be applied to other t 
but g’ and g” would be more complicated functions of 7 and Table 3 w 
have to be extended to negative values of k, if they should occur). Sec 
we resort to numerical integration to obtain a portion of HW(f). Be 
of the second item our results are either tabulated or are given as cu 
shown in Figs. 8 and 9, except when Q = O (noise only) in which case 
power spectrum of @’ is given by the series (8.7). 
The power spectrum of /y is assumed to be 


Wo f—f.)2/(%2 


o(f) e 8.1 
oV 2 

The mean square value of /y is equal to that of a noise current whose power 
spectrum has the constant value of fo/ (o/ 27) over a band of width /; 

oV le = 02.507. The value of w(f) is one quarter of its mid-band valu 
at the points f — fg tov/ 2 log, 4 = +01.665 (the 6 db points) and the 
distance between these points is 3.330¢. Integration of (8.1) shows that t 
mean square value of Jy is Yo in accordance with our customary notati 
The mid-band value of w(f) is fo/(or+v/ 2m). 

Assuming f, > o and evaluating the integrals (A2-1) of Appendix I] 
detining 6») and g gives 


b= vo, g= poe” = Pee 


ce : , ) ” fc 2. 2 en 2 
PE uu £T7ou, : a i lla (2rc) "(1 u) (29 
‘AA 49 : a 

5 es 2 42/2 
—(2xo), k= 2/6, = €" 
e 
where we have set 

u 2roT, u'i= 2xo (8.3 


and the primes on g and u denote differentiation with respect to 7. The cor- 


relation function is accordingly, from (7.16). 
Q(r) = 2n*o*(v1 — u*ye) (8.4 
If 6’(¢) be regarded as a noise current its power spectrum is 
a 


W(f) =4 | Q(t) cos 2afs dr (8.5 
“0 


Ti 


When noise alone is present, p is zero and (7.25) yields 


‘ 
Q(r) = —27°o? log, (1 — k?) = —2r’o? log, (1 — e™ ) (8.( 


i 
; 
t 
1 
{ 
i 
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t) In this case the power spectrum is, from (8.3), 
Wi 
7 
< 
Wy(f) = —4mo | cos(uf/o) log (1 — © **) du 
ery J 
UT (8.4) 


S x 
, ) : > i 
¢ £07 i ¢ 
n=l 


the series being obtained by expanding the logarithm and integrating term- 
wise. When this equation was used for computation it was found conven 
ient to apply the Euler summation formula to sum the terms in the series 
beyond the (N 1)st. Writing 6 for f*/ (407), the series in (8.7) becomes 


we b)'* erf [(b/N)'] + A 4 ( : =] 
d the + (7/0 er yi lV) ee ae é D 12\ > \ (8.8) 


tir 1 ( — 105, 105 6 21 6 b ) 
Y 790N 8 tN 2 N2:' N 


When 6 is zero the sum" of the series is 2.61237 --- . The values for p = 0 
in Table 4+ were computed by taking .\ 12 in (8.8). As b— « the domi 
nant term in (8.8) is seen to be the one containing erf (choose .V so that 
b = N22), Hence as / > 


Wey(f) ~ 41°07/f. (8.9) 


When both noise and the sine wave are present it is convenient to split the 
power spectrum into three parts. The first part, 117,(/), is proportional to 
Wy(f), the power spectrum with noise alone. The second part W2(f) is 
proportional to the form H/(f/) assumes when rms /y «< ( and the third 
part W’;(/) is of the nature of a correction term. This procedure is suggested 
when we subtract the leading terms in the expressions (7.26) and (7.27) 
(corresponding to k land k 0, respectively) from y;. Likewise we 
8.4 subtract the leading term in ye, (7.27), at k = 0 but do not bother to do so 
at the end k 1 because wy. approaches zero there. We therefore write 
i UVo = [yy t+ e ‘log (1 — hk?) —k(1—e*)/p—u Vo 

+ wk(1 — e °)?/p] —e * log 1 — &) + 1 — w)R(I e °)*/p (8.10) 


wid 

Y . 9 2 |Zma) 2 

= Z(u) — e “log (1 — k&) — > ee") 
bop 


4% “Theory and Application of Infinite Series,” Knopp, (1928), page 561. 
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where 7(u) denotes the function enclosed by the brackets in the first eq 
tion and the expressions for g”/g and k in (8.2) have been used in the repla 


























transform (8.5) leads to 


4iro | Z(u) cos (uf/o) du 
0 


ment of (1 u*)k, 
TABLE 4 
VALUES OF W;(f)/(42°o) 
p 1.0 2.0 0 

er 

0 0 03517 .03891 02444 — .OO1948 
1 0 03003 .03196 01830 — OO1814 
2 0 .01717 01486 003304 004052 
3 0 — .002436 .004014 .01252 .008225 
4 0 .008757 .01730 .02244 .01027 
6 0 .01478 .02157 .02167 .007665 
8 0 01018 .01366 £01237 003505 
10 0 005768 .007378 .006201 001437 
2 0 004027 004463 003552 0006439 

Vatures or W(f)/(42°o0 

0 7369 4118 2322 07529 003017 
1 7098 $294 2672 .1134 .02342 
2 6439 $516 3231 1784 05828 
3 5542 $225 3225 . 1947 .06852 
4 $623 3496 2654 .1580 01590 
6 3195 2178 1508 .07554 .01540 
8 2390 .1553 1019 04506 ~005325 
10 1908 1215 .07768 .03206 .002726 
2 1595 . 1003 .06306 .02511 001719 


Inserting (8.10) in the expression (8.4) for Q(7) and taking the Fourier 


W(f) Wilf) + Wal(f) + Ws(f) 
Wi(f) = e *Wr(f) 
—_— > p 2 = 
W2(f) ms. lie. | g”’ cos 2nfr dr 
. ~ s (8.11) 
(1 — ¢*) er eo 
e') (2xf)?! 
ov 2r 
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these equations Wy(f) is obtained from (8.7), and W2(f) by two-fold 
tegration by parts to reduce g” to g then evaluating the integral obtained 
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Fig. 8—Power spectrum of d@/dl 


Power spectrum of Jy is assumed to be 


Yolo 2r exp |[—(f — f,)?/(2e?)]. 

(ov ) )-1 xp | (f fy (2 ) 

In this expression fis a frequency near f The f in W(f) and in the abscissa is a mucl 
t d i 7% d 

wer frequency. W(f power spectrum of 6’ dé/dt, & being regarded as a rando 


noise current. Dimensions of W(f) df same as (d6@/dt)* or (radians)?/sec.? 


by substituting the expression (8.2) for g. That W(f) approaches W2(/) 
as p — * follows when expression (8.11) for W2(/) is compared with the 


limiting form (8.13) given below. 
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Insiead of dealing with W(f) it is more convenient to deal with (47 

































I(/) which is the sum of the three components 
e7P x 
l . o f2)(4n@2 
leo) Wf) =a Yn e : 
LV TF n=1 
(1 — CPP LEY’ _ 2, cre 
irc) WoC f) é | 
pv 2r o 
as 
1 : ’ , 
tro) Ws(f) Z(u) cos (uf/o) du 
T 0 
0.08 - , “3 
| | | | 
0.07 | | 4 | | | | | 
| | | 
0.06 —\}——_}—_}—_— + | | 
] | | | | 
a | 
te ee 2 | 4 } + r + —+ 
e | | | 
ke | | } | | 
ate, Eee EE + + + + ——— —— + + 
C | | | | 
= oe oe | | 
} Nee | 4 | 
EXACT FORM | | 
FOR p=5 | 
_LIMITING FORM | 
~ WITH 2&5 | | | | | | 
tes tas | t t +—4 
| aa | | | | | 
)6-DECIBEL} | 
| POINT | : 
3 5 6 
f/o 
| Fig. 9—Approach of W'(f) to limiting form. 
As p— x, W(f) > 49°o (pV 2r)* (f/a)? exp [ — f2/(20?)]. 
Phe integral involving Z(u) has been computed by Simpson’s rule, y; and 
v» being obtained from Table 3, with the results shown in the first sectior 
of Table 4. The value of 1:(f) may be computed directly, and W ,(f) ma 
be obtained from H’y(f). The values of these two functions together wi 
those of Is (f) enable us to compute the values of (42°0)~'IW’(f) given 
i Table 4+ and plotted in Fig. 8. 
i Since, as is shown by (8.9), Wy(f) varies as 1/f for large values of f, the 


areas under the curves of Fig. 8 become infinite. This agrees with the fa 
that the mean square value of 4’ is infinite. 

The values of (42°o)~'I(0) for p equal to 0, .5, 1, 2, and 5 are .7369, 4118 
7299 


2044 


07529, and .003017 respectively. When these values are plotted « 


| 
{ 
| 





whe 
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emi-log paper they tend to lie on a straight line whose slope suggests 


()) decreases as e “ when p becomes large 


Phe limiting form assumed by IV'(f) as p—> ~ is given by equatio 


1 
When the normal law expression (8.1) assumed in this section for the power 
spectrum of Jy is put in (7.2) we find that 
dro } 9 
W(f)- P Qe 12 


pv 27 \o 


Fig. 9 shows that for p 5 the limiting form (8.13) agrees quit 
the exact form computed above. 
Both (7.2) and (8.13) show that, for small values of /, the power spect: 
7 of 6’ varies as /* when p >> 1. This is in accord with Crosby’s* result 
that the voltage spectrum of the random noise in the output of a frequency 


modulation receiver is triangular when the carrier to noise ratio is large 
When this ratio becomes small hé tinds that the spectrum becomes re¢ 
tangular. Tig. 8 shows this effect in that the areas under the curves betwe 
| the ordinates at / Q and f Ao (where A is some number, generally less 
than unity, depending on the ratio of the widths of the i.f. and audio bands 


| . 

become rectangles, approximately, as p decreases. 

ni ; 

| APPENDIX | 

3 | 

| 

| 

| THE INTEGRAL Te (Rk, 2 
The integral!® 


Te(kx) e Io(ku) du \1-1 
0 


where /9(ku) denotes the Bessel function of imaginary argument and order 


zero, occurs in Sections 2 and 6. The following special cases are of interest 
Te(O, x) 1 é 
. Te(1, Xx) Xe ve \ T [ya \1-2 


‘ Te(k, «) 


lhe second of these relations is due to Bennett.” 


*M. G. Crosby, ‘* Frequency Modulation Noise Characteristics,” Proc. I. R. E. \ 5 


1937), 472-514. See also J. R. Carson and T. C. Fry,‘ Variable Electric Circuit Thi 
he vith Application to the Theory of Frequency Modulation,” B.S.T.J. Vol. 16 (193 
513-540. 
6 The notation was chosen to agree with that used by Bateman and Arc] G 
to Tables of Bessel Functions appearing in “Math. Tables and Aids to ( \ 
g 1944) pp. 205-308) to discuss integrals used by Schwarz (page 248 


"It is given in equation (62) of the referen cited in cont t 
1.2) in Section 1 
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The values in the table given below were computed by Simpson’s ruk 


numerical integration. 


where 


where ¢; = 


Tek, x) ~ (1 


le(k 


= 


» &) 


a. [ cous 


) 


ik (k 9) 


n=) 


(2n)! 


A, 


nin, 


ol : 
' e 
(2n)! 


When x is so large that /e(k, «) is nearly equal to Je(k, « ) we have 


+- 


aD 


Rk)? — [2k — Rk) *(2// we) e' dt 


V «(1 —k). 


However, this was not found to be especially usefu 


in checking the values given in the table. 


tN 


mn 


. 


oo) 


CD keh 


1 
329 
1 


8 
2 
5 
5 


SV he ee 
“TN 1 Ge 


550 
6321 
6988 
.7534 
7981 
8347 


8647 
.8892 
9093 
.9257 
9392 


9502 
959? 
9666 
.9727 


.9776 


.9817 
9830 
.9877 
9899 
9918 


9933 
9955 


— ne 


TABLE OF Je(k, x) 
2 4 
0 0 
.1813 .1814 
3298 3303 
24517 .4530 
.5516 .5545 
6337 6386 
7012 .7086 
.7567 7669 
.8025 .8157 
8401 8566 
8712 8910 
8968 9201 
.9179 9446 
9354 9655 
9499 9831 
9618 9982 
9718 1.0110 
9800 1.0220 
9868 1.0314 
9925 1.0394 
9971 1.0463 
0010 1.0522 
0043 1.0574 
.0070 1.0619 
0092 1.0657 
0111 1.0090 
0140 1.0743 


1 
1 


.9255 
.9607 
.9916 
0186 
.0424 


0635 
.O822 
0988 
.1136 
.1208 


1386 
. 1492 
.1587 
.1672 
1749 


.1818 
.1937 


ee ee ee ee pt pet ed pe 


.1816 
.3322 
.4586 
5661 


.6584 
7386 
8089 
.8712 
.9267 


9766 
0217 
0627 
. 1001 
.1345 


1061 
.1953 
2223 
. 2475 


2708 


. 2926 
3130 
3320 

3499 
3606 


3823 


4110 


“ty 


ae a ee ped peed feed bed gee 


The work was checked at several points by usit 


1.0 


0 


.1818 
3337 
.4629 


.5749 


67360 
.7620 
.8422 
.9157 
9839 


0476 
.1075 
.1642 
.2183 
. 2699 


3195 


aah 
. 3672 


4132 


.4578 


/ 


5010 


5430 
. 5839 
.6237 
.6625 
. 7005 


7376 


8095 








aoauw 
me oo 


oo a 


aes) 
= => 


10.0 


12.0 
13.0 
14. 
15.¢ 


are ¢ 
this 
and 
In n 
thes 
cent 
The 
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PaBLE—Continued 


0 ey 4 6 § 9 1 
5.8 .9970 1.0160 1.0783 1.2034 1.4364 1.8786 
6.2 .9980 1.0174 1.0814 1.2114 1.4590 1.9452 
6.6 .9986 1.0183 1.0837 1.2180 1.4792 2.0097 
7.0 .9991 1.0190 1.0854 1.2234 1.4972 2.0722 
7.4 .9994 1.0195 1.0867 1.2278 1.5134 2.1328 
7.8 .9996 1.0198 1.0876 1.2375 1.5279 2.1917 
8.2 .9997 1.0201 1.0885 1.2346 1.5409 2.2491 
8.6 9998 1.0202 1.0891 1.2371 1.5526 2.3050 
9.0 .9999 1.0203 1.0896 1.2393 1.5631 2.3597 
10.0 1.0000 1.0205 1.0902 1.2431 1.5852 1.9207 2.4910 
11.0 1.0000 1.0206 1.0907 1.2456 1.6024 1.9668 2.6157 
12.0 1.0000 1.0206 1.0909 1.2471 1.6158 2.0066 2.7347 
13.0 1.0000 1.0206 1.0910 1.2482 1.6263 2.0411 2.8487 
14.0 1.0000 1.0206 1.0910 1.2488 1.6346 2.0711 2.9584 
15.0 1.0000 1.0206 1.0911 1.2492 1.6412 2.0973 3.0641 

2 1.0000 1.0206 1.0911 1.2500 1.6667 2.2942 L 

x 
als 90 96 1. 

15.0 1.8773 2.0973 2.5810 3.0641 

16.0 1.8899 2.1201 2.637 3.1663 

17.0 1.9006 2.1403 2.6894 3.2653 

18.0 1.9095 2.1579 2.7381 3.3614 

19.0 1.9171 2.1737 2.7837 3.4548 

20.0 1.9235 2.1870 2.8263 3.5457 

0 1.9597 2.2942 3.5714 £ 


APPENDIX II 
SECOND MOMENTS ASSOCIATED WITH I, AND I, 


The in-phase and quadrature components of the noise current /, 
M 
Li) = Pm Cn COS [(wn ae OE Onl 
n=l 


M 
I,(t) = Do en sin [(wn — g)t — eal 

n=l 
are Closely related to the envelope R and phase angle @ of the total current, 
this relationship being being shown by the equations (3.4) and (3.5). /,(t 
and J,(¢) and their time derivatives may be regarded as random variables 
[In much of our work we have to deal with the probability distribution of 
these random variables. By virtue of the representation (3.3) and the 
central limit theorem’® this distribution is normal in the several variables 
The coefficients in the quadratic form occurring in the exponent are deter 


Section 2.10 of Reference A. 
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mined by the second moments of the variables.” Here we state these 
moments. Some of the moments have already been given in Sections 3.7 
and 3.8 of Reference A. For the sake of completeness we shall also give 
them here. The new results given below are derived in much the same wa\ 
as those given in Reference A. 

Let 


@ 


by = (2m) [ w6A(F- fa)” df 
0 


@ 


y= | wf) df = fo 
0 
(A2-1 


2 


g= | w(f) cos nf — fr di 


in | wif) dn 20lf ~ fed 
0 


and let g’, g”, h’, h” denote the first and second derivatives of g and h/ wit! 
respect to 7. For example, 
pw 


gl = —2e | w(f)(f — fa) sin 2x(f — far df 


Incidentally, in many of our cases w(f) is assumed to be symmetrical about 


fq. This introduces considerable simplification because },, b3, 05, --: 
h, h’, h” , reduce to zero. 


The following table gives values of },’s and g for two cases of frequent 


occurrence 
Ideal band pass filter Normal law filter 
centered on f, centered on f,, fg > o 
w(f) ws for fa < f < fo 7 Yo ge SFa)?/20 
— and zero elsewhere ov 2m 
bo wo( fs — fa) Wo 
bs rwo( fe — fa)°/3 drop 
bs two(fe — fa)®/5 48r'toty 
g (rr) wo sin t(fo — fa)t ow" 


If we write ,, 1., I for I(t), I.(é), 12), where the primes denote differ- 


1 Section 2.9 of Reference A. 
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an 
uw 


entiation with respect to ¢, and do the same for /,(t) and its derivatives we 
have, from Section 3.8 of Reference A, 


eS¢ 


> 


ivi =[=b, Idt,=0 
Il,=—-Il,=b, Ide = 1d, =0 
WeaT2=-10=-1 NM =b, T= 10 = 10 = 0 
ll,=-Il,=b, If, =I, = 0 
I2=F%=h 1,1, =0 (A2-2) 
4 When we deal with moments in which the arguments of the two variables 
are separated by an interval r as in (see the last of equations (3.7-11) of 
Reference A) 
I (t)I,(¢ + 7) = fh, 
it is convenient to denote the argument ¢ by the subscript 1 and the argu- 
rit ment!-+ 7 by 2. Then our example becomes 
Tele =h 
We shall need the following moments of this type. 
Teil e2 = Lala = g; Teil «2 - —T eel 1 =h 
" —_— -_- aa / r , ‘ , 
Lele = Inala = —Tele = —Iale =g (A2-3) 
ae 7 a ae ee ' ie , 
Tels _ T col a1 _ —Teul a a —Teol a1 =h 
er 
Tele = Ils = —8", Tala = —Talay = —h" 

It should be remembered that in these equations the primes on the /’s 
denote differentiation with respect to ¢ while the primes on g and / denote 
differentiation with respect to r. 

APPENDIX III 
EVALUATION OF A MULTIPLE INTEGRAL 

Several multiple integrals encountered during the preparation of this 
paper were initially evaluated by the following procedure. The integral 
was first converted into a multiple series by expanding a portion of the inte- 
grand and integrating termwise. It was found possible to sum these series 
when one of the factorials in the denominator was represented as a contour 

- integral. This reduced the multiple integral to a contour integral and some- 
. times the latter could be evaluated. 
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We shall illustrate this procedure by examining the integral 


7s 
I dé | dx x exp | - x + 2a cos 6 + 26x sin 6+ ¢ sin® (A3-1 
Expanding that part of the exponential which contains the trigonometrica 
terms and integrating termwise gives 


x 00 x 


ab cla T(t + m+ 3) 


I —_ 
p» nao cao MIEL(E + m + n)!T(m + 3) 


where we have used 
2" T(n + 3)n! = VY x(2n)! 
We next make the substitution 


1 1 e dt 
= : esos, (A3-2 
(€+m+n)! 2nide {htm 
where the path of integration C is a circle chosen large enough to ensure th. 
convergence of the series obtained when the order of summation and integra 
tion is changed. ‘The summations may now be performed: 


1 i +q2 - 2m m m 2 
a = | deel See eee 
al J 


m=) 


» pale 1/2 (A3-3 
1 | t mf: ai Cc) = gttatit iy 


3 faw¢=§ 


C encloses the pole at c + 6 and the branch point at c as well as the origin. 
When a? is zero the integral may be reduced still further. Let ¢ be com- 
plex and 6 such that the point c + 6° does not lie on the line joining 0 to: 
Deform C until it consists of an isolated loop about c + 6 and a loop about 
0 and ¢, the latter consisting of small circles about 0 and ¢ joined by two 
straight portions running along the line joining 0 toc. The contributions 
of the small circles about 0 and ¢ vanish in the limit. Along the portion 
starting at 0 and running to ¢, arg (t — c) = —a + argc, and along the por- 
tion starting at c and running to 0, arg (¢ — c) = w + argc. On both 
portions arg / = argc. Bearing this in mind and setting ¢ = c sin? @ on the 
two portions gives 
esin?6 


> 9 
cos e 


ano = wb(c + BY) 7 iy ey @  (AB4 
ee ie a il = 0 Tivaie \ 


The integral may be expressed in terms of the function 


Te(k, x) = fg e “Io(ku) du 
0 
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by noting that 


ai e a-—Bcost a7 1 l 
| ee [le tatmen a 
Jo a+ Bcosv 0 a+ Bcosv Jo 


l — 
2 2 —1/2 x A3-5 
= g(a — Bp) — 7/ e “ Io(Bt) dt 
0 
= wa? — B*)-!'? — (r/a)Te(B/a, a) 
Thus 
c/2 / F 9 b2+¢ Cc 
| ae = We Io(c/2) _— (wb aye Te (, P «) ( \ 3-6) 
2a 
where 
a b 7 ¢ Z \3 7 
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Noise in Resistances and Electron Streams 
By J. R. PIERCE 


5 in ‘HNICALLY correct results in a field are achieved initially in diverse 
and often confusing and complicated ways. Sometimes, such results 
are later brought together to give them a more unified form and a sounder 
basis; such critical summary and exposition is of great value. In quite 
another way, a worker who uses results established in a field will discover 
many plausible reasons for believing the results, and he will find eventually 
that an air of inevitability and “understanding” pervades the subject 
Such “understanding” is not to be confused with the process of rigorous 
proof carried out step by step, but it can help in organizing and making use 
of a body of related material. 

The field of ‘“‘noise’’, especially as it affects electron devices and communi 
cations in general, is one particularly troublesome to engineers. The sound 
work on the subject has commonly involved mathematics and especiall) 
statistical ideas unfamiliar to many who must deal! with the practical prob 
lems of noise. In early papers on noise, a great deal of heat was generated 
in acrimonious controversy between two schools, one of which assigned a 
uniform noise spectrum to certain noise sources, while the other held this to 
be inadmissible and got identical answers by more recondite means. 
Happily, a recent paper by S. O. Rice! clearly presents both approaches. 
Rice’s paper further provides a fine broad summary of noise problems to- 
gether with considerable original material. It does not extend far into the 
field of electronics.’ 

The reader who has sufficient time could achieve a profound ‘“‘understand- 
ing’’ of the circuit aspects of noise by reading Rice’s paper. The under- 
standing would involve familiarity with much mathematics useful in itself. 
To many engineers, however, this might prove a lengthy and painful process 

The writer proposes to present here a series of plausible arguments for 
believing certain facts about noise. Both simple circuit considerations and 
“electronic” effects (as, space charge reduction of noise) are included. The 
arguments presented are not intended to be original and it is not claimed 
that they are rigorous; they do seem to be easily understood, and to help in 
remembering and in using some important practical material. Starting 
points of the arguments, or “postulates”, have been chosen on the basis oi 
familiarity, not simplicity. No effort is made to point out all of the hidden 
assumptions in the arguments, but a few important ones are indicated. 
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An initial warning should be made that quantum effects treated in 
Nyquist’s original paper on Johnson noise, but afterwards much neglected, 
are entirely disregarded here. 


[. JoHNSON NotsE® 


In 1926, in an investigation of amplifiers with exceedingly high grid 
resistances, J. B. Johnson discovered that a resistance acts as a noise genera- 
tor having an open-circuit voltage with a mean square value 


v? = 4kTRB. (1) 
Here and subsequently, lower case letters 2 and i will be used in referring to 


noise voltages and currents. In (1), v? is the mean square value of noise 
voltage components of frequency lying in a small bandwidth B (sometimes 
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Fig. 1—Relations between noise power, noise voltage and noise current can be derived 
by assuming the noise source to be a voltage in series with a resistance 


called df or Af), k is Boltzman’s constant, and R is resistance. We easily 
see from Fig. 1 that the maximum noise power which can be made to flow 
from the resistance into a load (that which will flow into a matched load) is 


a ae (7) 
P= ~~ RTB. 2 
This ‘‘available noise power” is a convenient alternative formulation. 

If an impedance has a reactive as well as a resistive component, the open 
circuit noise is given by (1) where R is the resistive component; if an admit- 
tance has a conductance G the noise may be represented as an impressed 
current (that which flows when the admittance is short circuited) of 
magnitude 


12 = 4kTGB. (3) 


We see from (1) that if two resistances are connected in series, the total 
squared noise voltage is the sum of the squares of the noise voltages produced 
by the resistances separately, and from (3) we see that the noise currents of 
conductances connected in shunt also add by summing squares. This rule 
of addition holds for adding the noise of all independent sources. Of course, 
if noise from the same noise source reaches a point by different paths, the 
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voltage or current components near any frequency should be added dire: 
with due regard for phase. 

Johnson noise is related to many physically similar phenomena such as 
Brownian motion and the random fluctuations in position observed in te 
coils of very sensitive galvanometers. 

The simplest derivation of (1), (2) or (3) is that given by Nyquist‘ in a 
companion paper to Johnson’s. Consider a long lossless transmission line 
of length / terminated at each end in resistances equal to its characteris‘ 
impedance. Imagine line and terminations in thermal equilibrium at a 
temperature 7, as shown in Fig. 2. If electrical energy flows from the 
resistance at 1 to that at 2, then equal energy must then flow from 2 to 1, 
as any net gain or loss of energy would violate the second law of therm 
dynamics. 

Now, suppose that we suddenly close the switches at 1 and 2, short circuit 
ing the ends of the line. The line now becomes a resonator, having resonant 





f } 


Fig. 2—Two resistances terminating a transmission line act as generators of therma| 
noise power traveling along the line. 





rrPw- 


frequencies such that the line is ” half wavelengths long. The resonant 
frequencies will be 
f = n(c/2L). (4 


Here » is an integer and c is the velocity of light. The frequencies are 
separated by frequency intervals 


Af = (¢/2L). (5 


The energy which originally flowed right to left and left to right between th 
resistances 1s now reflected at the ends. It may be expressed as the thermal 
energy associated with the resonant modes of the line. According to 
statistical mechanics, there is an energy kT associated with each resonant 
mode. ‘The energy per unit bandwidth is obtained by dividing this by the 
frequency interval between modes, given by (5) and is 


w= kT/Af = kT/(c/2L). (6) 


Since it takes a wave a time L/c to pass completely through the line, this 
energy w represents the energy per unit bandwidth which flowed into the 
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ne from both resistances over a period Lc. If p is the power per unit 
bandwidth from one resistance, then 


2p(L/c) = w = kT/(c/2L) 
p = kT. 


Or, we may say that the power flow from a resistance into a matched load 
the available power) is, for a bandwidth B 


P = kTB. (8) 


Sometimes it may be desired to know the mean squared fluctuation voltage 
integrated over all frequencies. Carrying out such an integration for the 
voltage between a pair of terminals connected by a complicated network 
would seem to be a difficult procedure. However, if the pair of terminals is 
shunted by a capacitance, the integrated fluctuation voltage can be ob- 
tained by direct application of the principles of statistical mechanics. 

In a lumped network composed of capacitive, inductive and resistive 
elements each capacitance and each inductance constitutes a degree of 
freedom; that is, the electrical state of the network can be specified com- 
pletely by specifying the voltage across each capacitance and the current 
in each inductance**. According to statistical mechanics, the average 
stored energy per degree of freedom is k7/2. The stored energy in a 
capacitance is Cv?/ 2. Thus, the mean squared noise voltage of all frequen 
cies across a capacitance C must be 


v- kT _, (9) 


II 


Similarly, the mean squared noise current of all frequencies flowing in an 
inductance L is 


i? = kT/L. (10) 


We have conveniently thought of Johnson noise as generated in the 
resistances in a network. We need not change this concept and say that 
the voltage and current of (9) and (10) are generated in the capacitance or 
inductance any more than we would say that the thermal velocities of 
molecules are generated by the molecules’ mass. Relations (9) and (10) 
merely represent necessary consequences of the laws of statistical mechanics 
as, indeed, does (1). 

It is of some interest to illustrate the use of (9) and its connection with (1) 

* Strictly, such a lumped network is an unrealizable ideal. There are no pure capaci 
tances, inductances, or resistances. The conditions under which actual condensers, coils 
and resistors can be represented satisfactorily by these idealizations must be judged by 
measurement or calculation or by past experience or intuition 

** In enumerating the degrees of freedom, capacitances in series or shunt are lumped 
together as one element; the same holds true for inductances. 
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by a very simple example. Consider a resonant circuit consisting of 
capacitance C, an inductance L and a resistance Ro, all in parallel. The 
resistive component of the impedance across this circuit is 





—— . 
Ww Ww 

Q = RowC = Ro/woL (12 

wo = 1/V/LC. (13 


Here wo is the resonant frequency of the circuit and Q has its usual meaning 

From (1) we see that as Ro , the resistance at resonance (w = wo) is mack 
higher, the noise voltage for frequencies near resonance increases. How 
ever, if we regard wo and C in (12) as fixed, we see that as Ro is increase« 
the ( of the circuit is increased, the frequency range over which R is high is 
decreased, and R actually becomes lower far from resonance. (9) tells us 
that the mean square noise voltage integrated over all frequencies remains 
constant as Ro is changed. 

It is found that for a high Q circuit, the noise is much like a carrier oi 
frequency wo modulated by low-frequency noise. If we let the radian fre 
quency of this ‘noise modulation” be (w — wo), then the mean square 
amplituce of the noise modulation varies with frequency about as R given 
by (11) varies with (w — wo). 


II. ScrootrKy Norse ork SHot NOISE 


In 1918 Schottky® described the ‘‘Schrot-Effekt”: the noise in vacuum 
tubes due to the corpuscular nature of the electron convection current 
This is commonly known as “shot noise.’”” The magnitude of this noise is 
usually derived by means quite different from those used here. 

Johnson noise is necessarily associated with any electrical resistance, 
whatever its nature. Now, consider a close spaced planar diode shown in 
Fig. 3 consisting of two opposed emitting cathodes, each emitting a current 
Io. Suppose the whole diode is held at the same temperature. There are 
no batteries or other sources of power aside from thermal energy; the only 
electrical energy flow must then be Johnson noise, ascribable to the re- 
sistance of the diode. 

Assume that the cathodes both have the same uniform work function 
Then when the diode is short circuited, each electron emitted from cathode | 
will reach cathode 2, and each electron emitted from cathode 2 will reach 
cathode 1.* If cathode 2 were made negative, all the electrons from 2 would 


* It is here assumed that J» is small enough so that depression of potential due to spac« 
charge is avoided. 
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continue to reach 1, but some of the low-velocity electrons leaving 1 would 
be turned back from 2. 

It is well known® that if a Maxwellian velocity distribution is assumed 
for the electrons leaving 1, the electrons which can overcome the retarding 
field and reach 2 are found to constitute a current 


juke”, (14) 


Here J is the total current carried by electrons leaving 1 and V is the voltage 
of 2 with respect to 1, which has been assumed to be negative. 


. lo — 
a 
i} 


N 














v= 
L 











Fig. 3—An electronic resistance formed by two opposed cathodes at the same tem 
perature acts as a generator of thermal noise. 


By differentiating (14) we can obtain the diode conductance G at V = 0, 
and we find 
e 
G = Ip. (15) 
kT 
From (3) when the diode is short circuited and the voltage is zero we have 
a mean square noise current 


2 = 4kTGB = ae (4kTB) = 4ely B. (16) 


This noise is the sum of the noise due to two independent noise sources 
(the noise in the two currents J»). That due to either current Jo is* 
z= 2eT B. ( 17) 


*In this section, we are concerned with short transit angles only and no distinction 
need be made between the current induced in the circuit, i, and the electron convection 
current. 
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This is the expression for shot noise in a randomly emitted current, as 
temperature limited emission or in photo electric emission. 


Ill. Noise OTHER THAN SHOT NOISE: ELECTRON MULTIPLIERS 


Let us consider a class of systems in which the average output current is 
proportional to the average input current, in which an electron of charge, 
entering produces an output charge, me instantaneously, and in which the 
probability that any electron will produce m electrons is p,, . 

If the input current is 7), the average output current is 


[ nl (18 
n= > npn. (19 


It is easy to persuade ourselves that any frequency component of current, 
noise or signal, will produce an output current # times as great; this happens 
to be true,gand we will use the fact. 

Let us consider our device when it has randomly emitted electrons as an 
input. At the output we will see appear groups of 1, 2, 3 etc. electrons, 
each group caused by the entrance of a single electron. If Jo is the total 
input current, the output current consisting of groups of » electrons is 


T, = nlopn. (20) 

Each group carries a charge we. We may now use (17) to write the noise 

in the part of the current carried by groups of electrons, replacing the 
electronic charge, e, by the group charge, ne 

2% = l2(ne)(nIo p,)B. (21 


As there is no correlation between entering electrons, the total mean square 
output noise current will be the sum of the noise components carried by 
groups consisting of different numbers » of electrons. Summing (19) with 
respect to m we obtain 


1} = 2elo Bn* (22 
> 9 
n> = > n Py. (23 
n 


Now, the input has been taken as having shot noise. A part of the noise 
output is to be attributed to this input shot noise amplified by the device; 
that is, it will be # times the input shot noise. 


2 = n° 2elyB. (24 


The part of the noise output current due to the fact that an electron does 
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t produce # electrons, but may produce 0, 1, 2--- etc. electrons, must 
he the difference between (22) and (24), or 


Nm 
cs 


i? = ely Bin? — nv). 


[he quantity in parentheses is the mean square deviation in . 
If the input current has any noise components 7%,, then the total noise 
output component will be 


- = tj + hh 77. (26 


By applying (26) successively to stage after stage the noise output of a 
multistage electron multiplier can be eva'uated (if one knows (n? — n2)).? 

Wonder is sometimes expressed that current can be noisier than shot 
noise, in which the time of electron arrival is purely random. Obviously, 
we can have more than shot noise only if there is something non-random 
about the time of electron arrival, and the argument above discloses just 
what this is; it is the arrival of electrons in bunches.** We can easily see 
how erratic even large currents would be if electrons were bound together 
in groups having a total group charge of a coulomb, all the electrons in a 
group arriving simultaneously. Reverting to our shot noise formulas, 
we may illustrate this by assuming a perfect multiplier with a shot noise 
input, in which each input electron produces exactly V output electrons. 
Arguing from the shot noise equation (17) and replacing e by .Ve we should 
expect an output noise current 


i? = 2(Ne)I,B (27) 
where /; is the output current; we get exactly the same result by assuming 
the input noise current squared amplified by V? 

i? = (2eloB)N? 
= 2(Ne) (NIp)B (28) 
= 2(Ne) 1,B. 


*The mean square deviation js the sum with respect to of the square of the devia- 
tion from the mean value of n, 7. 


= (m — 7)? p, = Dn? p, — WNIT np, + NZ fn 


lhe summation in the first term is m?, that in the second term is % and that in the third 
term is unity. Hence 


= (nm — 7i)* pa = (n? nz 
** Anything, (such as transit time difference for elec ithin a | \ I 
Anything, (such as transit time difference for electrons within a bunch) which tends 


to break up the bunches will reduce the noise—and the signal as well. Such noise reduc 
tion involves a return to a more nearly random flow. 
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Conversely, we are led to wonder whether a current less bunched thar 
that produced by random emission might not have less noise. The most 
smoothly distributed current we can imagine is that of fo electrons per 
second emitted at evenly spaced intervals. Obviously, such a current wi! 
have a spectrum consisting of frequencies nfo, integral multiples of / 
Thus for f < fo, there will be no “noise” and similarly for fo< f< 2/, 
2fo< f< 3fo, etc. 

For a current of 10 ma, fo = 6.3 X 10"; thus, even for small currents ar 
evenly spaced emission would have no a-c components in the radio-fre 
quency range; this is a comforting thought in considering space-charg: 
reduction of noise, which is discussed in section 5. However, purely to 
satisfy our curiosity we may pursue the matter a little further. If we 
assume that each electron constitutes an instantaneous pulse of current, a 
simple harmonic analysis shows that the a-c current component of fre 
quency fo will have a mean square value 


i= = 2elofo. (29 


Thus, in each interval fo wide centered about a frequency 1 fo there will be a 
mean squared a-c current equal to that which would be associated with the 
same band for random emission with the same current. By making the 
emission regular we have not reduced the mean square “‘noise”’ current in a 
broacl frequency range; we have merely changed its frequency distribution 
from a uniform distribution to a distribution of sharp, high peaks. 


IV. PARTITION NOISE 


Consider a tetrode, shown in Fig. 4, with a cathode current J, , a screen 
current J, , and a plate current J,. 

The grid current is taken as zero. Suppose that the screen is very tine, 
so that every electron leaving the cathode has the same chance of striking 
the screen, regardless of its point of departure. We may now regard the 
function of the screen as that of a peculiarly simple electron multiplier, 
for which m can be zero (electron striking screen) or 1 (electron passing 
screen). 

The probability of an electron passing the screen is /,/Z.. According 
from (19) and (23), 


ly, 


n= I,/Te (30 
n? = I,/I.. (31 
Suppose we write the noise in the cathode current as 


2 = T2e/.B 32 


oe 


Here I, a factor less than unity, is introduced to account for the “spac¢ 
charge noise reduction” in space charge limited flow. 
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Now, by applying (25) and (26) we obtain for the noise in the plate 


current 
i, = 2el.B(1p/Ie — (I,/I-)*) + T*2el. BUIp/I.)* 
ie (33) 
i, = 2el, Bi — (1 — I)(/,/I.)). 
It is to be noted that if ©? = 1, that is, if the cathode current is random, 


the noise in the plate current is purely shot noise. The screen cannot make 
the plate current noisier than shot noise since it does not act to produce 
bunches of electrons. 

The noise in the screen current can be obtained by substituting J, for J, 
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Fig. 4—Electrons randomly hitting or missing the screen grid make a tetrode noisier 
than a triode. 


in (33). There is a correlation between the screen and plate noise currents; 
the total noise in the screen current plus the plate current must, of course, be 


THR = B= 2ehB (34) 


and not the sum of #2, and 22. 
Partition noise has been discussed by Thompson, North and Harris.® 


V. SpaAcE CHARGE REDUCTION OF NOISE 


In this section an approximate derivation of noise in a space charge lim- 
ited diode will be presented. The derivation leads to an expression valid 
for many practical tubes and illustrates the nature of the noise in space 
charge limited flow. 

Consider a parallel plane diode of unit area and spacing x, with an applied 
voltage Vo, as shown in Fig. 5. When the voltage is applied, the electron 
convection current in the diode rises to value J). Neglecting thermal 
velocities of electron emission, this current is such that the electronic 
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‘space chatge”’ associated with it causes the voltage gradient at the cat) od 


surface to be zero. A greater current would mean a negative gradie: 
the cathode and hence no emission; a smaller current would mean a pos (ivy 
gradient at the cathode and unlimited emission. On this basis Child’s 

is derived, which gives the current per unit area /9 in amperes in terms oj 
the voltage Vo and the spacing in centimeters x as 
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Fig. 5—Part of the electrons leaving the cathode of a diode are turned back befor 
reaching the potential minimum; others proceed to the anode. Ordinarily the greater 
amount of noise is associated with the space between the potential minimum and the anode 


From (35) we can obtain a useful relation for the conductance G 
G = 10/8Vo = (3/2)(10/Vo). (36 
The resistance R is 
R = 1/G = (2/3)(Vo/Io). (37 


In actual diodes, the electrons are emitted from the cathode with a 
thermal velocity distribution; a potential minimum of some negative voltage 
Vm is formed at some distance x, from the cathode surface. If the magni- 
tude of the emitted electron current is 7, and the actual current passing the 
potential minimum is J , then because of the Maxwellian velocity distribu 
tion we have 

Io ad | OP sia al 


11,600 »_/T 
= [,e aT 

















before 
rreater 
anode 


(36 


(37 


ith a 
tage 











nt 8 RON! ARO NRE nS em 


sea NeRnA AiR sh een 


Le hen 2 mee ND ea a OO NEI EP ee TNE thet A Hal 





NOISE IN RESISTANCES 169 


Ordinarily, the magnitude of V,, is very small compared with Vo; Jo is 
very small compared with /, and x,, is very small compared with x. 

Suppose V,, were held constant, say, by putting a conducting plane of 
potential V,, at x, . Then, tie electrons which pass this plane are quite 
independent of the low energy electrons which are turned back, and hence 
in the current passing x,, there will be pure shot noise. 


1'= Jel ob. (39) 

Now suppose we change V,,. Th» change in Jo will be, from (38), 
dIg = dVn/ Rm 40) 
Ry = (elo/ kT)“. +1) 


If we use a constant current instead of a constant voltage d-c supply, then 
V,, must fluctuate in such a way as to cause a current equal and opposite to 
(39), or, there must be a fluctuating voltage v3, such that 


v, = 2eloBRi, 
= (1/2)4kTR,,B. 


42 } 


Suppose we consider the noise fluctuation of the anode voltage of a space 
charge limited diode supplied from a constant-current source. If there 
were no fluctuations in the voltage drop between the potential minimum 
at % and the anode at x, (42) would give the noise voltage fluctuation of 
such an “open circuited” diode. Actually, much larger fluctuation voltages 
are observed, and we must conclude that they arise in the space between the 
potential minimum and the anode. As the current is constant in this region 
(by definition—we have assumed a constant-current supply) we are forced 
to conclude that such fluctuations are due to a variation of mean electron 
speed in this region. The field at x,, is necessarily zero. Hf, with a constant 
current, electrons travel more rapidly between x,, and the anode, there is 
less electronic charge everywhere in this region, the rate of change of field 
with distance, and hence, the field, are everywhere smaller, and the voltage 
between x,, and the anode at x will be smaller. 

It is somewhat involved to treat the problem of mult!-velocity flow 
exactly; this has been done by Rack* and others*!°"'; however, Rack has 
shown that an approximate treatment yields very nearly the correct result 
over a fairly wide range of conditions. In this approximation, the stream 
of electrons with many velocities and a fluctuating mean velocity is replaced 
by a stream in which all electrons have the same velocity, and this has a 
mean square fluctuation equal to that of the multi-velocity stream. 

Let us now measure x from the potential minimum. Suppose we con- 
sider an electron which passed the potential minimum (x = 0) at ¢ = 0. 











170 BELL SYSTEM TECHNICAL JOURNAL 


he field at the potential minimum is zero. The charge which has flowe 
in behind the electron at the time tis —¢Jo. Hence, from Gauss’s theore1 
the potential gradient is 


AV /dx = Igt/e (43 
where e is the dielectric constant of vacuum. We have for the acceleratio 
e lot 
ie anee (44 
m € 
Ii at the time ¢ = O (at the potential minimum), x = 0, = % 
, ellos, . nee 
4=——b + XK (45 
m 2e 
elo. . 
x= —t + Kl. (46 
m 6€ 


Now the voltage V between the potential minimum and any point x must 
be such that 


#—-8=25Y, (47) 
m 
1 sar > Io 2. 
AL be - e oes 
: + (é%) at adi (48 
“Mm 


At any fixed point x, if we vary %o by a small amount d%p, we find by dif 
ferentiating (46) 


_ ee ae 
dx (: Io P + in) (49) 
m 2e 
From (48) 
dV» -Mt(fhey in) dr + 2 Pate (50) 
e \m2e 2e 
Using (49) 
: 2k - 
dVo = —— TF da. (51) 
2e 


It now remains to evaluate ¢. For most cases, the thermal velocities at 
the potential minimum are so small compared with the velocities in most 
of the region between the minimum and the anode that we can take the 
value of ¢ for #9 = 0. Then, from (45) and (47) 


. e Io coat | — TP] 
S = a 2— Ve ‘ (52) 
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From (51) and (52) 


‘ 1/2 
dV, = ~2"(¢) Vo" dio. 53 


m 


Now, if (dz)? is the mean square fluctuation in velocity, the mean square 
fluctuation in voltage will be 


Ber 73: \¢ 
o? = 2(e/m) Vo (dx,)?. 54 
The assumptions leading to (54) are those leading to Child’s law, and thus 
we can use (37) in connection with (54), giving 
vp? = 3(e/m) “IoR (dx)? 55 
It now remains to evaluate (d%,)?, the mean square fluctuation in the 
velocity of the electrons passing the potential minimum; to do this, we 


return to (25). Suppose N is the number of input electrons per second 
The output current can then be written 


and we can call the fluctuation in it 


7? = (5nNe)*. 


on 


Equation (25) applies for no fluctuation in J» and hence for no fluctuation 
in N; e is a constant, and thus we may write (25) as 


(8a)? = — (nm? — a). (58 


We may generalize this to say that each electron has a probability p of 
producing some effect of magnitude m and the fluctuation in the magnitude 
of the effect is (67)?. Before, we said that an electron had a probability 
p of producing m secondaries. Now we will say instead that an electror 
has an uncorrelated probability p of having a velocity «, and obtain for the 


mean fluctuation in the velocity, (dx)? 
——._ 2B 
(dX)? = = (142 som tf 59 
N 
In a Maxwellian distribution, the number of electrons passing a plane 


perpendicular to the direction of motion per second having velocities lying 
in the range du at « is 


“ 3 
dn = Aue ™* '™7® du. (60 
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Here 7. is cathode temperature. We see @ and u? are 


@ 
2 mu?/2kT ¢) 
| ueé © du 
0 


i= “= = ] 
nu? /2kT, 
i ie ‘“ du 
0 


mai 


/ kT. 
4 / — (6 


u? mu?/2kT. du 
€ 
: 40 eT. 
u? = = 2 : (6 
7 m ery 
mu? /2kT. 
ile du 
/0 
Accordingly 
ue— aw = 3(4 — ow) (63 
m 
Combining (63) with (59) we obtain 
B RE 
(dx)? = = (4 — mr) , (64 
N m 


Combining (62) with (53) and remembering that Jy) = Ne we find the mean 
square open circuit noise voltage to be 


vy = 3(4 — w) kT.RB 
(.644) 4&7.RB. 


(65) 


This is the chief contribution to noise in a space charge limited diode. 
Usually R is substantially equal to the plate resistance of the diode (11 

does not include effects on the cathode side of the potential minimum 

Hereafter R will be treated as the total plate resistance of the diode. 


VI. Noise iN TRIODES AND PENTODES 


Consider the triode shown in Fig. 6. Here we have a cathode, a grid, 
and a plate. The input admittance of the tube is represented in the di 
agram by the grid-cathode capacitance C; and the grid-plate capacitance C2 
The resistance R, is a fictitious noise resistance which will be evaluated 
later. It is assumed to act between the input admittance of the tube and 
the controlling action of the grid; no current can flow in R,, because the grid 
as indicated in the diagram is presumed to present an open circuit. 

We will regard the cathode-grid region of the triode as an ‘‘equivalent 
diode.” The anode voltage of the diode is taken as 


Vo = (V, + V>/u). (66 
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NOISE IN RESISTANCES 
Here V, is the grid voltage and V, the plate voltage of the triode. 
ate voltage is held constant and yu is taken as constant 
aV_= dV,. 
Hence, under these conditions 
01o/8Vo = G = O1p/AV,. 


Here G is the conductance of the equivalent diode, the reciprocal 


which appears in (65), and is also the transconductance of the triode. 
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(67) 


(603) 


of R 


As we wish to calculate the noise with no a-c grid or plate voltage, and 


as these through (64) specify the plate voltage Vo of the equivalent diode, 
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C; GRID-CATHODE CAPACITANCE 
Co GRIDO-PLATE CAPACITANCE 


Fig. 6—Low-frequency noise in a triode can be ascribed to a fictitious noise resistance 


R, , acting into an open circuit to cause voltage fluctuations on the grid 


the equivalent diode may be regarded as short-circuited. Hence, the 


current will be 


= (.644) 4kT, GB. 
If we express this as shot noise reduced by a factor I’? we obtain 


2 = el, I’ B 


II 


kT .G 
(.644) 2kT ‘ 
ely 


Il 


noise 


(69) 


(70) 


Often, the noise expressed by (69) is ascribed as a fictitious noise resistance 


R, , at room Temperature 7, connected between the grid-cathode « 
tance and the ‘controlling action’ of the grid as shown in Fig. 6. 


apaci- 


This 


fictitious resistance looks into a complete open circuit; hence, it has a noise 


voltage 


vo? = 4kTR,RB 


(71) 
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and produces a noise plate current (for zero load resistance) 


i? = 4kTR,BG?. 


~ 


tm 


Comparing (69) with (72) we find 
R, = (.644/G) (T./To). 73 


Here T> is a reference temperature, usually taken as 290° K. The effi 
load impedance on signal from this fictitious resistance is treated by purely 
circuit means. 

In pentodes there is noise according to (69) and in addition there is par 
tition noise according to (33). By taking I from (70) and equating th 
noise current given by (33) to (72) the fictitious “noise resistance” of » 
pentode can be evaluated in terms of g, /,/I. and T./T. 
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Abstracts of Technical Articles by Bell System Authors 


Gross-linkage of Linear Polyesters by Free Radicals... W.O. BAKER. Re- 


| actions fundamental to the use of the new low pressure laminating or casting 


fer 





resins have been studied. The striking property of these plastics, which are 
usually based on some polyester and a vinyl monomer, is their rapid and 
easy curing, leading to unique ease of fabrication. This curing, the forma- 
tion of a permanent three-dimensional polymer network, or gel, is achieved 
by reaction with a source of free radicals, such as from an organic peroxide. 
These agents cause polymerization of the vinyl monomer, as was previously 
understood, but they also seem to incorporate the polyester into the network, 
even if the polyester contains little or no unsaturation. 

Investigation of a series of simple polyesters, the polyundecanoates, 
showed that free radicals, such as come from the decomposition of benzoy] 
peroxide, could cross-link or gel the linear, saturated, chains. Apparently 
the hydrogen atoms in methylene groups next to polar groups like the car- 
bony], i.e., the a-hydrogens, are removed by the free radicals. The resulting 
chain radical attacks an adjacent chain, and a cross-link is formed. The 
effects of cross-links thus produced on solubility, dilute solution viscosity, 
melt viscosity and, finally, stress relaxation of the cured solid were examined. 
Probably the similar activity of a-hydrogen atoms is important in the chemi- 
cal aging or weathering of plastics and rubbers. It is likewise significant for 
the vulcanization of many synthetic rubbers. 

Rubberlike Products from Linear Polyesters” B.S. Biccs, R. H. Erick- 
son and C. S. FuLLER. The polymers which result from the condensation 
of dibasic acids with propylene glycol are viscous gums which can be vulcan- 
ized to rubberlike products. In the unpigmented condition these rubbers are 
quite weak, but when reinforced with suitable pigments their strengta and 
elongation compare favorably with other synthetic rubbers. Because poly- 
esters of known structure and molecular weight can be easily synthesized, 
these polymers are useful for the study of the relations between structure and 
properties in rubberlike materials in general. Factors affecting tensile 
strength, oil resistance, brittle temperature, and stability are discussed 

Pulse Code Modulation’ H.S. Brack and J.O. Epson. A radically new 
modulation technique for multichannel telephony has been developed which 
involves the conversion of speech waves into coded pulses. This new tech- 

' Jour. Amer. Chemical Soc., May 1947. 


? Indus. & Engg. Chem., September 1947. 
> Telephony, August 30, 1947. 
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nique is called Pulse Code Modulation or simply PCM. An eight-cha ine 
system embodying these principles was developed and produced in port bi 
form for field operation. Other work carried on simultaneously by W. M 
Goodall (see B.S. T. J., July 1947) resulted in the development of an ex) eri. 
mental system using a different method of coding. 

In carrying out this new type of modulation, the speech wave applied to 
each channel is, in effect, transmitted sample by sample, and each sample is 
represented by a multi-unit code employing on-or-off pulses, hence the term 
PCM. 

This method appears to have exceptional possibilities from the standpoint 
of freedom from interference. Its full significance in connection with future 
radio and wire transmission may take some time to reveal. 


Stereoscopic Drawings of Crystal Structures.4 W.L. Bond. A method is 
presented for getting stereoscopic pairs of atomic structure views given the 
coordinates of the atoms and cell constants. 

Properties of Liquids at High Sound Pressure.” H. B. BricGs, J. B. Jouy- 
son and W. P. Mason. When sound of high amplitude is transmitted into 
a liquid by means of a mechanical driving device, the ultimate limitation to 
the power that can be transferred is cavitation or breakdown of the liquid 
under high internal stresses. A study of cavitation has resulted in establish- 
ing the following results. Under steady-state conditions, light liquids filled 
with air cavitate when the negative acoustic pressure reaches the atmos- 
pheric pressure. When liquids are degassed, their natural cohesive pressure 
becomes effective and they will withstand a negative acoustic pressure. It 
is found that the total negative pressure required to cause cavitation is equal 
to the sum of the cohesive pressure—tensile strength—and the ambient 
pressure. Viscous liquids have a higher cohesive pressure and a proportion- 
ality has been established between the logarithm of the viscosity and the 
cohesive pressure. The amount of power that a liquid can withstand 
increases markedly as the pulse length is shortened. 

An explanation of these phenomena is attempted on the basis of Eyring’s 
theory of viscosity, plasticity and diffusion. On this theory natural holes 
exist in the liquid into which molecules can jump, leaving holes behind them 
A jump occurs when the molecule has accumulated enough heat energy to 
surmount an activation potential barrier of energy value Eo. Cavitation 
appears to be the result of coalescing of the natural holes in the negative 
pressure phase of the cycle. Since a molecule has to jump from a hole in 
order that this can coalesce with another hole, the cavitation pressure is pro- 
portional to the activation energy which in turn is proportional to the log 


* The American Mineralogist, July-August 1947. 
5 Jour. Acous. Soc. Amer., July 1947. 
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irithm of the viscosity. The increased power-transmitting capacity for 
short pulse lengths isa result of the finite time taken for the small holes to 
grow in size to a large enough hole to cause rupture of the liquid 

Modulation in Gommunication® F. A. Cowan. The fundamentals in- 
volved in introducing signals into one medium and transmitting them 
through another are simplified in this review, so that the relationships be- 
tween the many varieties of modulations attempted or in contemporary use 
are formed into a cohesive whole. 

lir-borne Magnetometers.’ E. P. Fetcu,* W. J. Means,* T. SLonezew- 
ski,* L. G. PARRATT, L. H. RumBAvuGH and A. J. TICKNER.* Developed 
under the impetus of the submarine menace of World War II, the air-borne 
magnetometer has found many peacetime uses. Navy airplanes equipped 
with magnetometers for exploration of Antarctica were used in the recent 
United States Navy expedition. An expedition now is studying the Aleutian 
Alaskan volcanos and the Aleutian submarine trench. From there it will 
proceed to Hawaii and Bikini. 

The Generation of Centimeter Waves.* H.D. Hacstrum. The electronic 
devices used most extensively, recently, for the generation of centimeter 
waves are discussed. The physical form, operating capabilities, and the 
basic physical principles of operation of the triode, velocity-variation, and 
magnetron oscillators are presented. An attempt is made to show how these 
oscillators are related to one another. For a variety of reasons, particular 
emphasis is placed on the magnetron oscillator. 

Selective Demodulation.2 DoNALp B. Harris. A method of demodula- 
tion is proposed in which the output current of the demodulator is a linear 
function of the input voltage, while at the same time provision is made for 
producing the necessary product terms which will result in demodulation. 
Demodulation is brought about by integrating the product of the instan- 
taneous value of the modulated wave by the instantaneous value of a wave 
having the same frequency and phase as the carrier. Where this method of 
demodulation is used it is proposed that two carriers in quadrature on the 
same frequency may be employed, reducing the bandwidth to that re 
quired for single-sideband transmission. 

It is suggested that the required linear demodulation characteristics may 
be obtained through the use of “electron-coupled” demodulators. Theo- 
retical considerations indicate that, when demodulation of this type is em 
ployed, selectivity ahead of the demodulator may be dispensed with, the 

° Elec. Engg., September 1947 

? Elec. Engg., July 1947 

* Of the Bell System. 


* Proc. I. R. E., June 1947 
* Proc. I. R. E., June 1947 
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signal-to-noise ratio is improved, greater economy of spectrum sp 
obtained, the number of tubes required is materially reduced throug) th, 
use of a common intermediate-frequency amplifier for a number of channels 
and any impairment due to the instability of the carrier or oscillator fre. 
quency is reduced. 

As an example of the possible application of the principles outlined, , 
hypothetical eight-channel transmission system is described. 

The Physical Significance of Birkhoff’s Gravitational Equations.’ px 
BERT E. Ives. Birkhoff’s gravitational equations are put in terms of di ip 
place of the local time ds used by him. The transformed equations show that 
Lorentzian mass has been used, and to the Newtonian attractive force js 
added a force normal to the direction of motion, v?/c? times the component 
of the gravitational force normal to the motion. 

Attenuation and Scattering of High-Frequency Sound Waves in Metals and 
Glasses." W. P. Mason and H. J. McSxkimin. By using a pulse method, 
attenuation and velocity measurements have been made for aluminum and 
glass rods in the frequency range from 2 to 15 megacycles. The sound 
pulses are generated by crystals waxed to the surface of the rod. This wax 
joint limits the band width of the transmitted pulse and measurements are 
made using long pulses which approach steady state conditions. The re- 
flected pulses show evidence of several normal modes which can be minimized 
by using specially shaped electrodes. Longitudinal waves show delayed 
pulses of smaller magnitude that are caused by the longitudinal wave break. 
ing up into reflected longitudinal and shear waves at the boundary. This 
effect is small if the diameter of the rod is Z0 wave-lengths or more. 

The measured losses for aluminum rods show a component proportiona 
to the frequency and another component proportional to the fourth power 
of the frequency. The first component is the hysteresis loss found for mosi 
solid materials. The component proportional to the fourth power of th 
frequency is caused by Rayleigh scattering losses which are the result o/ 
differences in the elastic constants between adjacent grains caused by) 
changes in orientation. Calculated scattering losses agree quite well wit! 
the measured values. The fourth-power scattering law holds quite wel 
until the grain size is equal to one-third of a wave-length. For higher fre- 
quencies the scattering loss increases more nearly with the square of the 
frequency. Glasses and fused quartz have a loss directly proportional to 
the frequency, showing that any irregularities must be of very small size. 

The Growth of Auditory Sensation.” W.A. Munson. The integration o! 
sensation with respect to time was studied experimentally by means of tones 

'0 Phys. Rev., August 1, 1947. 


" Jour. Acous. Soc. Amer., May 1947. 
12 Jour. Acous. Soc. Amer., July 1947. 
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of short duration. Loudness tests were made on sounds persisting from 
0.005 to 0.2 second and covering a wide range of levels. The observed in 
crease in magnitude of a sensation as the duration time is increased is 
attributed to the integration characteristic of the central nervous system, 
and an equivalent electrical circuit is derived. The circuit analogy is then 
used in the computation of loudness as a function of the duration of the 
stimulus. 

The Physics of Electronic Semiconductors.’ G. L. Pearson. The band 
theory of solids is capable of explaining such fundamental properties of elec- 
tronic semiconductors as the dependency of specific resistance on impurity 
content, the negative temperature coefficient of resistance, the sign of the 
Hall and thermoelectric effects, and the direction of rectification. Measure- 
ments of the specific resistance and the Hall constant enable the calculation 
of density, mobility, and mean free path of the electric carriers as a function 
of temperature and impurity. 

{utomatic Frequency Control of Microwave Oscillaters.4* VINCENT C 
Ripreout. A method for the automatic frequency control of any type of 
tunable microwave oscillator is described. In this method a servomechan- 
ism is used which includes a wave-guide discriminator circuit, a mercury- 
contact relay, a 60-cycle amplifier, and a small two-phase induction motor. 

Tests made on a preliminary model of a circuit of this type used with a 
4000-megacycle oscillator showed that a stability of one part in 50,000 was 
obtainable. The manner in which such a control system may be used in a 
microwave repeater is described. 

Proposed Method of Rating Microphones and Loudspeakers for Systems 
Use.™ FRANK F. RoMANOW and MELVILLE S. HAWLEY. Proposed, is a 
method of rating microphones and loudspeakers whereby the over-all per- 
formance of a sound system may be determined by adding together the 
microphone and loudspeaker ratings and the gain of the interconnecting 
network. This sum gives the performance quite accurately for most sys- 
tems. However, in some combinations of elements correction terms must 
be added. The formulas for these correction terms are derived. 

The proposed microphone and loudspeaker ratings have the additional 
usefulness of being in a form which permits the comparison of instruments 
of different impedances. 

Sulfur Linkage in Vulcanized Rubber.® Mutton L. SELKER and A. R. 
Kemp. The reaction of 2-methyl-2-butene with sulfur at 141.6°C. was 
studied. Reaction time and concentration paralleled those common in 

'8 Elec. Engg., July 1947. 

4 Proc. I. R. E., August 1947. 


' Proc. 1. R. E., Waves and Electrons Section, September 1947 
'6 Indus. & Engg. Chem., July 1947. 
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rubber-sulfur vulcanization. The results offer further insight into th: vy 
canization problem. The products of the reaction are liquids of the po  ysu). 
fide type R-S,-R, where x varies from 2 to 6 and R is an alkyl or alieny 
group and two solids (CsH¢S; and a higher homolog). The polysulfides ap 
pear to be somewhat richer in hydrogen than is expected from reaction oj 
two CsHio molecules, whereas the solids are hydrogen-poor. The structure 
of an acid anhydride in the sulfur system showing thione-thiol tautomerism 
is proposed for CsH¢S;, which is therefore 2,5-dithione-3-methyltet; ihy- 
drothiophene. The color changes with reaction time, from yellow to red to 
black, parallel those of rubber-sulfur vulcanizates. As in rubber-sulfur 
vulcanization the sulfur reaction rate is directly proportional to time, al- 
though the absolute rate is twice that in the polymer system. Starting with 
equal mole quantities of olefin and sulfur, there is a considerable amount of 
unreacted olefin in the system when all of the sulfur has reacted. The 
shorter the reaction time, the higher the value of x in the polysulfid 
R-S,-R and the larger the percentage of residues R that are saturated. 

On Hearing in Water vs. Hearing in Air.” L. J. Stvtan. The paper deals 
with the ability of a submerged listener to hear sounds generated in the air 
above him, compared with their audibility when his head projects above the 
water. Ina theoretical discussion it is shown that at 1000 c.p.s. a loss of 
the order of 45-55 db might be expected in the in-water audibility relative to 
the in-air value. This involves a number of assumptions, e.g., that there is 
no appreciable noise created by the listener’s propulsion, and that the effect 
of hydrostatic pressure unbalance on the eardrum is negligible. A few 
measurements made at 1000 c.p.s. and 3000 c.p.s. yielded values which are 
not at variance with the theoretical analysis. 

Cathode Phase Inverter Design.’ C.W.VADERSEN. Part I of this paper 
covers the general analysis of the cathode coupled phase inverter and de- 
velops formulae that enable the designer to compute circuit elements with 
good accuracy. The theory developed shows that degeneration exists onl) 
in the driven side of the amplifier and is limited to 6 decibels. Balance is 
discussed in terms of the tube parameters and external resistances, its being 
shown that considerable stability is attainable. A form of the inverter 1: 
which the power output stage utilizes a transformer with an unbalanced 
plate winding is presented. This is shown to give a true power balance ina 
manner analogous to the unbalanced plate resistor form of the voltage 
amplifying inverter. 

Part II presents a graph of the general design equations and illustrates its 
use with several working examples. 


17 Jour. Acous. Soc. Amer., May 1947. 
18 Audio Engineering, June and July 1947. 
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